数字信号处理matlab实验资料
数字信号处理实验报告-FFT算法的MATLAB实现
数字信号处理 实验报告实验二：FFT 算法的MATLAB 实现一、实验目的通过本实验的学习，掌握离散傅立叶变换的理论，特别是FFT 的基本算法以及其在在数字信号处理中的应用。
二、实验内容题一：若x(n)=cos(n*pi/6)是一个N =12的有限序列，利用MATLAB 计算它的DFT 并画出图形。
题二：一被噪声污染的信号，很难看出它所包含的频率分量，如一个由50Hz 和120Hz 正弦信号构成的信号，受均值随机噪声的干扰，数据采样率为1000Hz ，对这污染信号进行傅立叶变换，以检查所包含的频率分量题三：调用原始语音信号mtlb ，对其进行FFT 变换后去掉幅值小于1的FFT 变换值，最后重构语音信号。
（要求有四幅语音信号的频谱图在同一图形窗口以便比较：分别是1、原始语音信号；2、FFT 变换；3去掉幅值小于1的FFT 变换值；4、重构语音信号）三、实验原理1、有限长序列x(n)的DFT 的概念和公式：⎪⎪⎩⎪⎪⎨⎧-≤≤=-≤≤=∑∑-=--=101010)(1)(10)()(N k kn N N n kn N N n W k x N n x N k W n x k x2、基2的FFT 算法四、实验条件（1）微机（2）MATLAB 编程工具五、用matlab 程序实现：实验一：clc;N=12;n=0:N-1;k=0:N-1;xn=cos(n*pi/6);W=exp(-j*2*pi/N);kn=n'*kXk=xn*(W.^kn)stem(n,Xk);xlabel('k');ylabel('Xk');实验二clc;fs=1000;N=1024;n=0:N-1;t=n/fs;x=sin(2*pi*50*t)+sin(2*pi*120*t)+rand(1,N);y=fft(x,N);mag=abs(y);f=n*fs/N;subplot(1,2,1),plot(f,mag);xlabel('频率/Hz');ylabel('振幅');title('N=1024');subplot(1,2,2),plot(f(1:N/2),mag(1:N/2)); xlabel('频率/Hz');ylabel('振幅');title('N=1024');实验三：clc;load mtlbN=1500;subplot(2,2,1)plot([1:N],mtlb(1:N));title('原始语音信号');y=fft(mtlb(1:N));subplot(2,2,2)plot([1:N],y);title('原始语音信号FFT变换');y(y<1)=0;subplot(2,2,3)plot([1:N],y);title('去掉幅值小于1的FFT变换值');subplot(2,2,4)plot([1:N],ifft(y));title('重构语音信号');六、实验心得：(1)通过这次实验，掌握离散傅立叶变换的理论，特别是FFT的基本算法以及其在在数字信号处理中的应用。
Matlab数字信号处理——综合实验报告
实验题目：用Matlab语言编程实现学生成绩管理系统一、实验目的1.掌握建立和执行M文件的方法，学会M函数的编写及其调用，掌握多种函数类型；2.学习Matlab中的关系运算和逻辑运算，掌握它们的表达形式和用法；3.掌握Matlab中的选择结构和循环结构；4.学习使用帮助文档查找更多的Matlab功能；5.掌握Matlab的基本绘图命令，学习运用Matlab绘制一维、二维、三维图形的方法。
二、实验要求1.使用函数（多个函数M文件）进行分析和设计；2.要使用到以下函数类型：主函数、子函数、私有函数。
三、实验内容功能内容：1.学生信息为：学号、姓名、三门课成绩（英语、数学、物理）；2.从键盘输入学生成绩，建立信息库；3.从文件读入学生信息；4.显示信息库中的所有学生的信息；5.将学生成绩信息库保存为文件；6.根据学号查询某一个学生的成绩统计信息：三门课成绩、总分、平均分；7.根据学号将某一个学生的信息从信息库删除；8.新增一个学生的信息；9.按总分成绩从高到低排序输出学生信息；10.图形输出所有学生成绩统计信息（饼图pie、柱状图bar）。
四、实验设计及分析1 主函数的设计及分析主函数主要用于显示学生管理系统的主界面，以及通过对其他几个私有函数的调用来实现相应的功能。
主函数名为zonghe()，首先利用disp()屏幕输出函数来提示用户如何使用该管理系统。
然后通过input()输入函数输入字符，通过判断输入的字符与哪个私有函数对应的字符匹配来调用相应的私有函数实现相应的功能。
其中将条件语句放入while循环中，使系统能循环运行直到用户输入‘H‘才跳出循环。
主函数zonghe()的源代码：function zonghe()disp('Welcome to the students'' information library!');disp('Input A to input students'' information');disp('Input B to display all students'' information');disp('Input C to search a stuedent''s information');disp('Input D to delete a stuedent''s information');disp('Input E to add a new stuedent''s information');disp('Input F to display stuedents'' information from high to low of total score');disp('Input G to graphical display all stuedents'' information'); disp('Input H to exit');x=input('Please input the FUNCTION NUMBER:','s');while x~='H'if x=='A'disp('Input students'' information.');ztest();elseif x=='B'disp('Display all students'' information.');displaytest();elseif x=='C'disp('Search a stuedent''s information.');search();elseif x=='D'disp('Delete a stuedent''s information.');shanchu();elseif x=='E'disp('Add a new stuedent''s information.');add();elseif x=='F'disp('Display stuedents'' information from high to low of total score.');paixu();elseif x=='G'disp('Graphical display all stuedents'' information.');tuxing();elseif x=='H'disp('Thank you for using the students'' information library!');break;elsedisp('Invalid function number. Program is terminated!');endx=input('\nPlease input the FUNCTION NUMBER:','s');endend2 七个私有函数的设计与分析2.1 function ztest()的设计与分析ztest()函数实现从键盘输入学生成绩建立信息库的功能。
数字信号处理 Matlab实验一 Matlab 基本功能和基础知识操作
温州大学物理与电子信息工程学院Matlab 仿真及其应用 实验报告实验一Matlab 基本功能和基础知识操作 [实验目的和要求]1、 熟练掌握Matlab 的启动与退出2、 熟悉Matlab 的命令窗口、常用命令、帮助系统3、 熟悉Matlab 的数据类型、基本矩阵操作、运算符和字符串处理[实验内容]1、 用逻辑表达式球下列分段函数的值 22201112,=0:0.5:2.52123t t y t t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩其中2、 求[100,999]之间能被32整除的数的个数3、 建立一个字符串向量，删除其中的小写字母。
4、 输入矩阵1234514789A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，求出此矩阵的行列式，逆和特征根，并找出A 中大于5和小于9的元素，用行列式表示。
5、 不采用循环的形式求出和式63230034ii i i S ===+∑∑6、 给定矩阵E=rand （4，4），计算C+E ，C*E ，C\E实验结果及分析：经过Matlab 软件的程序编辑和测试分析，得出以下实验结果： 详见程序代码、注释及屏幕截图：【题1】程序代码：t=0:0.5:2.5y=t.^2.*((t>=0)&(t<1))+(t.^2-1).*((t>=1)&(t<2))+(t.^2-2*t+1).*((t>=2)&(t<3)) 效果截图：【题2】程序代码：p=rem([100:999],32)==0;sum(p)效果截图：【题3】程序代码：ch='dfghjGUIJKVC',k=find(ch>'a'&ch<='z'),ch(k)=[]效果截图：【题4】程序代码：A=[1 2 3;4 5 14;7 8 9];[i,j]=find(A>5&A<9) %定位for n=1:length(i)m(n)=A(i(n),j(n))endDA=det(A) %行列式IA=inv(A) %逆矩阵EA=eig(A) %特征根效果截图：【题5】程序代码：E=rand(4,4); %产生随机数0~1 C=rand(4,4);B1=C+E;B2=C*E;B3=C/E;B1B2B3效果截图：【题6】程序代码：E=rand(4,4); %产生随机数0~1 C=rand(4,4);B1=C+E;B2=C*E;B3=C/E;B1B2B3效果截图：4、心得：通过本次Matlab课程实验，我已熟练Matlab的命令窗口、常用命令、帮助系统，并掌握Matlab的数据类型、基本矩阵操作、运算符和字符串处理。
数字信号处理实验报告-离散时间序列卷积和MATLAB实现
数字信号处理 实验报告实验一：离散时间序列卷积和MATLAB 实现一、实验目的通过本实验学会用MATLAB 对信号与系统进行分析，实现离散序列卷积和的计算。
二、实验仪器：（1）微机（2）MATLAB 编程工具 三：实验内容题一：已知序列f1(k)=⎩⎨⎧≤≤其它0201k f2(k)=⎪⎪⎩⎪⎪⎨⎧===其它332211k k k调用conv()函数求上述两序列的卷积和题二：编写计算两离散序列卷积和f(k)=f1(k)*f2(k)的实用函数dconv().要求该程序在计算出卷积和f(k)的同时，还绘出序列f1(k),f2(k)和f(k)的时域波形图。
function[f,k]=dconv(f1,f2,k1,k2)%f1(k),f2(k)及f(k)的对应序号向量分别为k1，k2和k 。
题三：试用MATLAB 计算如下所示序列f1(k)与f2(k)的卷积和f(k)，绘出它们的时域波形，并说明序列f1(k)与f2(k)的时域宽度与序列f(k)的时域宽度的关系。
f1(k)=⎪⎪⎩⎪⎪⎨⎧==-=其它110211k k kf2(k)=⎩⎨⎧≤≤-其它221k提示：可用上述dconv()的函数来解决。
题四：已知某LTI 离散系统，其单位响应h(k)=e(k)-e(k-4)，求该系统在激励为f(k)=e(k)-e(k-3)时的零状态响应，并绘出其时域波形图。
四、实验原理、方法和手段1、离散时间序列f1(k)和f2(k)的卷积和定义：f(k)=f1(k)*f2(k)=∑∞-∞=-•iikfif)(2)(12、在离散信号与系统分析中有两个与卷积和相关的重要结论：a、f(k)= ∑∞-∞=-•iikif)()(δ=f(k)* δ(k)即离散序列可分解为一系列幅度由f(k)决定的单位序列δ(k)及其平移序列之积。
b、对线性时不变系统，设其输入序列为f(k)，单位响应为h(k)，其零状态响应为y(k)，则有：y(k)= ∑∞-∞=-•iikhif)()(五、matlab程序实现卷积和计算：1、题一(conv.m)clc;k1=3;k2=3;k=k1+k2-1;f1=[1,1,1];f2=[0,1,2,3];f=conv(f1,f2);nf=0:k;stem(nf,f,'*r'); xlabel('n'); ylabel('f(n)'); grid on ;2、题二 (dconv.m)function[f,k]=dconv(f1,f2,k1,k2);f=conv(f1,f2);k(1)=k1(1)+k2(1);k(length(f))=k1(length(f1))+k2(length(f2));k=k(1):k(length(f));subplot(1,3,1); stem(k1,f1,'*r'); xlabel('k1'); ylabel('f1(k1)'); grid on ;subplot(1,3,2); stem(k2,f2,'*r'); xlabel('k2'); ylabel('f2(k2)'); grid on ;subplot(1,3,3); stem(k,f,'*r'); xlabel('k'); ylabel('f(k)'); grid on ;3、题三(san.m)clc;f1=[1,2,1];f2=[1,1,1,1,1];k1=-1:1;k2=-2:2;dconv(f1,f2,k1,k2);4、题四 (si.m) clc;f1=[1,1,1,1];f2=[1,1,1];k1=0:3;k2=0:2;dconv(f1,f2,k1,k2);六、实验心得：(1)通过这次实验，让我知道序列的卷积和可以用matlab编程实现，并且实现的方法比较简单，而且matlab本身有丰富的函数库，画图、求卷积和时只要调用相应的函数即可。
MATLAB实验报告模板
MATLAB实验报告模板实验题目：使用MATLAB进行数字信号处理实验目的：1. 学习MATLAB的基本操作和数字信号处理中常用的函数。
2. 掌握数字信号的离散化、采样、量化等处理方法。
3. 实现数字滤波器的设计和应用。
实验内容：1. 生成并绘制一个正弦波信号，包括频率、幅度和相位等参数。
2. 对信号进行采样，并绘制采样后的离散信号图像。
3. 对采样后的信号进行量化，并绘制量化后的信号图像。
4. 设计一个数字滤波器，并将滤波前后的信号图像进行对比。
实验步骤：1. 生成一个正弦波信号：```matlab% 信号频率为3Hz，幅度为2，相位为0t = 0:0.01:1;f = 3;A = 2;theta = 0;x = A * sin(2*pi*f*t + theta);```2. 绘制信号图像：```matlabplot(t, x);xlabel('时间');ylabel('幅度');title('正弦波信号');```3. 进行信号的采样：```matlab% 采样频率为20Hzfs = 20;Ts = 1/fs;n = 0:Ts:1;xs = A * sin(2*pi*f*n + theta);```4. 绘制采样信号图像：```matlabstem(n, xs);xlabel('时间');ylabel('幅度');title('采样信号');```5. 进行信号的量化：```matlab% 将信号量化为8位bits = 8;delta = (2 * A) / (2^bits);xq = floor(xs / delta + 0.5) * delta;```6. 绘制量化信号图像：```matlabstem(n, xq);xlabel('时间');ylabel('幅度');title('量化信号');```7. 设计数字滤波器：```matlab% 采用FIR滤波器，截止频率为4Hz fcut = 4;n = 100; % 滤波器阶数b = fir1(n, 2*fcut/fs);y = filter(b, 1, xq);```8. 绘制滤波前后的信号图像：```matlabfigure;subplot(2,1,1);stem(n, xq);xlabel('时间');ylabel('幅度');title('量化信号');subplot(2,1,2);stem(n, y);xlabel('时间');ylabel('幅度');title('滤波信号');```实验结果和分析：根据以上实验步骤，可以得到正弦波信号、采样信号、量化信号和滤波信号等图像。
matlab dsp实验报告
MATLAB DSP实验报告介绍本实验报告将详细介绍在MATLAB环境下进行数字信号处理（DSP）的实验步骤和相关方法。
我们将通过逐步思考的方式，帮助读者理解和学习DSP的基本概念和技术。
实验环境和工具在进行DSP实验之前，我们需要准备以下环境和工具：1.MATLAB软件：确保已安装并配置好MATLAB软件，可以在MATLAB Command窗口中输入命令。
2.信号处理工具包：在MATLAB中，我们可以使用信号处理工具箱（Signal Processing Toolbox）来进行DSP实验和分析。
确保该工具箱已被安装并加载。
实验步骤下面是进行DSP实验的一般步骤：步骤一：加载信号首先，我们需要加载待处理的信号。
这可以通过在MATLAB中使用load命令加载一个音频文件或生成一个模拟信号实现。
例如，我们可以加载一个名为signal.wav的音频文件：load signal.wav步骤二：信号预处理在进行DSP之前，通常需要对信号进行预处理。
这可能包括去噪、滤波、均衡等操作。
例如，我们可以使用滤波器对信号进行降噪：filtered_signal = filter(filter_coefficients, signal);步骤三：信号分析一旦信号经过预处理，我们可以开始进行信号分析。
这可能涉及频域分析、时域分析、谱分析等。
例如，我们可以通过计算信号的快速傅里叶变换（FFT）获得其频谱：spectrum = fft(filtered_signal);步骤四：特征提取在信号分析之后，我们可以根据需要提取信号的特征。
这些特征可能包括幅度、频率、相位等。
例如，我们可以计算信号的能量：energy = sum(abs(filtered_signal).^2);步骤五：信号重构在完成信号分析和特征提取后，我们可以根据需要对信号进行重构。
这可能包括滤波、修复损坏的信号等。
例如，我们可以使用滤波器对信号进行重构：reconstructed_signal = filter(filter_coefficients, filtered_signal);步骤六：结果评估最后，我们需要评估重构后的信号和原始信号之间的差异。
数字信号处理第三版用MATLAB上机实验
实验二：时域采样与频域采样一、时域采样1.用MATLAB编程如下:%1时域采样序列分析fs=1000A=444.128; a=222.144; w=222.144; ts=64*10^(-3); fs=1000;T=1/fs;n=0:ts/T-1; xn=A*exp((-a)*n/fs).*sin(w*n/fs); Xk=fft(xn);subplot(3,2,1);stem(n,xn);xlabel('n,fs=1000Hz');ylabel('xn');title('xn');subplot(3,2,2);plot(n,abs(Xk));xlabel('k,fs=1000Hz'); title('|X(k)|');%1时域采样序列分析fs=200A=444.128; a=222.144; w=222.144; ts=64*10^(-3); fs=200;T=1/fs;n=0:ts/T-1; xn=A*exp((-a)*n/fs).*sin(w*n/fs);Xk=fft(xn);subplot(3,2,3);stem(n,xn);xlabel('n,fs=200Hz'); ylabel('xn');title('xn');subplot(3,2,4);plot(n,abs(Xk));xlabel('k,fs=200Hz'); title('|X(k)|');%1时域采样序列分析fs=500A=444.128; a=222.144; w=222.144; ts=64*10^(-3); fs=500;T=1/fs;n=0:ts/T-1; xn=A*exp((-a)*n/fs).*sin(w*n/fs); Xk=fft(xn);subplot(3,2,5);stem(n,xn);xlabel('n,fs=500Hz');ylabel('xn');title('xn');subplot(3,2,6);plot(n,abs(Xk));xlabel('k,fs=500Hz'); title('|X(k)|');2.经调试结果如下图:20406080-200200n,fs=1000Hzxnxn2040608005001000k,fs=1000Hz|X (k)|51015-2000200n,fs=200Hzx nxn510150100200k,fs=200Hz |X(k)|10203040-2000200n,fs=500Hzx nxn102030400500k,fs=500Hz|X (k)|实验结果说明:对时域信号采样频率必须大于等于模拟信号频率的两倍以上,才 能使采样信号的频谱不产生混叠.fs=200Hz 时，采样信号的频谱产生了混叠，fs=500Hz 和fs=1000Hz 时，大于模拟信号频率的两倍以上，采样信号的频谱不产生混叠。
数字信号处理实验指导书--Matlab版
下：
由于 e jω 是频率的周期函数，所以系统的频率特性也是频率的周期函数，且周期为 2π ，因 此研究系统频率特性只要在[ −π ,π ]范围内就可以了。
三、实验内容及步骤
1、在实验上机前，认真复习离散卷积和差分方程的有关内容，仔细阅读本实验原理和 步骤。准备实验所需数据。
2、离散卷积计算。 1）在主界面下进入“实验五” 的“离散卷积”子系统，选定几组不同的 x(n)、h(n) 代入离散卷积程序计算 y(n)，记录实验结果。 3、差分方程迭代解法 1）¨在主界面下进入“实验五” 的“差分方程”子系统，本实验要求首先确定系统方 程系数个数和需要输出 y(n)的样点个数，然后依次输入 a[0]~a[N]、b[0]~b[N]这些系数的 数值和 y(n)的初始值 y(0)~y(N-1)，h(n)的初始值 h(0)~h(N-1)，运行后即得冲激响应 h(n) 和阶跃信号激励下的响应 y(n)的结果和图形显示。 2）自己选择一个离散时间系统，写出其差分方程，设好初始条件进行实验并记录实验 结果。 4、在主界面下进入“实验六”即“离散系统 Z 域分析”，本实验中给出了计算系统
明了系统结构、参数、特性三者之间的关系，即同一结构，参数不同其特性也不同。 例如，下图所示离散系统:
数学模型由下列差分方程描述： y(n)=ay(n-1)+x(n)
系统函数 系统频率特性
H (z) = z , | z |> a z−a
幅度特性为
H (e jω )
=
e jω e jω − a
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
求下列信号，并画图
y1(n)=x(n)*h1(n), y2(n)=x(n)*h2(n), y3(n)=x(n)*h3(n)
要求
验证线性卷积的线性性质 分析输出和输入之间相对位移与对应系统的h(n)之间的 关系 绘制x(n),h(n)和y(n)。给出分析结果
其FT为
j j 2 j 3 0.8 0.44 e 0.36 e 0.02 e H (e j ) 1 0.7e j 0.45e j 2 0.6e j 3
用MATLAB计算的程序如下： k=256; num=[0.8 -0.44 0.36 0.02]; den=[1 0.7 -0.45 -0.6]; w=0:pi/k:pi; h=freqz(num,den,w); subplot(2,2,1); plot(w/pi,real(h));grid title(‘实部’) xlabel('\omega/\pi');ylabel('幅度')
图1.1给出了卷积结果的图形，求得的结果存放在数 组c中为：{-2 -4 1 3 1 5 1 -3}。
例 2 用 MATLAB 计算差分方程的输出：
已知输入序列 ,求输出 解 MATLAB程序如下： N=41; a=[0.8 -0.44 0.36 0.22]; b=[1 0.7 -0.45 -0.6]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(a,b,x); stem(k,y) xlabel('n');ylabel('幅度')
相关结果的长度与信号长度的关系 值点判断两信号之间的相对时差
例6：利用相关运算 检测两个信号的相对时移
计算例5中y1(n)和y2(n)的互相关，根据互相关峰
计算h1(n)和h2(n)的互相关，根据峰值点位置判
断两信号之间的相对时差 对比上述两个互相关的关系，分析二者之间的关 系
要求
。
1.5
1
0.5
幅度
0 -0.5 -1 0
5
10
15
20 n
25
30
35
40
图 1.2 给出了该差分方程的前41个样点的输出，即该系统的 单位脉冲响应。表达式进原理图输入。
例 3 用 MATLAB 计算差分方程所对应 的系统函数的频率响应
解 例2差分方程所对应的系统函数为：
0.8 0.44 z 1 0.36 z 2 0.02 z 3 H ( z) 1 0.7 z 1 0.45 z 2 0.6 z 3
例1 用MATLAB 下列序列的离散卷积： {-2 0 1 –1 3}，{1 2 0 -1} 解 MATLAB程序如下： a=[-2 0 1 -1 3]; b=[1 2 0M=length(c)-1; n=0:1:M; stem(n,c); xlabel('n'); ylabel('幅度');
位冲激响应和阶跃响应，并绘出其图形。要求分
别用 filter、conv、impz三种函数完成，给出理
论计算结果和程序计算结果并讨论。
y[n] 0.75y[n 1] 0.125y[n 2] x[n] x[n 1] y[n] 0.25{x[n 1] x[n 2] x[n 3] x[n 4]}
数字信号处理实验
实验安排
实验一 序列的基本运算 实验二 快速傅立叶变换（FFT）及其应用 实验三 IIR数字滤波器的设计
实验四 FIR数字滤波器的设计
实验一序列的基本运算
一、实验目的
1．在理论学习的基础上，通过本实验，加深对序列基本运算 的理解；
2．熟悉利用matlab软件做简单的仿真实验；
二、实验原理及内容
subplot(2,2,2); plot(w/pi,imag(h));grid title('虚部') xlabel('\omega/\pi');ylabel('Amplitude') subplot(2,2,3); plot(w/pi,abs(h));grid title('幅度谱') xlabel('\omega/\pi');ylabel('幅值') subplot(2,2,4); plot(w/pi,angle(h));grid title('相位谱') xlabel('\omega/\pi');ylabel('弧度')
实部 6 4 10 6 4 5 2 0 -5 0 -10 -2
虚部
2 0 -2
0
0.5
0.5 1 / 幅度谱
1.5
1 2
Amplitude Amplitude
幅度
幅值
4 2 0
弧度
幅 度
0
0.5
0.5 1 / 相位谱
1.5
1 2
8 6
3 4 2 1 0 -2 0 -1 -4
0
0.5
0.5 1 /
理解matlab显示的相关函数中点与rxy(0)或者ryx(0)之间 的关系
分析相关函数rxy与ryx之间的关系 理解掌握max(rxy )的相对于matlab显示的相关函数中点 的位置关系与参与运算的两信号之间的时差关系 分析输出y1与y2互相关的最大值的位置place_y（相对 于中点），分析输出h1与h2互相关的最大值的位置 place_h（相对于中点） 讨论place_y与place_h的关系，可以得出怎样的结论？ 绘制各种相关函数曲线，给出分析结果
y=filter(p,d,x)实现差分方程的仿真 y=conv(x,h)计算卷积 y=impz(p,d,N)求系统的冲激响应
例5：计算信号的卷积
i=0:1:49 x(n)=sin(2*pi*i/50) h1(n)=[zeros(1,10),1,zeros(1,20)];
h2(n) =[zeros(1,20),1,zeros(1,10)];
1.5
1 2
0
0.5
分析范围： 分析范围： [0, [0, 2π] π]
0.5 1 /
1.5
1 2
end
要求
验证线性卷积的线性性质 分析输出和输入之间相对位移与对应系统的h(n)之间的 关系 绘制x(n),h(n)和y(n)。给出分析结果
例4（选做）
编制程序求解下列两个系统的单