南京邮电大学matlab软件设计超详细doc综述Word下载.docx
《南京邮电大学matlab软件设计超详细doc综述Word下载.docx》由会员分享,可在线阅读,更多相关《南京邮电大学matlab软件设计超详细doc综述Word下载.docx(49页珍藏版)》请在冰豆网上搜索。
y=sin(t)'
'
4),holdoff
3.令
,运行[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;
7,8,9]%一个3*3矩阵
设[U,S,V]=svd(A),得'
[U,S,V]=svd(A)%奇异值分解(SVD)USV:
Um*m酉矩阵Vn*n酉矩阵S对角阵
设tol=S(1,1)*3*eps,得'
tol=S(1,1)*3*eps
设a=sum(diag(S)>
tol)-rank(A),求a'
a=sum(diag(S)>
tol)-rank(A)
设b=S(1,1)-norm(A),求b'
b=S(1,1)-norm(A)
设c=sqrt(sum(diag(S*S)))-norm(A,fro),求c'
c=sqrt(sum(diag(S*S)))-norm(A,'
设d=S(1,1)/S(3,3)-cond(A),求d'
d=S(1,1)/S(3,3)-cond(A)
设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
设f=V(:
1)*null(A),求f'
f=V(:
*null(A)
设g=abs(A*null(A))<
tol,求g'
g=abs(A*null(A))<
tol
设h=U(:
2)==orth(A),求h'
h=U(:
2)==orth(A)
7,8,9],得
A=
123
456
789
设[U,S,V]=svd(A),得
U=
-0.21480.88720.4082
-0.52060.2496-0.8165
-0.8263-0.38790.4082
S=
16.848100
01.06840
000.0000
V=
-0.4797-0.7767-0.4082
-0.5724-0.07570.8165
-0.66510.6253-0.4082
设tol=S(1,1)*3*eps,得
tol=
1.1223e-014
tol)-rank(A),求a
a=
0
设b=S(1,1)-norm(A),求b
b=
设c=sqrt(sum(diag(S*S)))-norm(A,fro),求c
c=
3.5527e-015
设d=S(1,1)/S(3,3)-cond(A),求d
d=
-8
tol,求e
e=
1
1)*null(A),求f
f=
tol,求g
g=
2)==orth(A),求h
h=
11
结果分析:
由上面的结果得出如下的题目答案:
(1)sum(diag(S)>
tol)-rank(A)的结果是0;
(2)S(1,1)-norm(A)=0的结果是0;
(3)sqrt(sum(diag(S*S)))-norm(A,'
)的结果是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(:
*null(A)的结果是0;
(7)abs(A*null(A))<
tol的结果是1
1;
(8)U(:
2)==orth(A)的运行结果是11
4.求积分
,
。
pi/100:
2*pi;
f=abs(sin(cos(t)));
F=cumsum(f)*pi/100;
plot(t,F);
x轴'
ylabel('
y轴'
title('
描绘y=abs(sin(cos(t)))在(0,x)间积分曲线[x∈(0,2π)]'
gridon;
5.求方程
的解。
计算方程组x^2+y^2=1x*y=2的根xy'
[x,y]=solve('
x^2+y^2=1'
x*y=2'
计算方程组x^2+y^2=1x*y=2的根xy
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.在某激励条件下,二阶系统归一化响应可表示为
,其中
为阻尼系数,
请用不同的颜色或线型,在同一张图上,绘制
取值下系统在
区间内的响应曲线,并要求用
和
对它们相应的两条曲线进行醒目的文字标志。
b=0.2:
0.2:
2.0;
18;
color=['
g.'
bo'
cx'
m+'
y*'
kv'
r--'
g:
'
b-.'
];
fori=1:
10
p(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));
ifi==1
plot(t,y,color(i))
text(t(30),y(30),'
\leftarrowξ=0.2'
FontSize'
18)
holdon
elseifi==10
text(t(50),y(50),'
\leftarrowξ=2.0'
else
end
end
7.构建一个简单的全波整流模型,并用示波器分两路同时观察原信号和整流后的信号波形。
要求:
信源模块被重新命名为“输入正弦波”;
信宿模块被重新命名为“示波器”;
连接到信宿的信号线上分别标注“原信号”和“整流信号”;
在模型窗中添加注释文本。
(1)、原理图:
各模块参量设置:
[示波器]
Parameters
Numberofaxes:
2
其他默认
[其余模块]
默认选项
结果
(2)、原理图:
各模块参数设定
[Gain]
Gain:
-1
[其余模块]
8.利用SIMULINK及其标准模块设计一个低通滤波器,从受噪声干扰的多频率混合信号
中获取10Hz的信号。
在此
,而
[SineWave]
Frequency(rad/sec):
10
[SineWave1]
100
Phase(rad):
pi/2
[RandomNumber]
Variance:
0.2
Sampletime:
[Sum]
Listofsigns:
|+++
[AnalogFilterDesign]
Passbandedgefrequency(rad/s):
9.已知系统的状态方程为:
,请构建该系统的仿真模型,并用XYGraph模块观察
相轨迹。
原理图:
参数配置;
|--+
[integrator]
Initialcondition:
0.2
[integrator1]
Initialcondition:
10.CD74HC00仿真(数字电路)
实现方法:
这里就是四个与非门的实现,利用LogicalOperator模块做出四个与非门,还要再和Vcc和Gnd相与,再输出。
然后用Subsystem封装成子系统即可;
(1)、封装前的原理图
(2)、分装后的测试原理图
分析:
前两路分别为产生“1100”和“0101”的脉冲发生器的波形,后四路波形则分别为四个与非门的输出;
可以看到,四路输出信号与输入信号均符合与非逻辑的关系,符合7400的真值表,设计正确!
11.2FSK信号的生成与接收滤波
1)产生二进制[0,1]随机信号波,信号波特率为1000B
2)对此信号进行2FSK调制,f1=1500Hz,f2=3050Hz
3)加入高斯白噪声
4)在接收端进行带通滤波
5)抽样速率为20KHz.
6)观察滤波前后信号的波形和频谱。
利用BernoulliBinaryGenerator模块产生随机信号;
利用SineWave模块产生f1=1500Hz,f2=3050Hz的载波,结合Switch模块进行2FSK调制;
利用AWGNChannel模块构建高斯白噪声信道,其参数设置为:
信噪比:
10dB,输入信号功率:
1W;
利用DigitalFilterDesign设计带通滤波器,分别为1500Hz和3050Hz滤波
然后利用相干解调,在将两个载波相加,用符号函数判决出原信号。
(zoushuSimulink11.mdl、zoushuSimulink11b.mdl)
参数设计:
[SineWave]:
Frequency:
1500Hz
其他:
默认
[SineWave1]:
3050Hz
[BernoulliBinaryGenerator]
Sampletime:
1/1000
[Switch]:
Threshoid:
1/2
[AWGNChannel]
SNR:
10db
Power:
1W
[DigitalFilterDesign]
Bandpass;
Fs:
20kHz
Fstop1:
0.8KHz
Fpass1:
1.2KHz:
Fstop2:
2KHz
Fpass:
2.2.5KHz:
[DigitalFilterDesign1]
2.5KHz
Fpass2:
2.8KHz:
3.5KHz
4KHz:
调制的原理图:
调制加解调的原理图:
频谱图:
调制、解调的波形图:
原信号经过2FSK调制,经过噪声信道,再进行滤波、相干解调,最后又恢复出原信号,结果符合题目要求,系统设计正确!
滤波器等相关设备的参数也设置恰当!
实验成功!
12.创建一个简单的离散多速率系统:
单位阶跃信号经过具有不同速率的采样后分别用作两个离散传递函数的输入。
这两个离散传递函数有相同的有理分式
,但采样时间和时间偏置二元对分别设为[1,0.1]和[0.7,0]。
观察这两个离散传递函数的输出有什么不同;
用不同的颜色标帜不同采样速率系统。
(通信系统)
二、模拟数字电路仿真实验
(一)数字逻辑电路基础
熟悉常用逻辑单元的特性,学会运用基本逻辑单元(与、或、非、异或、R-S触发器、D触发器、J-K触发器等),修改参数、增减输入/输出端。
(二)组合逻辑电路仿真
1、设计二/四线译码器
利用LogicOperator模块实现非门和与非门;
利用PulseGenerator模块产生0101和0011序列码,用来检验设计出的子系统。
原理图:
封装前原理图:
结果及分析
分析:
当!
EN端置0时,此时,随着第二、三两路信号输入的不同,作为输出端的Y0,Y1,Y2,Y3依次输出低电平,符合真值表要求,证明设计的子系统工作正常。
EN端置1时,子系统没有被使能。
此时,作为输出端的Y0,Y1,Y2,Y3没有随着输入A0,A1作出相应变化,而是一直输出高电平,证明设计的子系统工作正常。
设计正确!
2、设计四选一数据选择器
能仿真测试,并设计成子系统元件
实现方法:
利用的模块同上;
地址信号用PulseGenerator产生;
D0到D3也由PulseGenerator模块产生,为了区分各路信号,将其信号周期依次设置为0.0625s,0.125s,0.25s和0.5秒。
原理图
结果及分析:
前两路分别为输入信号,组合产生00,01,10,11的地址信号;
第四到第七这四路信号为D0到D1,它们的频率各不相同;
最后一路为输出信号;
可以看出,随着输入的地址信号的变化,输入一次输出该地址信号对应的信号。
子系统设计成功!
!
(三)时序逻辑电路仿真
1、设计四位二进制计数器(带置位和清零)
利用JK触发器和逻辑门构成之。
利用比较系数的方法可得电路的驱动方程为
检验电路将P,T,LD,CR置为1,然后观察示波器输出。
封装前的原理图:
封装后:
从图中可以看出,脉冲在上升沿有效,计数器按0000~1111方式进行计数。
当P、T、L(置数)、CR(清零)有一个为低电平时的输出为全零。
P、T为0时,计数器不能工作;
L为0置数,D3~D0全零,输出也全零。
CR为零时清零,所以输出也全零。
三、数字信号处理仿真实验
(一)、利用Kaiser窗函数设计的FIR低通滤波器进行数字滤波
利用Kaiser窗函数,设计具有如下指标的FIR低通滤波器:
fs=20KHz,fpass=4kHz,fstop=5KHz,Apass=0.1dB,Astop=80dB。
归一化低通滤波器的主要性能指标有:
绝对指标:
[0,wp]为滤波器的通带,δp为可容许的通带波动;
[ws,π]为滤波器的阻带,δs为阻带波动。
相对指标:
Ap为用dB表示的通带波动;
As为用dB表示的阻带波动
相对指标与绝对指标的关系为:
Ap=-20lg((1-δp)/(1+δp))
As=-20lg(δs)
设计基本思路:
首先选择一个符合要求的理想滤波器(这里是理想低通滤波器);
由于理想滤波器的冲击响应是非因果的且无限长,为了能用FIR滤波器实现,必须用适当的窗函数来截取,从而得到线性相位和因果的FIR滤波器。
而滤波器的特性与窗函数的长度(也就决定了FIR滤波器的长度)密切相关。
一个截止频率为wc的理想低通滤波器的单位冲击响应为:
h(n)=sin[wc(n-a)]/[π(n-a)]
a为采样延迟,对应的频谱特性为:
H(w)=exp(-jaw),当|w|<
wc时。
用一个窗函数w(n)去截取后,得到:
hd(n)=h(n).w(n)
窗函数满足在[0,M]内关于a对称,在其他区域为0。
在频域,FIR滤波器的频率响应由理想滤波器的频响和窗函数的频响卷积得到。
对给定的阻带衰减,Kaiser窗提供了最大的主瓣宽度,从而提供最陡的过渡带。
窗函数的表达式为:
w(n)=I0{βsquare(1-(1-2n/M)2)}
其中I0为零阶Bessel函数。
Kaiser窗的优点在于它可以通过改变参数β和M来改变阻带衰减和过渡带宽。
在MATLAB的SimLink中,有相关的模块供调用。
题目:
1生成2KHz和8KHz的混合信号,使该信号通过上述滤波器,观察滤波前后的波形变化和频谱分布情况。
(基本要求)
首先用两个SineWave模块分别产生2kHz和8kHz的正弦信号,之后用add模块将二者混合。
用DigitalFilterDesign设计低通滤波器,具体设置如图所示:
、
结果、
图为滤波之后和之前的信号频谱,可以观察到,滤波之后,处于通带之外的8kHz的正弦分量消失了。
上图是示波器的波形,其中第一路为2kHz和8kHz的混合信号,第二路为滤波之后的信号,可见低通滤波器符合要求!
2在改信号中加入高斯白噪声,观察滤波前后的波形变化和频谱分布情况。
在上题基础上增加AWGN信道,加入高斯白噪声。
上图是加噪声后示波器的波形,其中第一路为2kHz和8kHz和噪声的混合信号,第二路为滤波之后的信号(还有低频的噪声)!
3、将上述滤波器改为中心频率为100KHz的带通滤波器,信号源为带外的90KHz和带内的100KHz的混合正弦信号,重做1,2问。
利用DigitalFilterDesign设计带通滤波器,具体设置为:
Fs=1000kHz,Fstop1=99kHz,Fpass1=99.5kHz,Fpass2=100.5kHz,Fstop2=101kHz。
无噪声的频谱分析:
无噪声的波形分析:
加入噪声的频谱分析:
加入噪声的波形分析:
(二)、DSB-SC信号的生成与解调
1)用离散(DSP)的方法生成DSB信号
2)载波频率为150KHz,音频为500Hz和2000Hz的混合音。
4)带通滤波,分别用正交解调和检波的方法实现解调。
5)抽样速率为600Hz。
6)观察各点的信号波形和频谱。
A.正交解调法
实现方法
本题目要求用离散的方法生成DSB信号,因此我们用编写M文件方式而不是用Simulink来实现。
生成正弦音频信号时,对时间t用向量表示,可以实现离散化。
DSB调制与解调分别使用amod,ademod函数。
加入高斯白噪声用awgn函数。
带通滤波用cheb1ord和fftfilt函数实现。
画频谱用fft函数实现。
程序代码(zoushuCode3002.m):
1/(600*1000-1):
0.01;
m1=sin(500*2*pi*t);
%500Hz正弦信号
m2=sin(2*pi*2000*t);
%2000Hz正弦信号
m=m1+m2;
%混合信号
figure
(1);
subplot(3,1,1);
plot(t,m);
原始信号'
z=amod(m,1.5*10^5,6*10^5);
%调制
subplot(3,1,2),plot(t,z);
DSB-SC信号图'
y=awgn(z,5);
%加入高斯白噪声
subplot(3,1,3),plot(t,y);
加入噪声后的DSB-SC图,SNR=5'
figure
(2);
f=(0:
length(fft(z))-1)*(600*1000-1)/length(fft(z));
plot(f,abs(fft(z)));
调制后的频谱图'
%带通滤波器参数
As=100;
Ap=1;
fs=6*10^5;
%采样频率
fp1=400;
fp2=1.6*10^5;
%通频带
fs1=300;
fs2=1.7*10^5;
%截止带
wp=[fp1fp2]/(fs/2);
ws=[fs1fs2]/(fs/2);
[n,wn]=cheb1ord(wp,ws,Ap,As);
b=fir1(