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

《南京邮电大学matlab软件设计doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、南京邮电大学通信学院软件课程设计实验报告模块名称：___MATLAB软件设计专业班级：__通信工程_17班姓名：____张步涛______学号：____实验日期：2013年6月17—28日实验报告日期：2013年7月1日-47-一、要求练习的实验部分1．在时间区间[0,10]中，绘出曲线。程序：t=0:0.1:10;y=1-exp((-0.5)*t).*cos(2*t);plot(t,y,'r-');shg结果：2.写出生成如图E2-1所示波形的MATLAB脚本M文件。图中虚线为正弦波，要求它的负半波

2、被置零，且在处被削顶。程序：t=linspace(0,3*pi,500);y=sin(t);a=sin(pi/3);z=(y>=0).*y;z=(y>=a).*a+(y

3、然后回答以下问题：（1）sum(diag(S)>tol)-rank(A)是多少？（2）S(1,1)-norm(A)=0是多少?（3）sqrt(sum(diag(S*S)))-norm(A,'fro')的结果是什么?（4）S(1,1)/S(3,3)-cond(A)的结果是什么?（5）S(1,1)*S(2,2)*S(3,3)-det(A)

4、)的运行结果是什么?程序：clear;clc;disp('设A=[1,2,3;4,5,6;7,8,9]，得')-47-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:Um*m酉矩阵Vn*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

5、)>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)

6、(2,2)*S(3,3)-det(A)

7、0.2496-0.8165-0.8263-0.38790.4082S=16.84810001.06840-47-000.0000V=-0.4797-0.7767-0.4082-0.5724-0.07570.8165-0.66510.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)-norm(A)，求bb=0设c=sqrt(sum(diag(S*S)))-norm(A,fro)，

8、求cc=3.5527e-015设d=S(1,1)/S(3,3)-cond(A)，求dd=-8设e=S(1,1)*S(2,2)*S(3,3)-det(A)tol)-rank(A)的结果是0；（2