南邮软件课程设计MATLAB软件设计

南京邮电大学通信学院

软件课程设计

实验报告

模块名称：___MATLAB软件设计

专业班级：__通信工程_

姓名：______ ____

学号：___ _09

实验日期：2012 年 6 月 11—21日实验报告日期： 2012 年 6 月 26 日

索引目录

一、要求练习的实验部分 (1)

1、练习题1 (1)

2、练习题2 (1)

3、练习题3 (2)

4、练习题4 (5)

5、练习题5 (6)

6、练习题6 (7)

7、练习题7 (8)

8、练习题8 (11)

9、练习题9 (13)

10、练习题10 (14)

11、练习题11 (17)

12、练习题12 (20)

二、模拟数字电路仿真实验 (21)

(一)数字逻辑电路基础 (21)

(二)组合逻辑电路仿真 (22)

(三)时序逻辑电路仿真 (26)

三、数字信号处理仿真实验 (28)

（一）、利用Kaiser窗函数设计的FIR低通滤波器进行数字滤波 (28)

（二）、DSB-SC信号的生成与解调 (35)

四、通信系统实验 (40)

模拟调制---DSB调制与解调 (40)

一、要求练习的实验部分

1．在时间区间 [0,10]中，绘出t e y t 2cos 15.0--=曲线。

程序：（zoushuCode1.m ）

t=0:0.1:10;

y=1-exp((-0.5)*t).*cos(2*t);

plot(t,y,'r-');

shg

结果：(图图1.fig)

2. 写出生成如图E2-1所示波形的MA TLAB 脚本M 文件。图中虚线为正弦波，要求它的负半波被置零，且在2

3处被削顶。 程序：（zoushuCode2.m ）

t=linspace(0,3*pi,500);

y=sin(t);

a=sin(pi/3);

z=(y>=0).*y;

z=(y>=a).*a+(y

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

结果：（图图2.fig ）

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) 的运行结果是什么 ?

程序：(zoushuCode3.m)

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')

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)

结果：（结果3.txt）

设 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