【南邮】matlab软件的设计告报告.docVIP

软 件 设 计 报 告 课程名称 MATLAB通信仿真 实习时间 2011/2012学年第二学期 指导单位 ***学院 指导教师 *** 学生姓名 *** 班级学号 *** 学院(系) *** 专 业 *** 一、实验时间： 4月16日 - 4月27日 32学时 二、实验设备与仪器： 1、硬件：PC机 一台 2、软件：Matlab R2009b软件 三、实验基本要求： 了解MATLAB软件包，掌握相关的计算、设计方法，对设计的系统在Simulink环境中进行仿真。通过完成所有的练习题，初步熟悉Matlab命令窗的使用，M文件的编写，Simulink环境中仿真系统的搭建、运行、调试和验证。在此基础上，完成通信系统设计题，并调试验证通过。 四、实验内容：Matlab练习题 1、在时间区间 [0,10]中，绘出曲线。 1) M文件代码[No1code.m] clc; clear; t=0:0.1:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,r-); grid on; title(y=1-exp((-0.5)*t).*cos(2*t)的图形曲线); xlabel(x轴); ylabel(y轴); 2) 绘制曲线[No1pic.bmp] 2、写出生成如图E2-1所示波形的MATLAB脚本M文件。图中虚线为正弦波，要求它的负半波被置零，且在处被削顶。 1) M文件代码[No2code.m] clc; clear; t=linspace(0,3*pi,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),hold off 2) 绘制曲线[No2pic.bmp] 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) 的运行结果是什么 ? 1) M文件代码[No3code.m] clear; clc; disp(设 A=[1,2,3;4,5,6;7,8,9]，得) A=[1,2,3;4,5,6;7,8,9] %一个3*3矩阵 disp(设 [U,S,V]=svd(A)，得) [U,S,V]=svd(A) %奇异值分解(SVD)USV:U m*m酉矩阵 V n*n酉矩阵 S 对角阵 disp(设 tol=S(1,1)*3*eps，得) tol=S(1,1)*3*eps disp(设 a=sum(diag(S)tol) - rank(A)，求 a) a=sum(diag(S)tol) - rank(A) disp(设 b=S(1,1) - norm(A)，求 b) b=S(1,1) - norm(A) disp(设 c=sqrt(sum(diag(S*S))) - norm(A, fro )，求 c) c=sqrt(sum(diag(S*S))) - norm(A,fro) disp(设 d=S(1,1)/S(3,3) - cond(A)，求 d) d=S(1,1)/S(3,3) - cond(A) disp(设 e=S(1,1)*S(2,2)*S(3,3) - det(A) tol ，求 e) e=S(1,1)*S(2,2)*S(3,3) - det(A) tol disp(设 f=V(:,1) *null(A)，求 f) f=V(:,1)*null(A) disp(设 g=abs(A*null(A)) tol，求 g) g=abs(A*null(A)) tol disp(设 h=U(:,1:2) == orth(A)，求 h) h=U(:,1:2) == orth(A) 2) 实验结果[No3result.txt] 设 A=[1,2,3;4,5,6