2025年南京邮电大学Matlab软件设计深度解析教程.doc

南京邮电大学通信学院

软件課程设计

试验汇报

模块名称：___MATLAB软件设计

专业班级：__通信工程

姓名：________

学号：___

试验曰期：6月17—28曰试验汇报曰期：年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);

holdon;

plot(t,z,b-)

xlabel(t),ylabel(z=f(t)),title(逐段解析函数)

legend(y=sin(t),z=f(t),4),holdoff

成果：

令，运行[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:Um*m酉矩阵Vn*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)

成果：

设A=[1,2,3;4,5,6;7,8,9]，得

A=

123

456

789

设[U,S,V]=svd(A)，得

U=

-0.21480.88720.4082

-0.52060.2496-0.8165

-0.8263-0.38790.4082

S=

16.848100

01.06840

000.0000

V=

-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=