1、y=sin(t),4),hold off3. 令,运行 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*n
2、ull(A) tol) - rank(A),求 aa=sum(diag(S)tol) - rank(A)设 b=S(1,1) - norm(A),求 bb=S(1,1) - norm(A)设 c=sqrt(sum(diag(S*S) - norm(A, fro ),求 cc=sqrt(sum(diag(S*S) - norm(A,设 d=S(1,1)/S(3,3) - cond(A),求 dd=S(1,1)/S(3,3) - cond(A)设 e=S(1,1)*S(2,2)*S(3,3) - det(A) tol ,求 ee=S(1,1)*S(2,2)*S(3,3) - det(A) tol
3、设 f=V(:,1) *null(A),求 ff=V(:*null(A)设 g=abs(A*null(A) tol,求 gg=abs(A*null(A) 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的结果是
4、1 1 ;(8)U(:2)=orth(A)的运行结果是1 14. 求积分,。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)grid on;5. 求方程的解。计算方程组 x2+y2=1 x*y=2 的根 x yx,y = solve(x2+y2 = 1x*y = 2计算方程组 x2+y2=1 x*y=2 的根 x yx = -1/2*(1/2*5(1/2)+1/2*i*3(1/2)3+1/4*5(1/2)+1/4*i*3(1
5、/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. 在某激励条件下,二阶系统归一化响应可表示为,其中为阻
6、尼系数,请用不同的颜色或线型,在同一张图上,绘制取值下系统在区间内的响应曲线,并要求用和对它们相应的两条曲线进行醒目的文字标志。b=0.2:0.2:2.0;18;color=g.bocxm+y*kvr-g:b-.;for i=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); if i=1 plot(t,y,color(i) text(t(30),y(30), leftarrow =0.2FontSize,18) hold on
7、 else if i=10 text(t(50),y(50), leftarrow =2.0 else endend7. 构建一个简单的全波整流模型,并用示波器分两路同时观察原信号和整流后的信号波形。要求:信源模块被重新命名为“输入正弦波”;信宿模块被重新命名为“示波器”;连接到信宿的信号线上分别标注“原信号”和“整流信号”;在模型窗中添加注释文本。(1)、原理图:各模块参量设置: 示波器ParametersNumber of axes: 2其他默认其余模块默认选项结果(2)、原理图:各模块参数设定GainGain: -1 其余模块8. 利用 SIMULINK及其标准模块设计一个低通滤波器,从
8、受噪声干扰的多频率混合信号中获取10Hz的信号。在此,而Sine WaveFrequency (rad/sec): 10Sine Wave1 100Phase (rad): pi/2Random NumberVariance: 0.2Sample time:SumList of signs: |+Analog Filter DesignPassband edge frequency (rad/s):9. 已知系统的状态方程为:,请构建该系统的仿真模型,并用XY Graph模块观察相轨迹。原理图:参数配置; |-+ integratorInitial condition:0.2integrato
9、r1 Initial condition:10. CD74HC00仿真(数字电路)实现方法:这里就是四个与非门的实现,利用Logical Operator模块做出四个与非门,还要再和Vcc和Gnd相与,再输出。然后用Subsystem封装成子系统即可;(1)、封装前的原理图(2)、分装后的测试原理图分析:前两路分别为产生“1100”和“0101”的脉冲发生器的波形,后四路波形则分别为四个与非门的输出;可以看到,四路输出信号与输入信号均符合与非逻辑的关系,符合7400的真值表,设计正确!112FSK信号的生成与接收滤波1) 产生二进制0,1随机信号波,信号波特率为1000B2) 对此信号进行2F
10、SK调制,f1=1500Hz, f2=3050Hz3) 加入高斯白噪声4) 在接收端进行带通滤波5) 抽样速率为20KHz.6) 观察滤波前后信号的波形和频谱。利用Bernoulli Binary Generator模块产生随机信号;利用Sine Wave模块产生f1=1500Hz, f2=3050Hz 的载波,结合Switch模块进行2FSK调制;利用AWGN Channel模块构建高斯白噪声信道,其参数设置为:信噪比:10dB,输入信号功率:1W;利用Digital Filter Design设计带通滤波器,分别为1500Hz和3050Hz滤波然后利用相干解调,在将两个载波相加,用符号函数
11、判决出原信号。(zoushuSimulink11.mdl、zoushuSimulink11b.mdl)参数设计: Sine Wave: Frequency:1500 Hz 其他:默认 Sine Wave 1:3050 Hz Bernoulli Binary Generator Sample time:1/1000 Switch: Threshoid:1/2 AWGN Channel SNR:10 db Power:1W Digital Filter Design Bandpass; Fs:20kHz Fstop1:0.8KHz Fpass1:1.2KHz: Fstop2:2KHz Fpass:
12、2.2.5KHz: Digital Filter Design 1 2.5KHz Fpass2:2.8KHz:3.5KHz4KHz:调制的原理图:调制加解调的原理图:频谱图:调制、解调的波形图:原信号经过2FSK调制,经过噪声信道,再进行滤波、相干解调,最后又恢复出原信号,结果符合题目要求,系统设计正确!滤波器等相关设备的参数也设置恰当!实验成功!12. 创建一个简单的离散多速率系统:单位阶跃信号经过具有不同速率的采样后分别用作两个离散传递函数的输入。这两个离散传递函数有相同的有理分式,但采样时间和时间偏置二元对分别设为1 , 0.1和0.7 , 0。观察这两个离散传递函数的输出有什么不同;用
13、不同的颜色标帜不同采样速率系统。(通信系统)二、模拟数字电路仿真实验(一) 数字逻辑电路基础熟悉常用逻辑单元的特性,学会运用基本逻辑单元(与、或、非、异或、R-S触发器、D触发器、J-K触发器等),修改参数、增减输入/输出端。(二) 组合逻辑电路仿真1、 设计二/四线译码器利用Logic Operator模块实现非门和与非门;利用Pulse Generator模块产生0101和0011序列码,用来检验设计出的子系统。 原理图:封装前原理图:结果及分析 分析:当!EN端置0时,此时,随着第二、三两路信号输入的不同,作为输出端的Y0,Y1,Y2,Y3依次输出低电平,符合真值表要求,证明设计的子系统
14、工作正常。EN端置1时,子系统没有被使能。此时,作为输出端的Y0,Y1,Y2,Y3没有随着输入A0,A1作出相应变化,而是一直输出高电平,证明设计的子系统工作正常。设计正确!2、 设计四选一数据选择器能仿真测试,并设计成子系统元件 实现方法:利用的模块同上;地址信号用Pulse Generator产生;D0到D3也由Pulse Generator模块产生,为了区分各路信号,将其信号周期依次设置为0.0625s, 0.125s, 0.25s和0.5秒。原理图结果及分析:前两路分别为输入信号,组合产生00,01,10,11的地址信号;第四到第七这四路信号为D0到D1,它们的频率各不相同;最后一路为
15、输出信号;可以看出,随着输入的地址信号的变化,输入一次输出该地址信号对应的信号。子系统设计成功!(三) 时序逻辑电路仿真1、 设计四位二进制计数器(带置位和清零) 利用JK触发器和逻辑门构成之。利用比较系数的方法可得电路的驱动方程为检验电路将P, T, LD, CR置为1,然后观察示波器输出。封装前的原理图:封装后:从图中可以看出,脉冲在上升沿有效,计数器按00001111方式进行计数。当P、T、L(置数)、CR(清零)有一个为低电平时的输出为全零。P、T为0时,计数器不能工作;L为0置数,D3D0全零,输出也全零。CR为零时清零,所以输出也全零。三、数字信号处理仿真实验(一)、利用Kaise
16、r窗函数设计的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)设计基本思路:首先选择一个符合要求的理想滤波器(这里是理想低通滤波器);由于理想滤波器的冲击响应是非
17、因果的且无限长,为了能用FIR滤波器实现,必须用适当的窗函数来截取,从而得到线性相位和因果的FIR滤波器。而滤波器的特性与窗函数的长度(也就决定了FIR滤波器的长度)密切相关。 一个截止频率为wc的理想低通滤波器的单位冲击响应为:h(n)=sinwc(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窗提供了最大的主瓣宽度,从而提
18、供最陡的过渡带。窗函数的表达式为:w(n)=I0square(1-(1-2n/M)2)其中I0为零阶Bessel函数。Kaiser窗的优点在于它可以通过改变参数和M来改变阻带衰减和过渡带宽。在MATLAB的SimLink中,有相关的模块供调用。题目:1 生成2KHz和8KHz的混合信号,使该信号通过上述滤波器,观察滤波前后的波形变化和频谱分布情况。(基本要求)首先用两个Sine Wave模块分别产生2kHz和8kHz的正弦信号,之后用add模块将二者混合。用Digital Filter Design设计低通滤波器,具体设置如图所示:、结果、图为滤波之后和之前的信号频谱,可以观察到,滤波之后,处
19、于通带之外的8kHz的正弦分量消失了。上图是示波器的波形,其中第一路为2kHz和8kHz的混合信号,第二路为滤波之后的信号,可见低通滤波器符合要求!2 在改信号中加入高斯白噪声,观察滤波前后的波形变化和频谱分布情况。在上题基础上增加AWGN信道,加入高斯白噪声。上图是加噪声后示波器的波形,其中第一路为2kHz和8kHz和噪声的混合信号,第二路为滤波之后的信号(还有低频的噪声)!3 、将上述滤波器改为中心频率为100KHz的带通滤波器,信号源为带外的90KHz和带内的100KHz的混合正弦信号,重做1,2问。利用Digital Filter Design设计带通滤波器,具体设置为:Fs=1000
20、kHz, Fstop1=99 kHz, Fpass1=99.5 kHz, Fpass2=100.5 kHz, Fstop2=101 kHz。无噪声的频谱分析:无噪声的波形分析:加入噪声的频谱分析:加入噪声的波形分析:(二)、 DSB-SC信号的生成与解调1) 用离散(DSP)的方法生成DSB信号2) 载波频率为150KHz,音频为500Hz和2000Hz的混合音。4) 带通滤波,分别用正交解调和检波的方法实现解调。5) 抽样速率为600Hz。6) 观察各点的信号波形和频谱。A.正交解调法实现方法本题目要求用离散的方法生成DSB信号,因此我们用编写M文件方式而不是用Simulink来实现。生成正
21、弦音频信号时,对时间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*105,6*105);%调制subplot(3,1,2
22、),plot(t,z);DSB-SC信号图y=awgn(z,5); %加入高斯白噪声subplot(3,1,3),plot(t,y);加入噪声后的DSB-SC图,SNR=5figure(2);f=(0:length(fft(z)-1)*(600*1000-1)/length(fft(z);plot(f,abs(fft(z);调制后的频谱图%带通滤波器参数As=100; Ap=1;fs=6*105; %采样频率 fp1=400; fp2=1.6*105; %通频带fs1=300; fs2=1.7*105; %截止带wp=fp1 fp2/(fs/2);ws=fs1 fs2/(fs/2);n,wn=cheb1ord(wp,ws,Ap,As);b=fir1(
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1