信号与系统
仿真实验报告
1、求三阶系统
8106)
65(5)(2
3
2+++++=s s s s s s H 的单位阶跃响应，并绘制响应波形图。 解题分析：
本题求单位阶跃响应，输入传递函数，用step 函数求解，即可绘制出响应波形图。
程序设计如下： >> t=0:0.01:15;
>> num=[5 25 30];%分子系数 >> den=[1 6 10 8];%分母系数 >> sys=tf(num,den);%传递函数 >> g=step(sys,t);%求单位阶跃响应 >> plot(t,g);%绘制单位阶跃响应
响应的波形图：
051015
0.511.522.533.54
4.5
2、一个因果线性移不变系统)2()()2(81.0)(--+-=n x n x n y n y ， 求：（1）)(z H ； （2）冲激响应)(n h ； （3）单位阶跃响应)(n u ；
（4）
)(ω
j e H ，并绘出幅频和相频特性。 解题分析：
由差分方程可以直接得到系统函数)(z H ，求Z 反变换可得到单位冲激响应，还有单位阶跃响应，及幅频和相频特性。
程序设计如下： >> clear;
>> num=[1 0 -1];%分子系数 >> den=[1 0 -0.81];%分母系数 >> figure(1);
>> subplot(2,1,1);
>> dimpulse(num,den,20); %冲激响应 >> subplot(2,1,2);
>> dstep(num,den,50);%单位阶跃响应 >> w=[0:1:512]*pi/512; >>figure(2);
>>freqz(num,den,w);%幅频和相频特性
响应的波形图：
02468101214161820
-0.5
0.5
1
Impulse Response
Time (sec)
A m p l i t u d e
05101520253035404550
0.5
1
Step Response
Time (sec)
A m p l i t u d e
幅频和相频特性：
0.1
0.2
0.30.40.50.60.70.80.9
1
-100
-50050100
Normalized Frequency (?π rad/sample)
P h a s e (d e g r e e s )
0.1
0.2
0.30.40.50.60.70.80.9
1
-300
-200-1000100
Normalized Frequency (?π rad/sample)
M a g n i t u d e (d B )
实验三
3、模拟信号)8cos(5)4sin(2)(t t t x ππ+= ，求64=N 的DFT 的幅值谱和相位谱。
解题分析：
输入模拟信号，进行抽样，然后求其DFT ，从而输出幅值谱和相位谱。
程序设计如下： >> N=64;%采样点数 >>T=1;%采样时间终点
>>t=linspace(0,T,N);%给出N 个采样时间ti(i=1:N) >>x=2*sin(4*pi*t)+5*cos(8*pi*t);%求各采样点样本值x >>dt=t(2)-t(1);%采样周期 >>f=1/dt;%采样频率
>>X=fft(x);%计算x 的快速傅里叶变换X >>F=X(1:N/2+1); >>f=f*(0:N/2)/N; >> figure(1); >> plot(t,x); >>figure(2);
>>plot(f,abs(F),'-*');%绘制幅值-频率图 >>figure(3);
>>plot(f,angle(F),'-*');%绘制相位-频率图 模拟信号：
0.10.20.30.40.50.60.70.80.91
-8-6
-4
-2
2
4
6
幅值谱：
05101520253035
20406080100120140160180
相位谱：
05101520253035
-4
-3-2-101234
实验四
4、将信号)240sin()(t t x π=做离散序列，比较原序列与经过FFT 和IFFT 变换后的序列，并做出说明。
解题分析：
先将信号进行离散化抽样，然后求FFT 和IFFT 变换，比较原序列和变换后得到的序列。
程序设计如下： >> clear
>> t=0:1/255:1; >> x=sin(240*pi*t);
>> y=ifft(fft(x));%先进行FFT 变换，后进行IFFT 变换 >> subplot(2,1,1);
>> plot(t,x);%绘制原波形 >> subplot(2,1,2);
>> plot(t,y);%绘制变换后的波形
00.10.20.30.40.50.60.70.80.91
-1
-0.500.5100.10.20.30.40.50.60.70.80.91
-1
-0.500.51
说明：由图像可知，变换后的波形与原波形一样。
实验五
5、已知带通滤波器的系统函数为2
2100)1(2)(++=
s s
s H ，激励信号
)100c o s ()c o s 1()(t t t x *+=，
求：（1）带通滤波器的频率响应； （2）输出稳态响应并绘制图形。
解题分析：
输入系统函数，求带通滤波器的频率响应；输入激励信号，求得稳态响应，将激励信号与响应信号进行对比。
程序设计如下： >> clear;
>> w=[99,100,101];%输入信号的频率值 >> X=[0.5 1 0.5];%输入信号的幅值 >> num=[2 0];
>> den=[1 2 10001];
>> H=polyval(num,j*w)./polyval(den,j*w); >> H=freqs(num,den,w);
>> subplot(2,1,1),plot(w,abs(H)),grid; >> subplot(2,1,2),plot(w,angle(H)),grid; >> t=0:0.01:2*pi;
>> x1=X(1)*cos(99*t)+X(2)*cos(100*t)+X(3)*cos(101*t);%激励信号 >> x21=abs(X(1)*H(1))*cos(99*t+angle(X(1)*H(1)));%99频率分量 >> x22=abs(X(2)*H(2))*cos(100*t+angle(X(2)*H(2)));%100频率分量 >> x23=abs(X(3)*H(3))*cos(101*t+angle(X(3)*H(3)));%101频率分量 >> x2=x21+x22+x23; %稳态响应 >> figure(2); >> plot(t,x1); >> figure(3); >> plot(t,x2);
>> figure(4);%相互对比画图 >> subplot(2,1,1),plot(t,x1); >> subplot(2,1,2),plot(t,x2);
频谱响应：
9999.299.499.699.8100100.2100.4100.6100.8101
0.70.8
0.9
1
99
99.299.499.699.8100100.2100.4100.6100.8101
-1-0.500.51
激励信号：
01234567
-2
-1.5-1-0.500.511.52
稳态响应：
01234567
-2
-1.5-1-0.500.511.52
激励信号与响应信号进行比较：
01234567
-2
-101201234567
-2
-1012
说明：通过带通滤波器，信号的幅值得到了调整。
心得体会
通过本次实验，掌握了如何已知系统函数求解系统的频率响应，学会了DFT、IDFT、FFT、IFFT 函数的使用；还学会了新的时间间隔的定义方法t=linspace(0,T,N)。除此之外，对于如何用Matlab 作出各种各样的图，有了更加全面的认识。通过绘图，可以直观的了解信号与系统的特性，对信号与系统进行分析。
现代数字信号处理仿真作业 1.仿真题3.17 仿真结果及图形： 图 1 基于FFT的自相关函数计算
图 3 周期图法和BT 法估计信号的功率谱 图 2 基于式3.1.2的自相关函数的计算
图 4 利用LD迭代对16阶AR模型的功率谱估计16阶AR模型的系数为： a1=-0.402637623107952-0.919787323662670i; a2=-0.013530139693503+0.024214641171318i; a3=-0.074241889634714-0.088834852915013i; a4=0.027881022353997-0.040734794506749i; a5=0.042128517350786+0.068932699075038i; a6=-0.0042799971761507 + 0.028686095385146i; a7=-0.048427890183189 - 0.019713457742372i; a8=0.0028768633718672 - 0.047990801912420i a9=0.023971346213842+ 0.046436389191530i; a10=0.026025963987732 + 0.046882756497113i; a11= -0.033929397784767 - 0.0053437929619510i; a12=0.0082735406293574 - 0.016133618316269i; a13=0.031893903622978 - 0.013709547028453i ; a14=0.0099274520678052 + 0.022233240051564i; a15=-0.0064643069578642 + 0.014130696335881i; a16=-0.061704614407581- 0.077423818476583i. 仿真程序（3_17）: clear all clc %% 产生噪声序列 N=32; %基于FFT的样本长度
现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 (一)Ansoft公司的仿真工具 现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 Ansoft的信号完整性工具采用一个仿真可解决全部设计问题： SIwave是一种创新的工具，它尤其适于解决现在高速PCB和复杂IC封装中普遍存在的电源输送和信号完整性问题。 该工具采用基于混合、全波及有限元技术的新颖方法，它允许工程师们特性化同步开关噪声、电源散射和地散射、谐振、反射以及引线条和电源/地平面之间的耦合。该工具采用一个仿真方案解决整个设计问题，缩短了设计时间。 它可分析复杂的线路设计，该设计由多重、任意形状的电源和接地层，以及任何数量的过孔和信号引线条构成。仿真结果采用先进的3D图形方式显示，它还可产生等效电路模型，使商业用户能够长期采用全波技术，而不必一定使用专有仿真器。 (二)SPECCTRAQuest Cadence的工具采用Sun的电源层分析模块： Cadence Design Systems的SpecctraQuest PCB信号完整性套件中的电源完整性模块据称能让工程师在高速PCB设计中更好地控制电源层分析和共模EMI。 该产品是由一份与Sun Microsystems公司签署的开发协议而来的，Sun最初研制该项技术是为了解决母板上的电源问题。 有了这种新模块，用户就可根据系统要求来算出电源层的目标阻抗；然后基于板上的器件考虑去耦合要求，Shah表示，向导程序能帮助用户确定其设计所要求的去耦合电容的数目和类型；选择一组去耦合电容并放置在板上之后，用户就可运行一个仿真程序，通过分析结果来发现问题所在。 SPECCTRAQuest是CADENCE公司提供的高速系统板级设计工具，通过它可以控制与PCB layout相应的限制条件。在SPECCTRAQuest菜单下集成了一下工具： （1）SigXplorer可以进行走线拓扑结构的编辑。可在工具中定义和控制延时、特性阻抗、驱动和负载的类型和数量、拓扑结构以及终端负载的类型等等。可在PCB详细设计前使用此工具，对互连线的不同情况进行仿真，把仿真结果存为拓扑结构模板，在后期详细设计中应用这些模板进行设计。 （2）DF/Signoise工具是信号仿真分析工具，可提供复杂的信号延时和信号畸变分析、IBIS 模型库的设置开发功能。SigNoise是SPECCTRAQUEST SI Expert和SQ Signal Explorer Expert进行分析仿真的仿真引擎，利用SigNoise可以进行反射、串扰、SSN、EMI、源同步及系统级的仿真。 （3）DF/EMC工具——EMC分析控制工具。 （4）DF/Thermax——热分析控制工具。 SPECCTRAQuest中的理想高速PCB设计流程： 由上所示，通过模型的验证、预布局布线的space分析、通过floorplan制定拓朴规则、由规
HFSS天线仿真实验报告 半波偶极子天线设计 通信0905 杨巨 U200913892 2012-3-7
半波偶极子天线仿真实验报告 一、实验目的 1、学会简单搭建天线仿真环境的方法，主要是熟悉HFSS软件的使用方法 2、了解利用HFSS仿真软件设计和仿真天线的原理、过程和方法 3、通过天线的仿真，了解天线的主要性能参数，如驻波比特性、smith圆图特性、方向图 特性等 4、通过对半波偶极子天线的仿真，学会对其他类型天线仿真的方法 二、实验仪器 1、装有windows系统的PC一台 2、HFSS13.0软件 3、截图软件 三、实验原理 1、首先明白一点：半波偶极子天线就是对称阵子天线。 2、 对称振子是中间馈电，其两臂由两段等长导线构成的振子天线。一臂的导线半径为a，长度为l。两臂之间的间隙很小，理论上可以忽略不计，所以振子的总长度L=2l。对称振子的长度与波长相比拟，本身已可以构成实用天线。 3、 在计算天线的辐射场时，经过实践证实天线上的电流可以近似认为是按正弦律分布。取图1的坐标，并忽略振子损耗，则其电流分布可以表示为： 式中，Im为天线上波腹点的电流；k=w/c为相移常数、根据正弦分布的特点，对称振子的末端为电流的波节点；电流分布关于振子的中心店对称；超过半波长就会出现反相电流。 4、 在分析计算对称振子的辐射场时，可以把对称振子看成是由无数个电流I(z)、长度为dz的电流元件串联而成。利用线性媒介中电磁场的叠加原理，对称振子的辐射场是这些电流元辐射场之矢量和。
电流元I(z)dz所产生的辐射场为 图2 对称振子辐射场的计算 如图2 所示，电流元I(z)所产生的辐射场为 其中 5、方向函数 四、实验步骤 1、设计变量 设置求解类型为Driven Model 类型，并设置长度单位为毫米。 提前定义对称阵子天线的基本参数并初始化 2、创建偶极子天线模型，即圆柱形的天线模型。 其中偶极子天线的另外一个臂是通过坐标轴复制来实现的。 3、设置端口激励 半波偶极子天线由中心位置馈电，在偶极子天线中心位置创建一个平行于YZ面的矩形面作为激励端口平面。 4、设置辐射边界条件 要在HFSS中计算分析天线的辐射场，则必须设置辐射边界条件。这里创建一个沿Z轴放置的圆柱模型，材质为空气。把圆柱体的表面设置为辐射边界条件。 5、外加激励求解设置 分析的半波偶极子天线的中心频率在3G Hz，同时添加2.5 G Hz ~3.5 G Hz频段内的扫频设置，扫频类型为快速扫频。
信号与系统仿真实验报告1．实验目的 了解MATLAB的基本使用方法和编程技术，以及Simulink平台的建模与动态仿真方法，进一步加深对课程内容的理解。 2．实验项目 信号的分解与合成，观察Gibbs现象。 信号与系统的时域分析，即卷积分、卷积和的运算与仿真。 信号的频谱分析，观察信号的频谱波形。 系统函数的形式转换。 用Simulink平台对系统进行建模和动态仿真。 3．实验内容及结果 3.1以周期为T，脉冲宽度为2T1的周期性矩形脉冲为例研究Gibbs现象。 已知周期方波信号的相关参数为：x(t)=∑ak*exp(jkω),ω=2*π/T,a0=2*T1/T,ak=sin(kωT1)/kπ。画出x(t)的波形图（分别取m=1,3,7,19,79,T=4T1）,观察Gibbs现象。 m=1; T1=4; T=4*T1;k=-m:m; w0=2*pi/T; a0=2*T1/T; ak=sin(k*w0*T1)./(k*pi); ak(m+1)=a0; t=0:0.1:40; x=ak*exp(j*k'*w0*t); plot(t,real(x)); 3.2求卷积并画图 (1)已知：x1(t)=u(t-1)-u(t-2), x2(t)=u(t-2)-u(t-3)求：y(t)=x1(t)*x2(t)并画出其波形。 t1=1:0.01:2; f1=ones(size(t1)); f1(1)=0; f1(101)=0; t2=2:0.01:3; f2=ones(size(t2)); f2(1)=0; f2(101)=0; c=conv(f1,f2)/100;
t3=3:0.01:5; subplot(311); plot(t1,f1);axis([0 6 0 2]); subplot(312); plot(t2,f2);axis([0 6 0 2]); subplot(313); plot(t3,c);axis([0 6 0 2]); (2）已知某离散系统的输入和冲击响应分别为：x[n]=[1,4,3,5,1,2,3,5], h[n]=[4,2,4,0,4,2].求系 统的零状态响应，并绘制系统的响应图。 x=[1 4 3 5 1 2 3 5]; nx=-4:3; h=[4 2 4 0 4 2]; nh=-3:2; y=conv(x,h); ny1=nx(1)+nh(1); ny2=nx(length(nx))+nh(length(nh)); ny=[ny1:ny2]; subplot(311); stem(nx,x); axis([-5 4 0 6]); ylabel('输入') subplot(312); stem(nh,h); axis([-4 3 0 5]); ylabel('冲击效应') subplot(313); stem(ny,y); axis([-9 7 0 70]); ylabel('输出'); xlabel('n'); 3.3 求频谱并画图 (1) 门函数脉冲信号x1(t)=u(t+0.5)-u(t-0.5) N=128;T=1; t=linspace(-T,T,N); x=(t>=-0.5)-(t>=0.5); dt=t(2)-t(1); f=1/dt; X=fft(x); F=X(1:N/2+1); f=f*(0:N/2)/N; plot(f,F)
实验6 数字滤波器的网络结构 一、实验目的： 1、加深对数字滤波器分类与结构的了解。 2、明确数字滤波器的基本结构及其相互间的转换方法。 3、掌握用MA TLAB 语言进行数字滤波器结构间相互转换的子函数及程序编写方法。 二、实验原理： 1、数字滤波器的分类 离散LSI 系统对信号的响应过程实际上就是对信号进行滤波的过程。因此，离散LSI 系统又称为数字滤波器。 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器（FIR ）和无限长单位脉冲响应滤波器（IIR ）。 一个离散LSI 系统可以用系统函数来表示： M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 也可以用差分方程来表示： N M k m k=1 m=0 y(n)+a y(n-k)=b x(n-m)∑∑ 以上两个公式中，当a k 至少有一个不为0时，则在有限Z 平面上存在极点，表达的是以一个IIR 数字滤波器；当a k 全都为0时，系统不存在极点，表达的是一个FIR 数字滤波器。FIR 数字滤波器可以看成是IIR 数字滤波器的a k 全都为0时的一个特例。 IIR 数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、直接Ⅲ型、级联型和并联型。 FIR 数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型等。本实验对线性相位型及频率采样型不做讨论，见实验10、12。 另外，滤波器的一种新型结构——格型结构也逐步投入应用，有全零点FIR 系统格型结构、全极点IIR 系统格型结构以及全零极点IIR 系统格型结构。 2、IIR 数字滤波器的基本结构与实现 （1）直接型与级联型、并联型的转换 例6-1 已知一个系统的传递函数为 -1-2-3 -1-2-3 8-4z +11z -2z H(z)=1-1.25z +0.75z -0.125z 将其从直接型（其信号流图如图6-1所示）转换为级联型和并联型。
Introduction Consider the proverb, “It takes a village to raise a child.” Similarly, multiple design team members participate in assuring PCB power integrity (PI) as a design moves from the early concept phase to becoming a mature product. On the front end, there’s the electrical design engineer who is responsible for the schematic. On the back end, the layout designer handles physical implemen-tation. Typically, a PI analysis expert is responsible for overall PCB PI and steps in early on to guide the contributions of others. How quickly a team can assure PCB PI relates to the effectiveness of that team. In this paper, we will take a look at currently popular analysis approaches to PCB PI. We will also introduce a team-based approach to PCB PI that yields advantages in resource utilization and analysis results. Common Power Integrity Analysis Methods There are two distinct facets of PCB PI – DC and AC. DC PI guarantees that adequate DC voltage is delivered to all active devices mounted on a PCB (often using IR drop analysis). This helps to assure that constraints are met for current density in planar metals and total current of vias and also that temperature constraints are met for metals and substrate materials. AC PI concerns the delivery of AC current to mounted devices to support their switching activity while meeting constraints for transient noise voltage levels within the power delivery network (PDN). The PDN noise margin (variation from nominal voltage) is a sum of both DC IR drop and AC noise. DC PI is governed by resistance of the metals and the current pulled from the PDN by each mounted device. Engineers have, for many years, applied resistive network models for approximate DC PI analysis. Now that computer speeds are faster and larger addressable memory is available, the industry is seeing much more application of layout-driven detailed numerical analysis techniques for DC PI. Approximation occurs less, accuracy is higher, and automation of How a Team-Based Approach to PCB Power Integrity Analysis Yields Better Results By Brad Brim, Sr. Staff Product Engineer, Cadence Design Systems Assuring power integrity of a PCB requires the contributions of multiple design team members. Traditionally, such an effort has involved a time-consuming process for a back-end-focused expert at the front end of a design. This paper examines a collaborative team-based approach that makes more efficient use of resources and provides more impact at critical points in the design process. Contents Introduction (1) Common Power Integrity Analysis Methods (1) Applying a Team-Based Approach to Power Integrity Analysis (3) Summary (6) For Further Information (7)
实验二电磁波发射天线的模拟仿真电动力学实验报告电磁波发射天线的模拟仿真 学院: 应用科学学院专业班级: 学生姓名: 某某某 学号: 指导教师: 完成时间: 2013年7月2号 一、实验目的 1(熟悉并了解CST 的软件环境。 2(通过实验掌握天线的实际画法及步骤。 3(了解电磁波发射天线的模拟仿真过程，进一步了解电磁波发射现象。 二、实验原理及要求 在CST微波工作室中，通常采用瞬态求解器来计算天线，典型的天线特性，如S参量(S参数)、主瓣方向、增益、效率等，都将被自动计算和显11 示。按照如下图的天线模型形自行设计可接受2GHz左右的电磁波信号的天线并仿真出结果，同时作出一定分析。(碳纳米管的半径为R，轴向方向沿z轴，长度为L，中间馈电端口缝隙为D) 三、实验步骤 1、选择天线模板 启动CST，在弹出的“Welcome”对话框中点击“OK” 按钮，创建一个新项目。然后会看到选择模板对话框，选择 Antenna(Horn,Waveguide)，并点击OK按钮。 2、设置单位
用鼠标左键单击主菜单上的
实验5 抽样定理 一、实验目的： 1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。 2、进一步加深对时域、频域抽样定理的基本原理的理解。 3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和插公式的编程方法。 二、实验原理： 1、时域抽样与信号的重建 （1）对连续信号进行采样 例5-1 已知一个连续时间信号sin sin(),1Hz 3 ππ=0001f(t)=(2f t)+6f t f ，取最高有限带宽频率f m =5f 0，分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。 程序清单如下： %分别取Fs=fm ，Fs=2fm ，Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2; f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f); axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3; fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end 程序运行结果如图5-1所示：
原连续信号和抽样信号 图5-1 （2）连续信号和抽样信号的频谱 由理论分析可知，信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。因此，我们对上述三种情况下的时域信号求幅度谱，来进一步分析和验证时域抽样定理。 例5-2编程求解例5-1中连续信号及其三种抽样频率（F s>2f m、F s=2f m、F s<2f m）下的抽样信号的幅度谱。 程序清单如下： dt=0.1;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm; t=-2:dt:2;N=length(t); f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); wm=2*pi*fm;k=0:N-1;w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2;N=length(n); f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); wm=2*pi*fs;k=0:N-1; w=k*wm/N;F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end 程序运行结果如图5-2所示。 由图可见，当满足F s≥2f m条件时，抽样信号的频谱没有混叠现象；当不满足F s≥2f m 条件时，抽样信号的频谱发生了混叠，即图5-2的第二行F s＜2f m的频谱图，，在f m=5f0的围，频谱出现了镜像对称的部分。
作者简介：曹宇（１９６９－），男，上海人，硕士，工程师． 第６卷第 ６期 ２００６年１２月泰州职业技术学院学报 ＪｏｕｒｎａｌｏｆＴａｉｚｈｏｕＰｏｌｙｔｅｃｈｎｉｃａｌＩｎｓｔｉｔｕｔｅＶｏｌ．６Ｎｏ．６ Ｄｅｃ．２００６摘要：针对高速数字电路印刷电路板的板级信号完整性，分析了ＩＢＩＳ模型在板级信号完整 性分析中的作用。利用ＡＤＳ仿真软件，采用电磁仿真建模和电路瞬态仿真测试了某个 实际电路版图，给出了实际分析结果。 关键词：信号完整性；ＩＢＩＳ；仿真；Ｓ参数 中图分类号：ＴＰ３９１．９文献标识码：Ａ文章编号：１６７１－０１４２（２００６）０６－００３０－０３ 信号完整性（ＳＩ，ＳｉｇｎａｌＩｎｔｅｇｒｉｔｙ）的概念是针对高速数字信号提出来的。以往的数字产品，其时钟或数据频率在几十兆之内时，信号的上升时间大多在几个纳秒，甚至几十纳秒以上。数字化产品设计工程师关注最多的是“数字设计”保证逻辑正确。随着数字技术的飞速发展，原先只是在集成电路芯片设计中需要考虑的问题［１］在ＰＣＢ板级设计中正在逐步显现出来，并由此提出了信号完整性的概念。 在众多的讲述信号完整性的论文和专著中［２，３］，对信号完整性的描述都是从信号传输过程中可能出现的问题（比如串扰，阻抗匹配，电磁兼容，抖动等）本身来讨论信号完整性，对信号完整性没有一个统一的定义。事实上，信号完整性是指信号在通过一定距离的传输路径后在特定接收端口相对指定发送端口信号的还原程度，这个还原程度是指在指定的收发参考端口，发送芯片输出处及接收芯片输入处的波形需满足系统设计的要求［４］。 １、板级信号完整性分析 １．１信号完整性分析内容的确定 信号完整性分析工作是一项产品开发全流程工作，从产品设计阶段开始一直延续到产品定型。ＰＣＢ板级设计同样如此。在系统设计阶段，产品还没有进入试制，需要建立相应的系统模型并得到仿真结果以验证设计思想和设计体系正确与否，这个阶段称前仿真；前仿真通过后，产品投入试制，样品出来后再进行相应的测试和仿真，这个阶段称后仿真。假如将每一块ＰＣＢ板视为一个系统，影响这个系统正常工作的信号问题涉及到所有的硬件和软件，包括芯片、封装、ＰＣＢ物理结构、电源及电源传输网络和协议。 对系统所有部分都进行仿真验证是不现实的。应根据系统设计的要求选定部分内容进行测试仿真。本文所提及的“板级信号完整性分析”仅针对芯片引脚和走线的互连状态分析。 当被传输的信号脉冲时间参量（如上升时间、传输时间等）已缩短至和互连线上电磁波传输时间处于同一个量级时，信号在互连线上呈现波动效应，应采用微波传输线或分布电路的模型来对待互连线，从而产生了时延、畸变、回波、相邻线之间的干扰噪声等所谓的“互连效应”［１］。 假设ＰＣＢ板上芯片引脚的输入输出信号都是“干净”的，那么只要考虑互连线路本身的互连效应。事实上，每个芯片引脚在封装时都有其独特的线路特性，这些特性是由其内部的晶体管特性决定的，同样的信号在不同引脚上的传输效率差异很大。因此，在分析信号传输的互连效应时必须考虑芯片内部的电路特性以提取相对准确的电路模型，并在此基础上作进一步的分析。这个模型就是在业界被广泛使用的ＩＢＩＳ模型。 １．２ＩＢＩＳ标准模型的建立 ＰＣＢ板级信号完整性的仿真及应用 曹宇，丁志刚，宗宇伟 （上海计算机软件技术开发中心，上海２０１１１２）
电磁场与电磁兼容 实验报告 学号： 姓名： 院系： 专业： 教师： 05月20日
半波对称振子天线阵最大辐射方向控制 实验工具 ?Expert MININEC Classic电磁场数值仿真软件 实验目的 根据要求的参数，利用仿真软件设计和分析自由空间或地面上的细、直线天线的电磁场数值，并完成以下要求： ?改变每幅天线馈电电流的相位控制最大增益的方向：要求的最大增益方向是：1. 00 ；2. 400；3. 800 （选择与自己学号后2位数最近的度数） ?根据运行结果指出： 1.增益方向性图； 2.最大增益； 3.最大增益方向。 实验参数 ?频率 f = 300MHz，波长λ = 1m ?四分之一波长单极子天线L=0.25λ，四个半波长对称振子排列在一条直线上，相邻两幅天线的间隔是四分之一波长 实验过程 ?建立几何模型：点—> 线，尺寸，环境，坐标等 半波对称振子放在 YOZ 平面内，相邻振子的间距是四分之一波长 0.25m。
图1 问题描述图2 –图4 几何模型 图3 图4 ?定义电特性：频率，电压，当前节点 ZENITH(DEG) 对应球坐标系中的θ， AZIMUTH (DEG) 对应球坐标系中的φ 图5 电特性—频率图6 馈电电流相位设置
图7 球坐标参数θ、ψ以及间隔设置 ?选择模式：辐射模式 ?求解项：近场 ?调试、运行 表格中出现“No detected violations ”表明设置正确 图8 选择运行平面图9 调试结果 ?显示结果 3D display 显示所设计天线的图形 天线增益方向性图中给出了最大增益值和最大增益方向、以及半功率增益带宽的计算结果。
交通仿真实验报告 篇一：交通仿真实验报告 目录 1 上机性质与目的.................................. 2 2 上机内容....................................... 2 3 交叉口几何条件、信号配时和交通流数据描述.......... 3 3.1 交叉口几何数据................................ 3 3.2 交叉口信号配时系统............................ 3 3.3 交叉口交通流数据.............................. 4 4 交叉口交通仿真.................................. 4 4.1 交通仿真步骤.................................. 4 4.2 二维输出..................................... 13 4.3 3D输出...................................... 14 5 仿真结果分析................................... 15 6 实验总结和体会 (15) 实验上机名称：信号交叉口仿真 1 上机性质与目的 本实验属于计算机仿真实验，借助仿真系统模拟平面信号交叉口场景，学生将完成从道路条件设计到信号相位配置等一系列仿真实验。 实验目的： 1. 了解平面信号交叉口在城市交通中的地位； 2. 了解平面信号交叉口的主要形式、规模等基本情况； 3. 了解交叉口信号相位配时及对交叉口通行能力的影响；
本人技术屌丝一枚，从事PCB相关工作已达8年有余，现供职于世界闻名的首屈一指的芯片设计公司，从苦逼的板厂制板实习，到初入Pcblayout，再到各种仿真的实战，再到今天的销售工作，一步一步一路兢兢业业诚诚恳恳，有一些相关领悟和大家分享。买卖不成也可交流。 1.谈起硬件工作，是原理图，pcb，码农的结合体，如果你开始了苦逼的pcblayout工作，那么将是漫长的迷茫之路，日复一日年复一年，永远搞不完的布局，拉线。眼冒金星不是梦。最多你可以懂得各种模块的不同处理方式，各种高速信号的设计，但永远只能按照别人的意见进行，毫无乐趣。 2.谈起EDA相关软件，形象的说，就普通的PROTEL/AD来说你可能只有3-6K，对于pads 可能你有5-8K，对于ALLEGRO你可能6-10K,你会哀叹做的东西一样，却同工不同酬，没办法这就是市场，我们来不得无意义的抱怨。 3.众所周知，一个PCB从业者最好的后路就是仿真工作，为什么呢？一；你可以懂得各种模块的设计原则，可以优化不准确的部分，可以改善SI/PI可以做很多，这往往是至关重要的，你可以最大化节约成本，减少器件却功效相同；二；从一个pcblayout到仿真算是水到渠成，让路走的更远； 三：现实的说薪资可以到达11-15K or more，却更轻松，更有价值，发言权，你不愿意吗？ 现在由于本人已技术转销售，现在就是生意人了哈哈，我也查询过各种仿真资料我发现很少，最多不过是Mentor Graphics 的HyperLynx ，candense的si工具，
但是他们真的太low了，精确度和完整性根本不能保证，最多是定性的能力，无法定量。真正的仿真是完整的die到die的仿真，是完整的系统的，是需要更高级的仿真软件，被收购的xxsigrity，xx ansys，hspicexx，adxx等等，这些软件才是真正的仿真。 本人提供各种软件及实战代码，例子，从基本入门到高级仿真，从电源仿真，到ddr仿真到高速串行仿真，应有尽有，，完全可以使用，想想以后的高薪，这点投入算什么呢？舍不得孩子套不住狼哦。 所有软件全兼容32位和64位系统。 切记本人还提供学习手册，你懂的，完全快速进入仿真领域。你懂的！ 希望各位好好斟酌，自己的路是哪个方向，是否想更好的发展，舍得是哲学范畴，投资看得是利润的最大化，学会投资吧，因为他值得拥有，骚年！ 注：本人也可提供培训服务，面面俱到，形象具体，包会！ 有购买和学习培训兴趣的请联系 QQ:2941392162
[键入公司名称] [键入文档标题] 通信0905 杨巨 U2 2012-3-7 半波偶极子天线仿真实验报告 一、实验目的 1、学会简单搭建天线仿真环境的方法，主要是熟悉HFSS软件的使用方法 2、了解利用HFSS仿真软件设计和仿真天线的原理、过程和方法 3、通过天线的仿真，了解天线的主要性能参数，如驻波比特性、smith圆图特性、方向图 特性等 4、通过对半波偶极子天线的仿真，学会对其他类型天线仿真的方法 二、实验仪器 1、装有windows系统的PC一台 2、HFSS13.0软件 3、截图软件 三、实验原理 1、首先明白一点：半波偶极子天线就是对称阵子天线。
对称振子是中间馈电，其两臂由两段等长导线构成的振子天线。一臂的导线半径为a，长度为l。两臂之间的间隙很小，理论上可以忽略不计，所以振子的总长度L=2l。对称振子的长度与波长相比拟，本身已可以构成实用天线。 在计算天线的辐射场时，经过实践证实天线上的电流可以近似认为是按正弦律分布。取图1的坐标，并忽略振子损耗，则其电流分布可以表示为：式中，Im为天线上波腹点的电流；k=w/c为相移常数、根据正弦分布的特点，对称振子的末端为电流的波节点；电流分布关于振子的中心店对称；超过半波长就会出现反相电流。 4、 在分析计算对称振子的辐射场时，可以把对称振子看成是由无数个电流I(z)、长度为dz的电流元件串联而成。利用线性媒介中电磁场的叠加原理，对称振子的辐射场是这些电流元辐射场之矢量和。 电流元I(z)dz所产生的辐射场为 图2 对称振子辐射场的计算 如图2 所示，电流元I(z)所产生的辐射场为 其中 5、方向函数 四、实验步骤 1、设计变量 设置求解类型为Driven Model 类型，并设置长度单位为毫米。 提前定义对称阵子天线的基本参数并初始化 2、创建偶极子天线模型，即圆柱形的天线模型。 其中偶极子天线的另外一个臂是通过坐标轴复制来实现的。 3、设置端口激励 半波偶极子天线由中心位置馈电，在偶极子天线中心位置创建一个平行于YZ面的矩形面作为激励端口平面。 4、设置辐射边界条件 要在HFSS中计算分析天线的辐射场，则必须设置辐射边界条件。这里创建一个沿Z轴放置的圆柱模型，材质为空气。把圆柱体的表面设置为辐射边界条件。 5、外加激励求解设置 分析的半波偶极子天线的中心频率在3G Hz，同时添加2.5 G Hz ~3.5 G Hz频段内的扫频设置，扫频类型为快速扫频。 6、设计检查和运行仿真计算 7、HFSS天线问题的数据后处理 具体在实验结果中阐释。 五、实验结果 1、回波损耗S11 回波损耗回波损耗是电缆链路由于阻抗不匹配所产生的反射，是一对线自身的反射，是天线设计需要关注的参数之一。 图中所示是在2.5 G Hz ~3.5 G Hz频段内的回波损耗，设计的偶极子天线中心频率约为3 G Hz，S11<-10dBd的相对带宽BW=（3.25-2.775）/3*100%=15.83%
MATLAB通信系统仿真实验报告
实验一、MATLAB的基本使用与数学运算 目的：学习MATLAB的基本操作，实现简单的数学运算程序。 内容： 1-1要求在闭区间[0，2π]上产生具有10个等间距采样点的一维数组。试用两种不同的指令实现。 运行代码：x=[0:2*pi/9:2*pi] 运行结果： 1-2用M文件建立大矩阵x x=[0.10.20.30.40.50.60.70.80.9 1.11.21.31.41.51.61.71.81.9 2.12.22.32.42.52.62.72.82.9 3.13.23.33.43.53.63.73.83.9] 代码：x=[0.10.20.30.40.50.60.70.80.9 1.11.21.31.41.51.61.71.81.9 2.12.22.32.42.52.62.72.82.9 3.13.23.33.43.53.63.73.83.9] m_mat 运行结果： 1-3已知A=[5，6;7，8]，B=[9，10;11,12]，试用MATLAB分别计算 A+B,A*B,A.*B,A^3，A.^3，A/B,A\B. 代码：A=[56;78]B=[910;1112]x1=A+B X2=A-B X3=A*B X4=A.*B X5=A^3 X6=A.^3X7=A/B X8=A\B
运行结果： 1-4任意建立矩阵A，然后找出在[10,20]区间的元素位置。 程序代码及运行结果： 代码：A=[1252221417;111024030;552315865]c=A>=10&A<=20运行结果： 1-5总结：实验过程中，因为对软件太过生疏遇到了些许困难，不过最后通过查书与同学交流都解决了。例如第二题中，将文件保存在了D盘，而导致频频出错，最后发现必须保存在MATLAB文件之下才可以。第四题中，逻辑语言运用到了ij，也出现问题，虽然自己纠正了问题，却也不明白错在哪了，在老师的讲解下知道位置定位上不能用ij而应该用具体的整数。总之第一节实验收获颇多。
仿真作业 姓名：李亮 学号：S130101083
4.17程序 clc; clear; for i=1:500 sigma_v1=0.27; b(1)=-0.8458; b(2)=0.9458; a(1)=-(b(1)+b(2)); a(2)=b(1)*b(2); datlen=500; rand('state',sum(100*clock)); s=sqrt(sigma_v1)*randn(datlen,1); x=filter(1,[1,a],s); %% sigma_v2=0.1; u=x+sqrt(sigma_v2)*randn(datlen,1); d=filter(1,[1,-b(1)],s); %% w0=[1;0]; w=w0; M=length(w0); N=length(u); mu=0.005; for n=M:N ui=u(n:-1:n-M+1); y(n)=w'*ui; e(n)=d(n)-y(n); w=w+mu.*conj(e(n)).*ui; w1(n)=w(1); w2(n)=w(2); ee(:,i)=mean(e.^2,2); end end ep=mean(ee'); plot（ep）; xlabel('迭代次数');ylabel('MSE');title('学习曲线'); plot(w1); hold; plot(w2); 仿真结果：
步长0.015仿真结果 0.10.20.30.4 0.50.60.7迭代次数 M S E 学习曲线
步长0.025仿真结果
步长0.005仿真结果 4.18 程序 data_len = 512; %样本序列的长度 trials = 100; %随机试验的次数 A=zeros(data_len,2);EA=zeros(data_len,1); B=zeros(data_len,2);EB=zeros(data_len,1); for m = 1: trials a1 = -0.975; a2 = 0.95; sigma_v_2 =0.0731; v = sqrt(sigma_v_2) * randn(data_len, 1, trials);%产生v(n) u0 = [0 0]; num = 1; den = [1 a1 a2]; Zi = filtic(num, den, u0); %滤波器的初始条件 u = filter(num, den, v, Zi); %产生样本序列u(n) %（2）用LMS滤波器来估计w1和w2 mu1 = 0.05; mu2 = 0.005; w1 = zeros(2, data_len);