用MATLAB编程判定系统的线性特性.doc

课 程 设 计 课程设计名称：数字信号处理课程设计 专 业 班 级 ： 学 生 姓 名 ： 学 号 ： 指 导 教 师 ： 课程设计时间： 电子信息工程 专业课程设计任务书 学生姓名 专业班级 学号 题 目 用MATLAB编程判定系统的线性特性 课题性质 其他 课题来源 自拟 指导教师 同组姓名 主要内容 已知一个系统为： ，输入3个不同的序列、、 用MATLAB编程求出、和，画图并判定该系统是否是线性系统（a,b可以为任意实数）。 任务要求 1. 掌握线性系统的特性。 2. 掌握用MATLABMATLAB画图 参考文献 1．程佩青著，《数字信号处理教程》，清华大学出版社，2001 2．Sanjit K. Mitra著，孙洪，余翔宇译，《数字信号处理实验指导书（MATLAB版）》，电子工业出版社，2005年1月 3．郭仕剑等，《MATLAB 7.x数字信号处理》，输入3个不同的序列、、 用MATLAB编程求出、和，画图并判定该系统是否是线性系统（a,b可以为任意实数）。 二　设计原理及设计思路 对系统输入两个序列x1[n]、 x2[n]对应的输出结果为y1[n]、y2[n],在某一权系数下，输入序列x3[n]为x1[n] 与x2[n]的加权和，其对应的输出为y3[n]，在相同的权系数下对y1[n]、y2[n]进行加权求和得出yt[n]，对y3[n]和yt[n]进行差异比较分析，通过结果d[n]来判定系统是否为线性系统，如果y3[n]和yt[n]差异很小，则可判定该系统为线性系统，如果y3[n]和yt[n]差异较大，则判定该系统为非线性系统。 三　程序流程图 a b 四　程序源代码 % Program P2_3 % Generate the input sequences clf; n = 0:40; a = 2;b = -3; x1 = cos(2*pi*0.1*n); x2 = cos(2*pi*0.3*n); x = a*x1 + b*x2; num = [2.2403 2.4908 2.2403]; den = [1 -0.4 0.75]; ic = [0 0]; % 初始条件 y1 = filter(num,den,x1,ic); % 计算输入x1[n]对应的输出y1[n] y2 = filter(num,den,x2,ic); % 计算输入x2[n]对应的输出y2[n] y = filter(num,den,x,ic); % 计算加权输入对应的输出y[n] yt = a*y1 + b*y2; %加权求和 d = y - yt; % 计算不同的输出 d[n] % 图的输出和信号的差别 subplot(3,1,1) stem(n,y1); ylabel(Amplitude); title(Output Due to Input: x_{1}[n]); subplot(3,1,2) stem(n,y2); ylabel(Amplitude); title(Output Due to Input: x_{2}[n]); subplot(3,1,3) stem(n,y3); ylabel(Amplitude); title(Output Due to Weighted Input: a \cdot x_{1}[n] + b \cdot x_{2}[n]); 五仿真结果图 在权系数为a=2、b= -3初始条件ic=[0 0]时，下图分别为：加权输入得到的输出y3[n]；在相同权系数下 输出y1[n]和y2[n]相加和yt[n]；以及y3[n]与yt[n]的差异比较结果d[n] 结果分析：因为数量级是在10- 14，可以看出y3[n]与yt[n]二者差异不大，故系统可以判定为线性系统。 换一组系数a=-3和b=2运行一下，运行结果如下图： 结果分析：在10-14数量级条件下y3[n]与yt[n]二者差别很小，故可以判定为线性系统。