南京邮电大matlab软件设计.doc.doc

PAGE 1 -南京邮电大学通信学院软件课程设计实验报告模块名称：___MATLAB软件设计 专业班级：__通信工程_17班 姓 名：____张步涛__ ____ 学 号：___ 实验日期：2013年 6 月 17—28日 实验报告日期： 2013 年 7 月 1 日一、要求练习的实验部分在时间区间 [0,10]中，绘出曲线。程序：t=0:0.1:10;y=1-exp((-0.5)*t).*cos(2*t);plot(t,y,r-);shg结果：写出生成如图E2-1所示波形的MATLAB脚本M文件。图中虚线为正弦波，要求它的负半波被置零，且在处被削顶。程序：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结果：令，运行 [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) 的运行结果是什么 ?程序：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*epsdisp(设 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)) toldisp(设 h=U(:,1:2) == orth(A)，求 h)h=U(:,1:2) == orth(A)结果：设 A=[1,2,3;4,5,6;7,8,9]，得A = 1 2 3 4 5 6 7 8 9设 [U,S,V]=svd(A)，得U = -0.2148 0.8872 0.4082 -0.5206 0.2496 -0.8165 -0.8263 -0.3879 0.4082S = 16.8481 0 0 0 1.0684 0 0 0 0.0000V = -0.4797 -0.7767 -0.4082 -0.5724 -0.0757 0.8165 -0.6651 0.6253 -0.4082设 tol=S(1,1)*3*eps，得tol = 1.1223e-014设 a=sum(diag(S)tol) - rank(A)，求 aa = 0设 b=S(1,1