%例一 正弦信號

%本程序中時間單位是毫妙
%頻率單位為KHz
%本例說明取樣點數與仿真精度的關系

%global dt df t f N
close all

figure(1)
 set(1,'Position',[10,50,300,200])
                %設定圖1的窗口位置及大小 
figure(2)
 set(2,'Position',[410,50,300,200])
                %設定圖2的窗口位置及大小 

k=input('取樣點數＝2^k, k=[10]');
if k==[], k=10; end
N=2^k

dt=0.02;        %ms
df=1/(N*dt);     %KHz
T=N*dt;          %截短時間
Bs=N*df/2;      %系統帶寬

f=[-Bs+df/2:df:Bs];   %頻域橫坐標
t=[-T/2+dt/2:dt:T/2];   %時域橫坐標

f0=1;   %f0=1KHz

 s=sin(2*pi*f0*t);
 S=t2f(s);    %S是s的傅氏變換
 a=f2t(S);
 P=S.*conj(S)/T;    %功率譜


 figure(1)
 as=abs(S);
 plot(f,as,'b');
 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([-2,+2,-1.5,1.5'])
 xlabel('t (ms)')
 ylabel('s(t) (V)')

 P1=sum(s.*s)*dt/T  %時域計算的功率
 P2=sum(P)*df	    %頻域計算的功率