MATLAB的基本概念.pptVIP

人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * 2.2.3 仿真的一般过程和步骤 系统仿真一般分为3个步骤，及仿真建模、仿真实验和仿真分析。这3个步骤往往需要循环执行多次之后才能够获得令设计者满意的仿真结果。 （1）仿真模型、（2）仿真实验、（3）仿真分析 上述三个步骤在实际分析时应该是一个循环操作，在循环过程中，得到的结果并不十分满意，设计者还需要不断修改仿真模型，不断对取得的实验数据进行分析，直到达到预期目的为止。 2.3 MATLAB在信号与系统中的应用 2.3.1 连续信号系统的分析 这里讨论用MATLAB表示和分析连续信号系统的问题，严格地说，分析连续信号系统应该用符号推理方法，而非数值方法。因为数值只能给出采样数据，只有采样数据到达一定密度才能近似模拟连续信号系统。 【例题2-4】假定相对于采样点密度而言，信号变化足够慢的情况下。用MATLAB描述连续信号系统。 ① 单位冲激响应 ② 单位阶跃响应 ③ 复指数函数 解： 建立仿真模型 单位冲激响应 单位冲激响应函数 可以看作是宽度为 ，数学模型中用 表示，幅度为 的矩形脉冲，即 表示在 处的脉冲。 单位阶跃响应 在 处的阶跃可以写成 。即 ③ 复指数函数 若 ，它是实指数函数，如 ，则为虚指数函数，其余部分为余弦函数，虚部为正弦函数。本例题中令 MATLAB程序 clear;clc; t0=0;tf=5;dt=0.005;t1=1; t=[t0:dt:t1]; st=length(t); n=floor((t1-t0)/dt); %t1的采样信息 %单位冲激响应信号 %在t1处有一个持续时间为dt,面积为1的脉冲信号，其余时间为0 x1=zeros(1,st); x1(n)=1/dt; %给出t1处的脉冲信号 subplot(2,2,1),stairs(t,x1) grid on %绘制图形,图为第一行左1 axis([0,5,0,25]) %设置图形坐标 %单位阶跃响应信号 %从t0到tf,在t1前信号为0，从t1处信号跳变为1 x2=[zeros(1,n-1),ones(1,st-n+1)]; %产生阶跃信号 subplot(2,2,2),stairs(t,x2),grid on %绘制图形，图为第一行左2 axis([0,5,0,1.5]) %设置图形坐标 %复指数信号 alpha=-0.5;w=10; x3=exp((alpha+j*w)*t); %产生复指数信号 subplot(2,2,3),plot(t,real(x3)),grid on %绘制图形，图为第二行左1 subplot(2,2,4),plot(t,imag(x3)),grid on %绘制图形，图为第二行左2 图2-3 例题2-4运行结果 从上图中可以看出，复数指数信号可以分解为余弦信号和正弦信号，它们分别是复数信号的实部和虚部。图中第二行的两个衰减振荡信号就代表 这两个分量的相位差是 2.3.2 MATLAB在模拟通信链路中的应用 【例题2-5】卷积编码器和解码器在模拟通信链路中的应用。 解： MATLAB程序及仿真图形 EbNo = 4.5:.5:7; linEbNo = 10.^(EbNo(:).*0.1); M = 4; codeRate = 1/2; constlen = 7; k = log2(M); codegen = [171 133]; tblen = 32; % 反馈长度 trellis = poly2trellis(constlen, codegen); dspec = distspec(trellis, 7); expVitBER = bercoding(EbNo, conv, hard, codeRate, dspec); semilogy(EbNo, expVitBER, g); xlabel(Eb/N