南京邮电大学通信学院
实验日期:
软件课程设计 实验报告
模块名称:_____MATLAB_________ 专业班级:_____ ________ 姓 名:_____ ________ 学 号:____ _______
年 月 日 实验报告日期: 2015 年月 日 2015
设计题目 1.组合逻辑电路仿真 2. DSB-SC信号的生成与解调 3. 模拟调制---DSB调制与解调 任务要求 1. 掌握MATLAB部分工具包的使用方法,侧重于simulink环境。 2.掌握DSP设计、通信系统性能设计、模拟和数字电路设计中的基本设计方法。 实验设备微型计算机,MATLAB2012 及软件 同组人员 学号及姓无 名 参考文献 1. 2.
一、要求练习的实验部分
1.在时间区间 [0,10]中,绘出y?1?e?0.5tcos2t曲线。 实现方法 程序
t=0:0.1:10;
y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); shg
结果及结果分析
1.510.50012345678910
2. 写出生成如图E2-1所示波形的MATLAB脚本M文件。图中虚线为正弦波,要求它的
负半波被置零,且在
3处被削顶。 2程序:
t=linspace(0,3*pi,500); y=sin(t); a=sin(pi/3); z=(y>=0).*y;
z=(y>=a).*a+(y
xlabel('t'),ylabel('z=f(t)'),title('逐段解析函数') legend('y=sin(t)','z=f(t)',4),hold off
结果:
逐段解析函数10.80.60.40.2 z=f(t)0-0.2-0.4-0.6-0.8-1 012345t678y=sin(t)z=f(t)910
?123???3. 令A?456,运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题: ????789??(1)
(2) (3) (4) (5) (6) (7) (8)
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:U m*m酉矩阵 V n*n酉矩阵 S 对角阵
disp('设 tol=S(1,1)*3*eps,得') tol=S(1,1)*3*eps
disp('设 a=sum(diag(S)>tol) - rank(A),求 a')%disg(s)取s主对角线元素
a=sum(diag(S)>tol) - rank(A)% rank(A)函数返回结果集分区内指定字段的值的排名,指定字段的值的排名是相关行之前的排名加一
disp('设 b=S(1,1) - norm(A),求 b')%norm(A)求矩阵A的范数 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 =
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 0 0 0 1.0684 0 0 0 0.0000 V =
-0.4797 -0.7767 -0.4082 -0.5724 -0.0757 0.8165