yuxian.m

实现余弦函数及其傅立叶变换后的波形

global dt df N t f T %全局变量
close all
%k=input('取样点数=2^k,k=[10]');k=13;
%if k==[],k=10;end

figure(1);clf;
%set(1,'position',[10,50,300,300])%设定窗口1位置及大小

%figure(2)
%set(2,'position',[350,50,350,300]);%设定窗口2位置及大小

k=20;
N=2^k      %采样点数
dt=0.01;    %时域采样间隔
df=1/(N*dt)%频域采样间隔
T=N*dt;     %截短时间
Bs=N*df/2  %系统带宽
t=linspace(-T/2,T/2,N);
f=linspace(-Bs,Bs,N);2

f0=1;
s=cos(2/3*pi*f0*t);
S=t2f(s);
a=f2t(S);

%figure(1)
as=abs(S);  %求模
subplot(1,2,1);plot(f,as,'g')
grid;
axis([-2*f0,+2*f0,min(as),max(as)]);
xlabel('f(kHZ)');
ylabel('S(f)(V/kHZ)');
title('余弦信号傅式变换曲线');

%figure(2)
subplot(1,2,2);plot(t,a,'b')
grid
axis([-4,+4,-1.5,1.5']);
xlabel(('t(ms)'));
ylabel('s(t)(V)');
title('余弦信号时域波形')