找回密码
 注册
关于网站域名变更的通知
查看: 620|回复: 1
打印 上一主题 下一主题

FRFT程序运行时说Not enough input arguments. 是怎么回事???

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2019-8-2 14:14 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
如题:FRFT程序运行时说Not enough input arguments. 是怎么回事???
! C4 Z/ H6 k" S6 p2 b% S  X. Y' y, }1 c2 _, ~/ u9 M

function Faf = frft(f, a)
( s1 q- X# S/ l  K8 I) i1 f* Q, t4 g% The fast Fractional Fourier Transform- j, {7 z0 I. z3 g. S; i9 C
% input: f = samples of the signal
3 ?# H3 [  d& J% a = fractional power
# f* A! T! n) E( B7 S: i  {" D) e5 X% output: Faf = fast Fractional Fourier transform

error(nargchk(2, 2, nargin));

f = f(;
- \3 y. ~  }) t  J/ f3 }8 tN = length(f);  y, Y: V9 O0 o9 E% t
shft = rem((0:N-1)+fix(N/2),N)+1;+ k, g1 z: M$ S; M' H; m
sN = sqrt(N);
2 h* u# Q6 Y/ m/ Z/ n, F: Q4 r6 Da = mod(a,4);

% do special cases; C% R  ^' Z- C3 s* X
if (a==0), Faf = f; return; end;
/ C# o3 ]3 M5 N, y( wif (a==2), Faf = flipud(f); return; end;( Y% j0 H6 P% ^9 U% E. l
if (a==1), Faf(shft,1) = fft(f(shft))/sN; return; end . `% W  m2 X* |) v! B/ n
if (a==3), Faf(shft,1) = ifft(f(shft))*sN; return; end

% reduce to interval 0.5 < a < 1.5) C& S% Q4 F! l# \6 d6 d- ]$ R! E
if (a>2.0), a = a-2; f = flipud(f); end! _7 W- T$ K8 x  f) _- Y6 h
if (a>1.5), a = a-1; f(shft,1) = fft(f(shft))/sN; end
6 p6 B) B  X$ c2 r( Nif (a<0.5), a = a+1; f(shft,1) = ifft(f(shft))*sN; end

% the general case for 0.5 < a < 1.5
- c9 g- j( b3 @5 e( E# @/ balpha = a*pi/2;
9 R+ w. m: P( d' v" W: T5 |tana2 = tan(alpha/2);
/ H* L2 X3 @8 K" M7 K7 i0 }6 l; ksina = sin(alpha);8 |$ R& E. r* R) |$ Z
f = [zeros(N-1,1) ; interp(f) ; zeros(N-1,1)];

% chirp premultiplication
4 U/ z5 l, W# ~0 U& |- Ichrp = exp(-i*pi/N*tana2/4*(-2*N+2:2*N-2)'.^2);
3 B! v& w, }5 T$ X# p; S5 y' M) {f = chrp.*f;

% chirp convolution9 v/ v9 j5 \4 v* A0 {. b
c = pi/N/sina/4;
7 k: _, s/ }; D! i3 F0 b0 WFaf = fconv(exp(i*c*(-(4*N-4):4*N-4)'.^2),f);
' x) Y0 L0 A1 q) R% d  WFaf = Faf(4*N-3:8*N-7)*sqrt(c/pi);

% chirp post multiplication
1 C( ]% I% b8 ~& p3 VFaf = chrp.*Faf;

% normalizing constant
, Q$ w3 _+ u  N  K8 EFaf = exp(-i*(1-a)*pi/4)*Faf(N:2:end-N+1);

%%%%%%%%%%%%%%%%%%%%%%%%%
  e! \+ z% F: P5 mfunction xint=interp(x)
# A+ `, t9 t. O8 C% sinc interpolation

N = length(x);
5 h, `- e7 T! Jy = zeros(2*N-1,1);
; g! v3 X$ p- [. Y4 j' K* ry(1:2:2*N-1) = x;6 u  G" u$ Z0 c4 V5 {) x# O) Q& R
xint = fconv(y(1:2*N-1), sinc([-(2*N-3)2*N-3)]'/2));$ S5 K) C: O$ {4 G. [0 z: {
xint = xint(2*N-2:end-2*N+3);

%%%%%%%%%%%%%%%%%%%%%%%%%- r# E% F% L+ n5 Y1 [/ y4 I) \7 X
function z = fconv(x,y)/ L7 n3 ?2 c) Z! S
% convolution by fft

N = length([x(;y(])-1;
* I! |$ M9 u* Z1 H6 |" S! LP = 2^nextpow2(N);
& @& F, P3 q& h+ Oz = ifft( fft(x,P) .* fft(y,P));
) v$ D; V! K6 p& z7 B9 jz = z(1:N);


9 {  U' \1 i: m' g5 y* Y4 v- K# t+ _; n
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-10-29 02:03 , Processed in 0.125000 second(s), 23 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表