实验一报告.docx

1、实验一报告离散信号的MATLAB产生和图形显示实验报告目录实验一 21、实验内容 22、编程原理、思路和公式 23、程序脚本，并注释 24、仿真结果、图形 35、结果分析和结论 46、遇到的问题、解决方法及收获 4实验二 51、实验内容 52、编程原理、思路和公式 53、程序脚本，并注释 54、仿真结果、图形 75、遇到的问题、解决方法及收获 7实验三 81、实验内容 82、编程原理、思路和公式 83、程序脚本，并注释 84、仿真结果、图形 115、遇到的问题、解决方法及收获 13实验一1、实验内容编写程序产生序列其中画出序列的图形，说明信号时域振荡程度的变化，理解数字信号的“低频”和“高频”

2、的概念。2、编程原理、思路和公式1、Matlab自带余弦函数，只需调用函数即可，即。2、因为有不同的的值，为了简便，用一个一行八列的矩阵来表示的所有的取值，然后用一个for循环来调用。3、程序脚本，并注释clc;clear %清除原有的东西n=0:1:99; %定义n的范围w0=0,pi/8,pi/4,pi/2,pi,7*pi/4,15*pi/8,2*pi; % 表示w0的矩阵for (i=1:4) %for循环语句来调用不同的w0 x=cos(w0(i)*n);%定义序列 subplot(4,1,i); %一行四列方便观察 stem(n,x,k.),grid on;%画图 xlabel(n)

3、;%横轴为n ylabel(x(n);%纵轴为x（n） title(x(n)=cos(w0*n);%图片的名字为x(n)=cos(w0*n)endfigure %另作一幅图for (i=5:8) %for循环语句来调用不同的w0 x=cos(w0(i)*n); %定义序列 subplot(4,1,(i-4); %一行四列方便观察 stem(n,x,k.),grid on;%画图 xlabel(n);%横轴为n ylabel(x(n);%纵轴为x（n） title(x(n)=cos(w0*n);%图片的名字为x(n)=cos(w0*n)end4、仿真结果、图形图一：NO1仿真图5、结果分析和结论

4、从图中可以看出，从0频率越来越高，频率越来越低，在时频率最高。6、遇到的问题、解决方法及收获问题及解决：1、最开始挨个给赋值，但是那样程序很长很冗杂，然后我就想到用for循环来做2、最初解决循环时，直接用去乘n，但是是1*7的矩阵，n是1*100的，两者没法相乘，后来用for循环，挨个调用，这样就可以了。3、作在同一个图里，纵轴很短，不好看，所以用figure指令，把图画在了两张图上。收获：1、调用一个数组的各个元素，需要用for循环2、同一个程序里画多张图，用figure指令3、把默认的画图里的圈圈换成点用.指令，换颜色也可以，如k.是黑色，r.是红色实验二1、实验内容画出下列序列的图形，进

5、行计算并用图形显示结果。计算2、编程原理、思路和公式x(n)是三角函数，可以直接调用，但是y(n)是矩形序列和单位冲击序列的和，没有现成的序列，所以先用子函数写好y(n)的两个函数，然后调用函数。y的取值范围是-11到16，但是要求x是从-5到25的，两序列相加是需要统一长度，需要长度统一和补零3、程序脚本，并注释function R_N = s2z(ny1,N,m ,ny,L )R_N=zeros(1,L);%给R_N赋初值R_N(-ny1+1-N+1-m:-ny1+1-m)=1;%定义R_10(-n-m）endfunction delta = sy2z( ny1,L,m ,ny )delt

6、a=zeros(1,L); %给delta赋初值delta(-ny1+17)=1;%定义单位冲击函数endclear all; %清除全部nx1=-5; nx2=25; nx=nx1:nx2; %规定nx的范围x=sin(pi/3*nx)+3*cos(nx*pi/4);%x的函数subplot(3,1,1); %三行一列x（n）在第一个stem(nx,x,k.); %画图xlabel(n);ylabel(x(n); title(x(n)=sin（pi/3*n)+3*cos(n*pi/4);%x(n)的横纵轴及图片的名字 ny1=-12; ny2=25; ny=ny1:ny2; %规定ny的范围

7、m=2; %R_N(-n-m)m的取值L=length(ny);%ny的长度N=10; %R_N的长度R_N = s2z(ny1,N,m ,ny,L );%调用矩形序列 delta = sy2z( ny1,L,m ,ny );%调用冲击序列y=R_N+delta;%y的函数subplot(3,1,2); %三行一列y（n）在第2个stem(ny,y,k.); %画图xlabel(n);ylabel(y(n);title(y(n)=R_1_0(-n-2)+delta(n-16);%y(n)的横纵轴及图片的名字 nz1=min(nx1,ny1); %nz的起点为ny和的最小值nz2=max(nx2

8、,ny2); %nz的终点为ny和的最大值nz=nz1:nz2;%统一x和y的长度为nz xa=zeros(1,abs(nz1-nx1) x zeros(1,abs(nz2-nx2);%补零定义新的函数xaya=zeros(1,abs(nz1-ny1) y zeros(1,abs(nz2-ny2);%补零定义新的函数yaz=xa+ya; %两函数相加subplot(3,1,3); %三行一列z（n）在第3个 stem(nz,z,k.); %画图xlabel(n);ylabel(z(n); title(z(n)=x(n)+y(n);%z(n)的横纵轴及图片的名字4、仿真结果、图形5、遇到的问题、

9、解决方法及收获问题及解决：1、R_N(-ny1+1-N+1-m:-ny1+1-m)=1，括号内的是数组内元素是第几个，开始误以为是n的取值，后来查书改正。2、X和y两数组相加必须长度相等，开始没有注意，频繁报错，后来意识到后用老师PPT给的补零相加的方式相加。3、函数定义时忘记了函数名前后不许一致，在错误提示下改正。4、画图时自带的那个带蓝圈圈很大，图看起来很奇怪，查matlab的书k.改成黑色的点，这样看起来方便多了。收获：1、两数组相加必须是长度相等。2、数组定义只能指出是第几个元素，而不能指出n的取值。实验三1、实验内容画出四种理想滤波器的幅频特性图（一个周期）和单位冲激响应波形图（根据

10、公式画图）。其中角频率间隔，时间间隔2、编程原理、思路和公式1、根据四种理想滤波器的幅频特性函数直接画图，用for循环写数组，用if条件语句规定0和1的点，然后计算其单位冲击响应函数，然后用函数直接画图。2、为了好看，画成四张图，每种滤波器及其单位冲击响应画在同一张图上3、画单位冲击响应时，因为零点没有意义，所以需要另外定义赋值零点3、程序脚本，并注释w=-pi:(2*pi/256):pi;%规定w的取值范围及间隔n1=-50;n2=50;n=n1:1:n2;%规定n的取值范围及间隔wc1=0.3*pi;wc2=0.8*pi;w1=0.5*pi;w2=0.7*pi;w3=0.65*pi;w4=

11、0.8*pi;%一些频率值W=length(w);%W为w的长度 for (i=1:W)%for循环表示理想低通滤波器幅频特性 if abs(w(i)=wc1 x1(i)=1;%当|w|=wc2&abs(w(i)=pi x2(i)=1;%当wc2=|w|=w1&abs(w(i)=w2 x3(i)=1;%当w1=|w|=0&abs(w(i)=w4&abs(w(i)=pi x4(i)=1;%当|w|=w3|w4=|w|=w4时x4=1 else x4(i)=0;%当条件不满足时x4=0 end%if语句结束end;%for语句结束subplot(2,1,1);plot(w,x4);axis(-4

12、4 0 1.5);%规定横纵轴的取值范围xlabel(w(rad);ylabel(h_l_p(exp(j*w); title(理想带阻滤波器幅频特性) y4=(-sin(w4*pi*n)+sin(pi*n)+sin(w3*pi*n)./(pi*n);%理想带阻滤波器的单位冲击响应y4(-n1+1)=0.85;%补充n=0的点的值subplot(2,1,2);stem(n,y4,.);xlabel(n);ylabel(h_l_p(n); title(单位冲击响应h_l_p(n)=(-sin(0.8*pi*n)+sin(pi*n)+sin(0.65*pi*n)./(pi*n);) 4、仿真结果、图形5、遇到的问题、解决方法及收获问题、解决1、矩阵相除用“./”，开始忘记了，直接用了“/”使程序出错。2、用if语句时直接用y=1,这只表示一个数，无法表示一个数组，而y是一个数组，所以用for循环来表示一个数组。3、开始忘记补零，使n=0时，y=0。后来补上就好了收获1、学会了用for循环和if条件语句2、学会了画连续函数3、学会了矩阵的相除4、全面的了解四种理想的滤波器及其单位冲击响应