%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%階比仿真，2011年3月21日
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear?all
close?all
clc
Fs=256;Ts=1/Fs;N=2048;t=0*Ts:Ts:（N-1）*Ts;
%采樣時間為8秒
x=sin（2*pi*t.^2）;%0.8*randn（size（t））;
y1=（abs（fft（x）））*2/N;f1=（0:N/2-1）*Fs/N;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fc=128;%理想低通濾波器截止頻率
wc=2*fc*pi;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%轉(zhuǎn)速信號，1024轉(zhuǎn)
n=2*60*t;%轉(zhuǎn)速
i=1;
a=zeros（32）;
for?nn=1:32
??%??if?nn==1
??%??a（1）=0;
????%else
????a（nn）=sqrt（nn*2）;
????%end
???%?i=i+1;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
b=[0;2*pi;4*pi];
kk=1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%計算等角度時刻
y=zeros（112）;
for?i=1:2:length（a）-3
????A=[1a（i）a（i）^2;1a（i+1）a（