说明.txt

％本程序中时间单位是毫秒 ％频率单位为kHz ％本例说明取样点数与仿真精度的关系 跟在％符号后面的是程序注释

%本程序中时间单位是毫秒
%频率单位为kHz
%本例说明取样点数与仿真精度的关系
global dt df t f N
定义变量dt，df，t，f，N为全局变量：一般由M文件定义的变量是局部变量，这些变量只在该函数内有效，它们与其他函数的变量和工作区中的变量是相互独立的。如果在函数，M文件或工作区中表明了某个变量为全局变量，那么它们就可为其他函数所共享。

figure(1)
    set(1,'Position',[10,50,500,200])
                %设定图1的窗口位置及大小 
   figure(2)
   set(2,'Position',[410,50,300,200])
                %设定图2的窗口位置及大小 
    打开图形窗口，并对图形窗口的大小和位置进行设定。如果没有该语句，打开的图形窗口的大小和位置是MATLAB默认的，并且所有的图形窗口重叠显示。

k=input('取样点数＝2^k, k=[10]');
   屏幕输入语句，程序运行到此，暂停，等待用户输入k的数值。单引号内的字符串会在工作区窗口显示。
   if k==[ ], k=10; end
    条件执行语句，如果用户没有在工作区输入数值，而是直接回车，即输入一个空向量，那么系统自动将10这个数赋值给k。（因为版本的关系，当可k＝[ ]时，程序不给k赋值为10，我们仍需在工作区输入数值）
   N=2^k

N=2^k
设置总取样点数N的值为2的整次幂
dt=0.02;          %ms
设置时域取样间隔
df=1/(N*dt)      %kHz
频域取样间隔，为了保证时域和频域的取样点数相
等，df 要根据dt和N设定。
T=N*dt            %截短时间
信号截短时间，N和dt确定后T即确定。
Bs=N*df/2       %系统带宽
N和df确定后Bs即确定

f=[-Bs+df/2:df:Bs];    %频域横坐标
t=[-T/2+dt/2:dt:T/2];   %时域横坐标
信号截短时间为T，定义域为[-T/2,T/2]，横坐标以  
-T/2+dt/2为起始点,T/2为终止点，dt为步长。
思考：为什么起始点要加上dt/2？
提示：此时的横坐标实际是仿真的抽样时刻。
该语句也可替换成linspace(-T/2,T/2,N)，即将截短
信号等间隔分成(N-1)份。
linspace(a,b,N)=[a:(b-a)/(N-1):b]