qiuzhi.m

采用分段设置采样速率的方式

function [A,P,A1,P1,A0,P0,y]=qiuzhi(S,N,num,den)
%sysd=c2d(sys,1/Fs,'tustin');
%[num,den]=tfdata(sysd,'v');
y=filter(num,den,S);
ham=hamming(N+1);%加窗函数
%S=S.*ham';
%y=y.*ham';
Y = fft(S,N); %做FFT
Y0=fft(y,N);
A1 =abs(Y(1:N/2)); %取模
A1=A1/(N/2);   %换算成实际的幅度加%%
A1(1)=A1(1)/2;
P1=[1:N/2];
for i=1:N/2
 P1(i)=phase(Y(i)); %计算相位
 P1(i)=P1(i)*180/pi; %换算为角度
end;
A0 = abs(Y0(1:N/2)); %取模
A0=A1/(N/2);   %换算成实际的幅度
A0(1)=A1(1)/2;
for i=1:N/2   
 P0(i)=phase(Y0(i)); %计算相位
 P0(i)=P0(i)*180/pi; %换算为角度
end;
A=[1:N/2];
P=[1:N/2];
for i=1:N/2
 Y(i)=Y0(i)/Y(i);
end;
A=abs(Y(1:N/2)); %取模
for i=1:N/2   
 P(i)=phase(Y(i)); %计算相位
 P(i)=P(i)*180/pi; %换算为角度
end;