周期非周期信号傅里叶频谱相谱分析

《周期非周期信号傅里叶频谱相谱分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、dt=0.01;t=-10:dt:10;f=u(t+1/2)-u(t-1/2);dw=0.1w=-4*pi:dw:4*pi;F_w=f*exp(-j*t'*w)*dt;%傅里叶变换Mod_F=abs(F_w);%计算幅度谱Phase_F=angle(F_w);%计算相位谱figure;subplot(2,2,1:2)plot(t,f)axis([-2,2,-2,2])gridontitle('f(t)')subplot(2,2,3)plot(w,Mod_F)gridontitle('

2、F(w)

3、')axistightsubplot(2,2,4

4、)plot(w,Phase_F)gridontitle('phi(w)')axistightdt=0.01;t=-10:dt:10;f=stepfun(t,0)-stepfun(t-1,0);dw=0.1w=-4*pi:dw:4*pi;F_w=f*exp(-j*t'*w)*dt;%傅里叶变换Mod_F=abs(F_w);%计算幅度谱Phase_F=angle(F_w);%计算相位谱figure;subplot(2,2,1:2)plot(t,f)axis([-2,2,-2,2])gridontitle('f(t)')subplot(2,2,3

5、)plot(w,Mod_F)gridontitle('

6、F(w)

7、')axistightsubplot(2,2,4)plot(w,Phase_F)gridontitle('phi(w)')axistightdt=0.01;T1=2;w1=2*pi/T1;t=-T1/2:dt:T1/2;tau=1;A=1;f=A*(u(t+tau/2)-u(t-tau/2));N=10;n=-N:N;Fn=f*exp(-j*t'*w1*n)*dt/T1;Mod_F=abs(Fn);%计算幅度谱Phase_F=angle(Fn);%计算相位谱subplot(2

8、,2,1:2)plot(t,f)title('f(t)')subplot(2,2,3)stem(n,Mod_F)title('

9、F(n)

10、')axistightsubplot(2,2,4)stem(n,Phase_F)title('phi(n)')axistightdt=0.01;T1=2;w1=2*pi/T1;t=-T1/2:dt:T1/2*2;tau=1;A=1;f=A*(u(t+tau/2-0.5)-u(t-tau/2-0.5));N=10;n=-N:N;Fn=f*exp(-j*t'*w1*n)*dt/T1;Mod_F=abs(Fn)

11、;%计算幅度谱Phase_F=angle(Fn);%计算相位谱subplot(2,2,1:2)plot(t,f)title('f(t)')subplot(2,2,3)stem(n,Mod_F)title('

12、F(n)

13、')axistightsubplot(2,2,4)stem(n,Phase_F)title('phi(n)')axistight