南京邮电大学通信学院软件课程设计实验报告模块名称：___MATLAB软件设计专业班级：__通信工程姓名：____ ____学号：___实验日期：2013年 6 月 17—28日实验报告日期： 2013 年 7 月 1 日一、要求练习的实验部分1．在时间区间 [0,10]中，绘出t e y t 2cos 15.0--=曲线。

程序： t=0:0.1:10;y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); shg 结果：2. 写出生成如图E2-1所示波形的MA TLAB 脚本M 文件。

图中虚线为正弦波，要求它的负半波被置零，且在23处被削顶。

程序：t=linspace(0,3*pi,500); y=sin(t); a=sin(pi/3); z=(y>=0).*y;z=(y>=a).*a+(y<a).*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 结果：3. 令⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A ，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： （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) < tol 的运行结果是什么？ （6） V(:,1)'*null(A) 得到什么结果 ?（7） abs(A*null(A)) < tol 得到什么结果 ? （8）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) < toldisp('设 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 34 5 67 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 00 1.0684 00 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 =设 b=S(1,1) - norm(A)，求 bb =设 c=sqrt(sum(diag(S*S))) - norm(A, fro )，求 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 ，求 ee =1设 f=V(:,1) *null(A)，求 f f =设 g=abs(A*null(A)) < tol ，求 g g =1 1 1设 h=U(:,1:2) == orth(A)，求 h h =1 1 1 1 1 1 结果分析：由上面的结果得出如下的题目答案：（1）sum(diag(S)>tol) - rank(A) 的结果是0； （2）S(1,1) - norm(A) = 0 的结果是0；（3）sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是3.5527e-015； （4）S(1,1)/S(3,3) - cond(A) 的结果是-8；（5）S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是1； （6）V(:,1)'*null(A)的结果是0； （7）abs(A*null(A)) < tol 的结果是1 1 1 ；（8）U(:,1:2)==orth(A)的运行结果是1 1 1 1 1 1 4. 求积分dt t x⎰]sin[cos ，π20≤≤x 。

程序：t=0:pi/100:2*pi; f=abs(sin(cos(t))); F=cumsum(f)*pi/100; plot(t,F);xlabel('x 轴'); ylabel('y 轴');title('描绘y=abs(sin(cos(t)))在(0,x)间积分曲线[x ∈(0,2π)]') grid on; 结果：5. 求方程2,122==+xy y x 的解。

程序： clear; clc;disp('计算方程组 x^2+y^2=1 x*y=2 的根 x y') [x,y] = solve('x^2+y^2 = 1','x*y = 2') 结果：计算方程组 x^2+y^2=1 x*y=2 的根 x y x =[ -1/2*(1/2*5^(1/2)+1/2*i*3^(1/2))^3+1/4*5^(1/2)+1/4*i*3^(1/2)] [ -1/2*(1/2*5^(1/2)-1/2*i*3^(1/2))^3+1/4*5^(1/2)-1/4*i*3^(1/2)] [ -1/2*(-1/2*5^(1/2)+1/2*i*3^(1/2))^3-1/4*5^(1/2)+1/4*i*3^(1/2)][ -1/2*(-1/2*5^(1/2)-1/2*i*3^(1/2))^3-1/4*5^(1/2)-1/4*i*3^(1/2)] y =[ 1/2*5^(1/2)+1/2*i*3^(1/2)] [ 1/2*5^(1/2)-1/2*i*3^(1/2)] [ -1/2*5^(1/2)+1/2*i*3^(1/2)] [ -1/2*5^(1/2)-1/2*i*3^(1/2)]6. 在某激励条件下，二阶系统归一化响应可表示为)sin(11)(θββξ+-=-t e t y t ，其中ξ为阻尼系数，21ξβ-=，)1(2ξξθ-=arctg 。

请用不同的颜色或线型，在同一张图上，绘制0.2:2.0:2.0=ξ取值下系统在]18,0[∈t 区间内的响应曲线，并要求用2.0=ξ和0.2=ξ对它们相应的两条曲线进行醒目的文字标志。

程序：b=0.2:0.2:2.0; t=0:0.1:18;color=['r','g.','bo','cx','m+','y*','kv','r--','g:','b-.']; for i=1:10p(i)=sqrt(1-b(i).*b(i));q(i)=atan(sqrt(1-b(i).*b(i))./b(i));y=1-1./p(i).*exp(-b(i).*t).*sin(p(i).*t+q(i)); if i==1plot(t,y,color(i))text(t(30),y(30),' \leftarrow ξ=0.2','FontSize',18) hold onelse if i==10plot(t,y,color(i))text(t(50),y(50),' \leftarrow ξ=2.0','FontSize',18) hold on elseplot(t,y,color(i)) hold on end end end 结果：7.构建一个简单的全波整流模型，并用示波器分两路同时观察原信号和整流后的信号波形。

要求：信源模块被重新命名为“输入正弦波”；信宿模块被重新命名为“示波器”；连接到信宿的信号线上分别标注“原信号”和“整流信号”；在模型窗中添加注释文本。

(1)、原理图：各模块参量设置：[示波器]ParametersNumber of axes: 2其他默认[其余模块]默认选项结果(2)、原理图：各模块参数设定[Gain]Gain: -1其他默认[其余模块]默认选项结果8.利用SIMULINK及其标准模块设计一个低通滤波器，从受噪声干扰的多频率混合信号π，而)(tx中获取10Hz的信号。