matlab实验步骤.doc

matlab实验步骤 篇一：matlab实验报告 软件课程设计 实验日期： 实验报告 模块名称：__MATLAB_____________ 专业班级：__ __________ 姓 名：__ _____________ 学 号：__ ____________ 年月日 实验报告日期： 年月 日 一、要求练习的实验部分 1. 题目：在时间区间 [0,10]中，绘出y?1?e?0.5tcos2t曲线。 实现方法： 在[0,10]区间中每隔0.1取一点，构造t数组，然后利用exp和cos函数计算出y，最后用plot函数画出图形。 程序： t=0:0.1:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,#39;r-#39;); shg 结果及结果分析： 2. 题目：写出生成如图E2-1所示波形的MATLAB脚本M文件。图中虚线为正弦波，要 求它的负半波被置零，且在 3 处被削顶。 2 图 E2-1题15所要表示的波形（逐段解析函数） 实现方法： 将t在区间[0,3*pi]用向量表示，分成500点，乘法要用.*，通过y=0和a=sin(pi/3)将函数y=sin(t)截断比较再相加得到f(t)； 对于ygt;0的部分进行截取，代码为z=(ygt;=0).*y； 对于ylt;a的部分进行截取，代码为z=(ygt;=a).*b+(ylt;a).*y 程序： t=linspace(0,3*pi,500); y=sin(t); a=sin(pi/3); z=(ygt;=0).*y; z=(ygt;=a).*a+(ylt;a).*z; plot(t,y,#39;:r#39;); hold on; plot(t,z,#39;b-#39;) xlabel(#39;t#39;),ylabel(#39;z=f(t)#39;),title(#39;逐段解析函数#39;) legend(#39;y=sin(t)#39;,#39;z=f(t)#39;,4),hold off 结果及结果分析： ?123???3. 题目：令A?456，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问????789?? 题： （1） sum(diag(S)gt;tol) - rank(A) 是多少 ？ （2） S(1,1) - norm(A) = 0 是多少 ? （3） sqrt(sum(diag(S*S))) - norm(A,#39;fro#39;) 的结果是什么 ? （4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ? （5） S(1,1)*S(2,2)*S(3,3) - det(A) lt; tol 的运行结果是什么？ （6） V(:,1)#39;*null(A) 得到什么结果 ? （7） abs(A*null(A)) lt; tol 得到什么结果 ? （8） 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.4082 S =16.8481 00 0 1.0684 0 0 00 V =-0.4797-0.7767-0.4082 -0.5724-0.0757 0.8165S(1,1) - norm(A) = 0-0.6651 0.6253-0.4082 键入 tol=S(1,1)*3*eps 得 tol =1.1223e-014 键入 sum(diag(S)gt;tol) - rank(A) 得 ans =0 键入 S(1,1) - norm(A) =0 得 ans =1 键入 sqrt(sum(diag(S*S))) - norm(A,#39;fro#39;) 得 ans =3.5527e-015 键入 S(1,1)/S(3,3) - cond(A) 得 ans =-8 键入 S(1,1)*S(2,2)*S(3,3) - det(A) lt; tol 得 ans=1 键入 V(:,1)#39;*null(A) 得 ans=0 键入 abs(A*null(A)) lt; tol 得 ans =1 1 1 键入 U(:,1:2)==orth(A) 得 ans =1 1 1 1 1 1 结果及结果分析： （1）sum(diag(S)gt;tol