MATLA例题.ppt

MATLAB例题 本讲义中黑体字为源程序语句，蓝体字为老师讲义 例一 正弦信号的波形及频谱 %本程序中时间单位是毫秒 %频率单位为kHz %本例说明取样点数与仿真精度的关系 跟在％符号后面的是程序注释 global dt df t f N 定义变量dt，df，t，f，N为全局变量：一般由M文件定义的变量是局部变量，这些变量只在该函数内有效，它们与其他函数的变量和工作区中的变量是相互独立的。如果在函数，M文件或工作区中表明了某个变量为全局变量，那么它们就可为其他函数所共享。 例一 正弦信号的波形及频谱 close all 关闭所有的图形窗口，这条语句是为了让前面执行的程序所打开的窗口全部关闭。 例一 正弦信号的波形及频谱 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] 例一 正弦信号的波形及频谱 f0=1; %f0=1kHz s=sin(2*pi*f0*t); S=t2f(s); %S是s的傅氏变换 a=f2t(S); %傅氏反变换 t2f，f2t为老师给定的自定义函数，运行程序 时必须和m文件放在同一个目录下。 例一 正弦信号的波形及频谱 figure(1) as=abs(S); 对傅氏变换求模，得到正弦信号的频谱。 plot(f,as,b,LineWidth,2.5); 绘图指令：在图形窗口1中绘图，横座标为f， 纵座标为as，LineWidth代表使用宽线型， 2.5为线宽指标，’b’表示线型颜色，不特意指 定颜色及线型则为系统默认颜色和线型。该指令 画出正弦信号的频谱 例一 正弦信号的波形及频谱 grid 在图形窗口中画出格线 axis([-2*f0,+2*f0,min(as)-10,max(as)+10]) 标定坐标范围：方括号中前两个值是图形横座标范 围，后两个值是纵座标范围。 xlabel(f (KHz)) ylabel(S(f) (V/KHz)) 在x轴和y轴上加标签，内容为单引号中的字符串， 目的是对x，y轴加以说明。 例一 正弦信号的波形及频谱 figure(2) plot(t,a,b,LineWidth,2.5) grid axis([-2,+2,-1.5,1.5]) xlabel(t (ms)) ylabel(s(t) (V)) 这组指令画出恢复后的正弦信号，表示从截 短后的信号完全能恢复出原信号 例