南邮matlab软件设计实验报告.doc

PAGE2 / NUMPAGES32 通信与信息工程学院 2014 /2015 学年 第2学期 软件设计 实验报告 模 块 名 称 MATLAB软件设计 专 业 通信工程 学 生 班 级 通工2班 学 生 学 号 学 生 姓 名 指 导 教 师 孟田甜 设计题目 ..练习题 ..模电仿真 ..数电仿真 ..数字信号处理仿真 ..通信系统仿真 任务要求 熟悉MATLAB软件，熟练掌握计算、设计方法，对设计的系统使用Simulink工具中进行仿真。通过完成所有的练习题，熟悉M文件的编写，Simulink环境中仿真系统的搭建、运行、调试和验证。在此基础上，完成通信系统设计题，并调试验证通过。 实验设备及软件 PC机一台 MATLAB7.0.1 同组人员学号及姓名 无 参考文献 1. 《matlab7.0基础教程》，孙祥编著，清华大学出版社 2. 《MATLAB 通信工程仿真》，张德丰编著，机械工业出版社 报告内容 一、要求练习的实验部分 题目：在时间区间 [0,10]中，绘出曲线。 实现方法：plot画图 程序：clc clear t=0:0.01:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,r-); xlabel(t); ylabel(y); title(绘制曲线); 结果及结果分析： 对源程序的改进：1.加了坐标轴标示，标题 2.将步长细化 题目：写出生成如图E2-1所示波形的MATLAB脚本M文件。图中虚线为正弦波，要求它的负半波被置零，且在处被削顶。 图 E2-1 题15所要表示的波形（逐段解析函数） 实现方法：plot画图 程序 t=0:0.01*pi:10; 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; 结果及结果分析 对源程序的改进：1.使用了常规的t生成方法。 2.因为原程序图形上下都顶到了，不好看。所以加上规定坐标轴值范围的语句 题目： 令，运行 [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) 的运行结果是什么 ? 实现方法：matlab矩阵操作 程序： clc clear 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) 的结果是0 S(1,1)/S(3,3) - cond(A) 的结果是0 S(1,1)*S(2,2)*S(3,3) - det(A) tol的结果是1 V(:,1)*null(A) 的结果是-1.3124e-016 abs(A*null(A)) tol 的结果是1；1；1 U(:,1:2) = = orth(A) 的结果是1，1 ；1，1 ；1，1 题目： 求积分，。提示：abs,sin,cos,cumsum,disp,plot,grid 实现方法：使用int函数 程