南邮MATLAB软件设计.docVIP

通信与信息工程学院 2015 / 2016 学年 第 2 学期 软件设计 实验报告 模 块 名 称 MATLAB软件设计 专 业 通信工程 学 生 班 级 B130102 学 生 学 号 学 生 姓 名 郭智超 指 导 教 师 孟田甜 设计题目 练习题部分 设计部分： （1）M文件实现部分 （2）Simulink实现部分 任务要求 熟悉MATLAB软件，熟练掌握计算、设计方法，对设计的系统使用Simulink工具中进行仿真。掌握MATLAB部分工具包的使用方法，侧重于simulink环境。通过完成所有的练习题，熟悉M文件的编写，Simulink环境中仿真系统的搭建、运行、调试和验证。在此基础上，完成通信系统设计题，并调试验证通过。 实验设备及软件 PC机一台 MATLAB7.0 同组人员学号及姓名 无 参考文献 《通信原理》第7版，樊昌信、曹丽娜编著，国防工业出版社 《数字信号处理实验指导书》MATLAB版，Sanjit K. Mitra著， 电子工业出版社 《基于MATLAB/Simulink系统仿真权威指南》，王江、付文利等编著，机械工业出版社 报告内容 一、要求练习的实验部分 1.题目：在时间区间 [0,10]中，绘出曲线。 程序代码： t=0:0.02:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,b); xlabel(time); ylabel(amplitude); title(figure 1); grid 实验结果： 题目：写出生成如图E2-1所示波形的MATLAB脚本M文件。图中虚线为正弦波，要求它的负半波被置零，且在处被削顶。 程序代码： t=linspace(0,10,500); %产生线性间隔的向量（在0-10之间取500点） y=sin(t); a=sin(pi/3); z=(y=0).*y; z=(y=a).*a+(ya).*z; plot(t,y,:r); hold on; %保持当前图形 plot(t,z,b-) xlabel(t); ylabel(z=f(t)); title(逐段解析函数); legend(y=sin(t),z=f(t),4); %确定向量长度 axis([0 10 -1.2 1.2]); %限定坐标轴范围 hold off; grid 实验结果： 3.题目： 令，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： sum(diag(S)tol) - rank(A) 是多少 ？ S(1,1) - norm(A) = 0 是多少 ? sqrt(sum(diag(S*S))) - norm(A,fro) 的结果是什么 ? S(1,1)/S(3,3) - cond(A) 的结果是什么 ? S(1,1)*S(2,2)*S(3,3) - det(A) tol 的运行结果是什么？ V(:,1)*null(A) 得到什么结果 ? abs(A*null(A)) tol 得到什么结果 ? U(:,1:2) = = orth(A) 的运行结果是什么 ? 程序： A=[1 2 3;4 5 6;7 8 9] [U,S,V]=svd(A);tol=S(1,1)*3*eps sum(diag(S)tol) - rank(A) S(1,1) - norm(A) sqrt(sum(diag(S*S))) - norm(A,fro) S(1,1)/S(3,3) - cond(A) S(1,1)*S(2,2)*S(3,3) - det(A) tol V(:,1)*null(A) abs(A*null(A)) tol U(:,1:2) == orth(A) 运行结果： sum(diag(S)tol) - rank(A) 的结果是0 S(1,1) - norm(A) = 0的结果是0 sqrt(sum(diag(S*S))) - norm(A,fro) 的结果是3.5527e-015 S(1,1)/S(3,3) - cond(A) 的结果是 -8 S(1,1)*S(2,2)*S(3,3) - det(A) tol的结果是1 V(:,1)*null(A) 的结果是 0 abs(A*null(A)) tol 的结果是1；1；1 U(:,1