%實驗一
%余弦信號的傅立葉變換
global dt df N t f T %全局變量
close all
k=input('取樣點數=2^k,k=[10]');
if k==[],k=10;end
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);
f0=1;
s=cos(2/3*pi*f0*t);
S=t2f(s);
a=f2t(S);
figure(1)
set(1,'position',[10,50,300,200])%設定窗口位置及大小
figure(2)
set(2,'position',[350,50,300,200]);%設定窗口位置及大小
figure(1)
as=abs(S);  %求模
plot(f,as,'g')
grid;
axis([-2*f0,+2*f0,min(as),max(as)]);
xlabel('f(kHZ)');
ylabel('S(f)(V/kHZ)');
figure(2)
plot(t,a,'b')
grid
axis([-4,+4,-1.5,1.5']);
xlabel(('t(ms)'));
ylabel('s(t)(V)');