%% 通过调用MyFT函数对矩形脉冲信号进行傅立叶变换;
%**************************************************************************
%A:幅值,width:脉宽,dt:时间间隔,T:时间长度,df:频率间隔;
%N:数据点数,fs:采样频率
%2014年3月15日
%作者:zhenbao0501
%**************************************************************************
clear all;clc;
A=10;width=0.5;dt=0.01;T=6;fs=50;N=1024;df=1/(N*dt);
figure;
[t,xt]=MySquare(A,width,dt,T);
[f,Xf]=MyFT(t,xt,fs,df);
subplot(2,2,1);
plot(t,xt,'m') %绘制矩形波脉冲波信号曲线
line1='矩形脉冲时域信号';
line2=['幅值为=',num2str(A),',脉宽为=',num2str(width),'s',...
',时间长度为=',num2str(T),'s'];
title({line1,line2}); %图形标题
xlabel('t(s)'); %X轴注释(t(s))
ylabel('x(t)'); %Y轴注释(x(t))
axis tight; %使坐标轴与数据长度一致
ylim([-1,A+1]); %设定y轴上下限加宽
subplot(2,2,2);
plot(f,Xf,'r'); %绘制矩形脉冲信号幅值频谱曲线
title('矩形脉冲FT双边谱');
xlabel('f(Hz)'); %X轴单位注释(f(Hz))
ylabel('|X(f)|'); %Y轴注释
%% 通过调用MyFT函数对正弦信号进行傅立叶变换;
%**************************************************************************
%A:幅值,width:脉宽,dt:时间间隔,T:时间长度,f:频率,df:频率间隔;
%2014年3月15日
%小组成员:郭金龙(组长),刘倩楠,孙守保
%**************************************************************************
clear all;
A=5;dt=0.001;f=20;fs=50;fai=0;N=1000;df=1/(N*dt);
[t,xt]=MySin(A,f,dt,fai,N);
[f,Xf]=MyFT(t,xt,fs,df);
subplot(2,2,3);
plot(t,xt,'b'); %绘制正弦信号曲线
line1='正弦信号时域曲线';
line2=['幅值A=','{',num2str(A),'}',',频率f=20','Hz'];
title({line1,line2}); %图形标题
xlabel('t(s)'); %X轴注释(t(s))
ylabel('x(t)'); %Y轴注释(x(t))
axis tight; %使坐标轴与数据长度一致
ylim([-A-1,A+1]); %设定y轴上下限加宽
%绘制FT单边谱
fn=length(f);
fa=ceil(fn/2);
subplot(2,2,4);
plot(f(fa:fn),2*Xf(fa:fn),'g'); %绘制信号幅值频谱曲线
title('正弦信号FT单边谱'); %图形标题说明
xlabel('f(Hz)'); %X轴注释(f(Hz))
ylabel('|X(f)|'); %Y轴注释(|X(f)|)
%%
评论0