exampl_1.m

这是正弦信号的图象

%例一 正弦信号

%本程序中时间单位是毫妙
%频率单位为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	    %频域计算的功率