matlab软件设计实验报告(各源程序均运行无误).docVIP

南京邮电大学通信学院 软件课程设计 实验报告 模块名称：__ MATLAB软件设计___ 专业班级：_____ 通信2班______ 姓 名：___ _ ______ 学 号：____ __ _____ 实验日期： 2012 年 3月5 日 实验报告日期： 2012年 3月8 日 在时间区间 [0,10]中，绘出曲线。 程序 close all clc clear t=0:0.1:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,r-); shg 运行结果： 写出生成如图E2-1所示波形的MATLAB脚本M文件。图中虚线为正弦波，要求它的负半波被置零，且在处被削顶。 程序： close all 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 运行结果： 令，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： 程序： close all 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) V(:,1)*null(A) U(:,1:2)=orth(A) sum(diag(S)tol) - rank(A) 是多少 ？ ans = 0 S(1,1) - norm(A) = =0 是多少 ? ans = 1 sqrt(sum(diag(S*S))) - norm(A,fro) 的结果是什么 ? ans = 0 S(1,1)/S(3,3) - cond(A) 的结果是什么 ? ans = 0 S(1,1)*S(2,2)*S(3,3) - det(A) tol 的运行结果是什么？ ans = 1 V(:,1)*null(A) 得到什么结果 ? ans = -5.5511e-017 abs(A*null(A)) tol 得到什么结果 ? ans = 1 1 1 U(:,1:2) = = orth(A) 的运行结果是什么 ? ans = 1 1 1 1 1 1 求积分，。 提示：abs,sin,cos,cumsum,disp,plot,grid 程序： close all clc clear i=0; for x=0:0.01:2*pi i=i+1; g=@(t)abs(sin(cos(t))); y(i)=quad(g,0,x); end plot(x,y,r) grid on 运行结果： 求方程的解。 提示：solve 程序： close all clc clear [x,y]=solve(x^2+y^2=1,x*y=2) 运行结果： x = ((15^(1/2)*i)/2 + 1/2)^(1/2)/2 - ((15^(1/2)*i)/2 + 1/2)^(3/2)/2 ((15^(1/2)*i)/2 + 1/2)^(3/2)/2 - ((15^(1/2)*i)/2 + 1/2)^(1/2)/2 (1/2 - (15^(1/2)*i)/2)^(1/2)/2 - (1/2 - (15^(1/2)*i)/2)^(3/2)/2 (1/2 - (15^(1/2)*i)/2)^(3/2)/2 - (1/2 - (15^(1/2)*i)/2)^(1/2)/2 y = ((15^(1/2)*i)/2 + 1/2)^(1/2) -(1/2*15^(1/2)*i + 1/2)^(1/2) (1/2 - (15^(1/2)*i)/2)^(1/2) -(1/2 - 1/2*15^(1/2)*i)^(1/2) 在某激励条件下，二阶系统归一化响应可表示为，其中为阻尼系数，，。请用不同的颜色或线型，在同一张图上，绘制取值下系统在区间内的响