搜档网
当前位置:搜档网 › 信号与系统Matlab课程设计报告

信号与系统Matlab课程设计报告

信号与系统Matlab课程设计报告
信号与系统Matlab课程设计报告

南通大学电子信息学院信号与系统课程设计报告

班级:信131

姓名:罗开元

学号:1311011022

指导老师:张小美

设计时间:2015.01.13

2014—2015学年第一学期

一、连续信号的时域分析

二、

1. 信号的产生 (1)阶跃函数

function [t,y1]=jieyue(t1,t2,t0) dt=0.01;

ttt=t1:dt:t0-dt; tt=t0:dt:t2; t=t1:dt:t2; n=length(ttt); nn=length(tt); u=zeros(1,n); uu=ones(1,nn); y1=[u,uu]; return 冲激函数

function [t,y2]=chongji(t1,t2,t0) dt=0.01; t=t1:dt:t2; n=length(t); y2(1:n)=0;

y2(1,(t0-t1)/dt+1)=1/dt;

(2)调用上述函数产生信号)2-t ε(,)(4-t δ,-t e )(t ε,-6s ≤t ≤6s,并画出波形。

Command Window

subplot(3,1,1); [t1,y1]=jieyue(-6,6,2);

stairs(t1,y1);

axis([-6 6 0 1.5]);

subplot(3,1,2);

[t2,y2]=chongji(-6,6,4);plot(t2,y2); subplot(3,1,3);

[t3,y3]=jieyue(-6,6,0);

y3=exp(-(t3)).*y3;plot(t3,y3);

波形如下图所示:

(3)根据f(t)画出f(2t)和f(1-0.5t)的波形

t=-3:0.01:3;

y=tripuls(t,4,0.6);

subplot(3,1,1);

plot(t,y);

title('f(t)');

xlabel('(a)');

y1=tripuls(2*t,4,0.6); subplot(3,1,2);

plot(t,y1);

title('f(2t)');

xlabel('(b)');

t1=2-2*t;

y2=tripuls(1-0.5*t1,4,0.6); subplot(3,1,3);

plot(t1,y2);

title('f(1-0.5*t)');

xlabel('(c)');

得到波形如下图所示:

已知信号f(t)=(1+t/2)*(u(t+2)-u(t-2)),用matlab求f(t+2),f(t-2),f(-t).f(2t),-f(t),并绘出时域波形。

Command Window

syms t;

f=sym(‘(t/2+1)*(Heaviside(t+2)- Heaviside(t-2))’);

subplot(2,3,1);ezplot(f,[-3,3]);

y1=subs(f,t,t+2);

subplot(2,3,2);ezplot(y1,[-5,1]);

y2=subs(f,t,t-2);

subplot(2,3,3);ezplot(y2,[-1,5]);

y3=subs(f,t,-t);

subplot(2,3,4);ezplot(y3,[-3,3]);

y4=subs(f,t,2*t);

subplot(2,3,5);ezplot(y4,[-2,2]);

y5=-f;

subplot(2,3,6);ezplot(y5,[-3,3]);

仿真结果如下图所示:

2、信号的卷积

(1)f (t )=f1(t)*f2(t)

function [k,f]=myconv(f1,f2,k1,k2,p) f=conv(f1,f2)*p; k0=k1(1)+k2(1);

k3=length(f1)+length(f2)-2; k=k0:p:k0+k3*p;

(2)求x1=[])()(2-t -t e t εε-和x2=δ(t+3)+δ(t-3)的卷积x1(t)*x2(t),并验证卷积的性质。

Command Window

[t1,f11]=jieyue(-8,8,0);[t1,f12]=jieyue(-8,8,2); f1=f11-f12; x1=exp(-t1).*f1;

[t2,f21]=chongji(-8,8,-3);[t2,f22]=chongji(-8,8,3);

x2=f21+f22;

subplot(3,1,1);plot(t1,x1); subplot(3,1,2);plot(t2,x2); [t3,f]=myconv(x1,x2,t1,t2,0.01) subplot(3,1,3);plot(t3,f); 仿真结果:

(3)已知???≤≤-<<=32,12

0,1t f 1

t t )(,t 22e t f -=)(,0<=t<=10的卷积f(t)=f1(t)*f2(t)的时域波形图。

Command Window

t11=0; t12=3; t21=0; t22=10;

t1=t11:0.001:t12;

ft1=-sign(t1-2);

t2=t21:0.001:t22;

ft2=exp(-2*t2);

t=t11+t21:0.001:t12+t22; ft=conv(ft1,ft2);

ft=ft*0.001;

subplot(2,2,1);

plot(t1,ft1);

title('f1(t)');

subplot(2,2,2);

plot(t2,ft2);

title('f2(t)');

subplot(2,2,3);

plot(t,ft);

title('f1(t)*f2(t)');

仿真结果:

例:已知两个信号f1(t)=u(t-1)-u(t-2),f2(t)=u(t)-u(t-1),求f(t)=f1(t)*f2(t)的时域波形图。

Command Window

t1=1;

t2=2;

t3=0;

t4=1;

t=0:T:t2+t4;

x1=ones(size(t)).*((t>t1)-(t>t2));

x2=ones(size(t)).*((t>t3)-(t>t4));

y=conv(x1,x2)*T;

subplot(3,1,1),plot(t,x1);

ylabel('x1(t)');

subplot(3,1,2),plot(t,x2);

ylabel('x2(t)');

subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1));

ylabel('y(t)=x1*x2');

xlable('---t/s');

仿真结果:

3.连续系统的响应

(1)已知系统的微分方程为y''(t)+5y'(t)+6y(t)=f(t),求系统的单位冲激响应h(t)和单位阶跃响应g(t);

CommandWindow

b=[1];a=[1,5,6];

subplot(2,1,1);impulse(b,a);

subplot(2,1,2);step(b,a);

仿真结果:

(2)对于上述系统,请画出激励f(t)分别为)(t e -t ε、)(t e -2t ε、)()(t t 2cos ε、

)(t t 2ε时系统的零状态响应的波形,分析与理论计算的结果是否相符。

Command Window

b=[1]; a=[1,5,6]; t=0:0.1:10; f1=exp(-t); f2=cos(2*t); f3=t.^2; f4=(exp(-2*t));

subplot(2,2,1);lsim(b,a,f1,t);grid;title('f(t)=exp(-t)'); subplot(2,2,2);lsim(b,a,f2,t);grid;title('f(t)=cos(2*t)'); subplot(2,2,3);lsim(b,a,f3,t);grid;title('f(t)=t.^2'); subplot(2,2,4);lsim(b,a,f4,t);grid;title('f(t)=exp(-2*t)');

仿真结果:

(3)已知系统的微分方程为y ’’(t)+3y ’(t)+2y(t)=3f(t)+f ’(t),初始条件:y(0+)=1,y ’(0-)=2,求:

1)系统的零输入响应)

(t y x ; 2)激励为f(t)=)(t e -3t 时,系统的零状态响应)(t y f 和全响应y(t),分析与理论计算的结果是否相符。

二.离散系统的时域分析

(1)已知离散系统的差分方程为:

y(k)+1/3y(k-2)=1/6f(k)+1/2f(k-1)+1/2f(k-2)+1/6f(k-3),

画出单位序列响应、单位阶跃响应的波形。

Command Window

k=0:1:32;

a=[1,0,1/3,0];

b=[1/6,1/2,1/2,1/6];

subplot(2,1,1);hk=impz(b,a,k);stem(k,hk,'k');

subplot(2,1,2);gk=dstep(b,a,k);stem(k,gk,'r');

仿真结果:

(2)已知离散系统的差分方程为:

y(k)-1.5y(k-1)+0.5y(k-2)=f(k)满足初始条件y(-1)=4,y(-2)=10,

用filtic和filter子函数求系统的激励为f(k)=(0.25) ^k*u(k)时的零输入、零状态以及完全响应。

Command Window

a=[1,-1.5,0.5];

b=[1];

N=20;

k=0:N-1;

f=0.25.^k;

f0=zeros(1,N);

y01=[4,10];

fi=filtic(b,a,y01);

y0=filter(b,a,f0,fi);

fi0=filtic(b,a,0);

y1=filter(b,a,f,fi0);

y=filter(b,a,f,fi);

y2=((1/2).^k+1/3*(1/4).^k+2/3).*ones(1,N);

subplot(2,3,1);stem(k,f);title('输入信号f(k)');

subplot(2,3,2);stem(k,y0);title('零输入响应');

subplot(2,3,3);stem(k,y1);title('零状态响应');

subplot(2,3,4);stem(k,y);title('用filter求完全响应');

subplot(2,3,5);stem(k,y2);title('用公式求完全响应');

仿真结果:

(3)已知离散系统的差分方程为y (k )+3y(k-1)+2y(k-2)=f(k),

(a )画出单位阶跃响应、单位序列响应的波形;

(b )画出激励)()(k 2k f k ε=时的系统零状态响应波形。

Command Window

n=0:1:100; a=[1,3,2]; b=[1,0,0];

subplot(3,1,1);hn=impz(b,a,n);stem(n,hn,'k'); subplot(3,1,2);gn=dstep(b,a,n);stem(n,gn,'r'); f2=2.^n;

subplot(3,1,3);y2=filter(b,a,f2);plot(n,y2,'.'); 仿真结果:

三.信号与系统的频域分析

1、门函数的频谱

,画出τ=10秒时门函数在-2(rad/s)<=w<=2(1)产生宽度为τ的门函数)

g

(t

τ

π(rad/s)频率范围内频谱,记录最大值,观察第一过零点位置:Command Window

dt=0.1;

N=500;

door_width=10;

tao=door_width/2;

t1=-(N-dt):dt:-tao+dt;

t2=-tao:dt:tao;

t3=tao+dt:dt:N-dt;

t=[t1,t2,t3];

f=[zeros(1,length(t1)),ones(1,length(t2)),zeros(1,length(t3))];

w=-2*pi:0.1:2*pi;

F1=f*exp(-j*t'*w)*dt;

plot(w,real(F1));grid;

仿真结果:

(2)改变τ的值为5秒和20秒,重复(1)步骤;

5=τ 20=τ

2、傅里叶变换及其性质的验证 (1)求f(t)=t

2e

-的傅里叶变换,并绘出f(t)及其傅里叶变换的波形图。

Command Window

syms t ;

ft=exp(-2*abs(t)); Fw=fourier(ft);

subplot(2,1,1);ezplot(ft); subplot(2,1,2);ezplot(Fw); 仿真结果:

(2)以门函数)(t g 为分析对象,验证傅里叶变换的时移性质、频移性质; 时移性质:

Command Window

dt=0.1; N=500; door__width=10; tao=door__width/2; t1=-(N-1):dt:-tao+1; t2=-tao:dt:tao; t3=tao-1:dt:N-1;

妇产科学病例分析题

《妇产科学》病例分析题 《妇产科学》病例分析题 1、孕妇35岁,G1P0,孕36周,浮肿(++),血压21.3/14.6kPa (160/110mmH蛋白尿(++),一天前出现头痛,恶心。入院诊断为重度妊高征,先兆子痫,(1)试写出此孕妇3个主要的护理诊断。(2)该孕妇目前应接受哪些护理措施。(3)目前用硫酸镁治疗,写出硫酸镁使用的注意事项。答案:(1)知识缺乏:缺乏对妊高征处理的相关知识有母儿受伤的危险:与子痫发作时病人意识丧失有关焦虑:与担心高血压及其对母儿的影响有关(2) 1、卧床休息,以左侧卧位为宜 2、指导病人摄入足够的水和富含纤维素的食品,及足够蛋白质的摄入 3、硫酸镁的用药护理 4、重度妊高征患者,保持病情稳定、预防子痫发生,为分娩做好准备①将病人安排在安静、光线较暗的病室,医护活动尽量集中;②准备下列物品:呼叫器、放好床档防坠床、急救药物用物、产包。(3)硫酸镁的治疗浓度和中毒浓度相近,故在使用过程中严密观察其毒性作用,控制入量,滴速以1g/h,不超过2g/h。同时监测以下指标: 1、膝腱反射必须存在; 2、呼吸不少于16次/分钟; 3、尿量不少于25ml/h;尿少提示肾脏排泄功能受到抑制,镁离子蓄积,故应随时备好10%葡萄糖酸钙,及时给予解

毒。 2、一孕妇,26岁,孕27W(夜晚无异常否认性生活)晨起时发一内裤及床上有鲜血,站立时也随着出血感,立即来院,入院后给予保守治疗,问病史曾有2次怀孕,自然流产,强烈要求保住此胎,胎儿安好,给母亲对症,支持期待疗法,阴道间断、少量出血,孕32W晚2Am,突然阴道大量出血,蹲痰盂,500mL,仍有活动出血;请作出诊断(相应的检查)及处理?(12分)诊断:前置胎盘答案:进一步检查:(1)观察生面体征。(2)B超观察胎儿、胎盘、羊水。(3)胎监监护胎儿处理:(1)促进胎儿肺成熟;(2)绝对卧床、左侧卧位、吸氧;(3)静脉给予补液,立即准备剖宫产手术,结束妊娠。 3、某患者,LMP:2004年3月10日,停经50天时,阴道流血,量少,有早孕反应,于2004年5月3日入院。无痛或有轻微下腹疼痛,伴腰痛及下坠感。妇科检查子宫颈口闭,子宫大小与停经月份符合,B超示:宫内可见一妊娠囊,有胚芽和胎心音波动。请给该孕妇作出诊断、治疗。(8分)答案:诊断:先兆流产治疗:(1)卧床休息,吸氧;(2)营养饮食;(3)禁止性生活(4)黄体酮20mg im QD;VitE20mg tid(5)观察胎儿宫内情况。 4:女性,15岁,13岁初潮,近3个月月经紊乱,经量较多,肛诊:子宫正常大小,双侧附件区软,未触及包块。(1)该患最可能的诊断是 什么?(2)拟给该患调整月经周期,宜

信号与系统课程设计报告材料

课程设计报告 课程名称信号与系统课程设计指导教师 设计起止日期 学院信息与通信工程 专业电子信息工程 学生 班级/学号 成绩 指导老师签字

目录 1、课程设计目的 (1) 2、课程设计要求 (1) 3、课程设计任务 (1) 4、课程设计容 (1) 5、总结 (11) 参考文献 (12) 附录 (12)

1、课程设计目的 “信号与系统”是一门重要的专业基础课,MATLAB作为信号处理强有力的计算和分析工具是电子信息工程技术人员常用的重要工具之一。本课程设计基于MATLAB完成信号与系统综合设计实验,以提高学生的综合应用知识能力为目标,是“信号与系统”课程在实践教学环节上的必要补充。通过课设综合设计实验,激发学生理论课程学习兴趣,提高分析问题和解决问题的能力。 2、课程设计要求 (1)运用MATLAB编程得到简单信号、简单信号运算、复杂信号的频域响应图; (2)通过对线性时不变系统的输入、输出信号的时域和频域的分析,了解线性时不变系统的特性,同时加深对信号频谱的理解。 3、课程设计任务 (1)根据设计题目的要求,熟悉相关容的理论基础,理清程序设计的措施和步骤; (2)根据设计题目的要求,提出各目标的实施思路、方法和步骤; (3)根据相关步骤完成MATLAB程序设计,所编程序应能完整实现设计题目的要求; (4)调试程序,分析相关理论; (5)编写设计报告。 4、课程设计容 (一)基本部分 (1)信号的时频分析 任意给定单频周期信号的振幅、频率和初相,要求准确计算出其幅度谱,并准确画出时域和频域波形,正确显示时间和频率。 设计思路: 首先给出横坐标,即时间,根据设定的信号的振幅、频率和初相,写出时域波形的表达式;然后对时域波形信号进行傅里叶变化,得到频域波形;最后使用plot函数绘制各个响应图。 源程序: clc; clear; close all; Fs =128; % 采样频率 T = 1/Fs; % 采样周期 N = 600; % 采样点数 t = (0:N-1)*T; % 时间,单位:S x=2*cos(5*2*pi*t);

MATLAB课程设计报告

华东交通大学MATLAB程序设计报告书 课题名称:基于MATLAB的粒子群优化算法的实现 姓名: 学号:20160280800014 专业:控制科学与工程 2016年 11月 20日

基于MATLAB的粒子群优化算法的实现 一、课程选题目的 本次课程设计的课题为《基于MATLAB的粒子群优化算法的实现》,主要为学会运用MATLAB对实际算法编程,加深对粒子群优化算法的理解,并为今后熟练使用MA TLAB进行系统的分析仿真和设计奠定基础。数值计算分析可以帮助更深入地理解理论知识,并为将来使用MA TLAB进行各领域数值分析分析和实际应用打下基础。 此次课程主要是为了进一步熟悉对MATLAB软件的使用,以及学会利用MA TLAB对数值运算这种实际问题进行处理,将理论应用于实际,加深对它的理解。 二、粒子群优化算法原理 优化是科学研究、工程技术和经济管理等领域的重要研究工具。它所研究的问题是讨论在众多的方案中寻找最优方案。例如,工程设计中怎样选择设计参数,使设计方案既满足设计要求又能降低成本;资源分配中,怎样分配有限资源,使分配方案既能满足各方面的基本要求,又能获得好的经济效益。在人类活动的各个领域中,诸如此类,不胜枚举。优化这一技术,正是为这些问题的解决,提供理论基础和求解方法,它是一门应用广泛、实用性很强的科学。近十余年来,粒子群优化算法作为群体智能算法的一个重要分支得到了广泛深入的研究,在路径规划等许多领域都有应用。 2.1 粒子群优化算法的起源 粒子群优化(PSO)算法是由Kennedy和Eberhart于1995年用计算机模拟鸟群觅食这一简单的社会行为时,受到启发,简化之后而提出的。 设想这样一个场景:一群鸟随机的分布在一个区域中,在这个区域里只有一块食物。所有的鸟都不知道食物在哪里。但是他们知道当前的位置离食物还有多远。那么找到食物的最优策略是什么呢。最简单有效的方法就是追寻自己视野中目前离食物最近的鸟。如果把食物当作最优点,而把鸟离食物的距离当作函数的适应度,那么鸟寻觅食物的过程就可以当作一个函数寻优的过程。鱼群和鸟群的社会行为一直引起科学家的兴趣。他们以特殊的方式移动、同步,不会相互碰撞,整体行为看上去非常优美。生物学家CargiReynolds提出了一个非常有影响的鸟群聚集模型。在他的模拟模型boids中,每一个个体遵循:避免与邻域个体相冲撞、匹配邻域个体的速度、试图飞向感知到的鸟群中心这三条规则形成简单的非集中控制算法驱动鸟群的聚集,在一系列模拟实验中突现出了非常接近现实鸟群聚集行为的现象。该结果显示了在空中回旋的鸟组成轮廓清晰的群体,以及遇到障碍物时鸟群的分裂和再度汇合过程。由此受到启发,经过简化提出了粒子群优化算法。 2.2粒子群优化算法的原理 在粒子群优化算法中,每个优化问题的潜在解都是搜索空间中的一只鸟,称之为“粒子”。所有的粒子都有一个由被优化的函数决定的适应值,每个粒子还有一个速度决定他们飞翔的方向和距离。然后粒子们就追随当前的最优粒子在解空间中搜索。优化开始时先初始化为一群随机粒子(随机解)。然后通过迭代找到最优解。在每一次迭代中,粒子通过跟踪两个极值来更新自己。第一个极值就是整个种群目前找到的最优解。这个极值是全局极值。另外也可以不用整个种群而只是用其中一部分作为粒子的邻居,那么在所有邻居中的极值就是局部极值。第二个极值是粒子本身所找到的最优解,称为个体极值。这是因为粒子仅仅通过跟踪全局极值或者局部极值来更新位置,不可能总是获得较好的解。这样在优化过程中,粒子在追随全局极值或局部极值的同时追随个体极值则圆满的解决了这个问题。这就是粒子群优化

2015年妇产科临床护理案例分析题

2015年妇产科临床护理案例分析题 天宇考王卫生资格高级职称考试题库包含:章节练习、综合复习题、模拟试卷、考前冲刺、历年真题等。 ------------------------第1题------------------------ 题型: 【共用题干单选题】 大问: 男性,30岁。夜间发作性腹部烧灼样痛二月余,进食后能迅速缓解,昨起排柏油样便3次,今晨起床晕倒。体检:T 37℃,P 120次/min,R 24次/min,BP 80/50mmHg。神志恍惚,皮肤苍白,四肢厥冷。 小问: 目前的情况应考虑发生 答案选项: A.穿孔复选 B.幽门梗阻复选 C.肠系膜动脉栓塞复选 D.失血性休克复选 E.肝性昏迷复选 正确答案: D 答案解析:略 ---------------------------------------------------------- ------------------------第2题------------------------ 题型: 【共用题干单选题】 大问: 男性,30岁。夜间发作性腹部烧灼样痛二月余,进食后能迅速缓解,昨起排柏油样便3次,今晨起床晕倒。体检:T 37℃,P 120次/min,R 24次/min,BP 80/50mmHg。

该病情的观察重点是 答案选项: A.生命体征复选 B.意识复选 C.尿量复选 D.瞳孔变化复选 E.四肢活动度复选 正确答案: A 答案解析:略 ---------------------------------------------------------- ------------------------第3题------------------------ 题型: 【共用题干单选题】 大问: 男性70岁。慢性咳嗽、咳痰12年,近两年来活动后出现气短,偶有踝部水肿,门诊以“慢性支气管炎合并慢性阻塞性肺气肿”收入院 小问: 护理评估时胸部阳性体征可表现为 答案选项: A.扁平胸复选 B.语颤减弱复选 C.语频增强复选 D.心浊音界扩大复选 E.呼气期短复选 正确答案: B

信号与系统课程设计报告

信号与系统课程设计报告 实验题目:信号的运算与处理 内容简介: 设计一个信号,对其进行信号运算和处理,利用Matlab仿真。 课设方式: 利用电子技术、电路理论和信号与系统的知识学习验证信号的运算和处理,如延时、相加、微分、抽样等。自已设计信号及运算方式,并利用Matlab仿真。 分析计算结果。 课程设计要求: 独立完成; 完成信号设计(任意信号均可)及其某种运算(任意运算均可,也可多做几种,或做组合运算)的验证; 学会利用Matlab仿真;提交课程设计报告。 例如: 设计一个信号为f(t)=3sin2t 对其做微分运算得到f/(t) , 用MATLAB 编程实现计算过程,画出f(t)和f/(t)

本次课程设计本人选的信号运算是: 设计一个信号为y1=y(x)=sin2x,对其作微分运算得到dy1,用MATLAB对其实现运算过程,后画出y1,dy1,y1+dy1的图像 实验步骤(操作过程) 1、 首先打开MATLAB软件,在其命令窗口直接输入以下程序,对y(x)进 行微分运算。得到dy1 clear >> syms x y1; >> y1=sin(2*x); >> dy1=diff(y1,'x') dy1 =2*cos(2*x) 运算过程如下图所示: 2、 接着便是对其进行验证,点击fire,新建一个文件,输入以下程序(绘制出y1=sin2x, dy1=2cos2x, 以及y1+ dy1=sin2x+2cos2x。的波形)

3、保存文件,后缀名为.m,随后按F5执行输出输出图形。实验结果如下图所示 、

结果分析 如图所示绿色波形为y1=sin2x,蓝色为dy1=2cos2x,红色波形为y1+dy1。仿真结果与运算结果一致。 实验心得体会(调试过程) 总的来说,这次课程设计难度并不是太高,而我选取的正玄信号也是较为简单常用的一种函数,对其进行微分运算之后,得到了余弦函数,其仿真结果波形也如上所示,与预期一致。在设计过程中,还是出现了几个小问题的,一个是变量的定义,之前没有定义x,直接取范围结果出错了,还有一个是注意各种函数的调用以及运算格式,还是希望能在之后再接再厉,掌握好matlab软件!(附上调试过程图片) 左边为文件、历史窗口,底下是命令窗口,最右下角为实验仿真波形,中间为运算程序,绘图画图程序。

通信原理课程设计报告(基于Matlab)

2DPSK调制与解调系统的仿真 设计原理 (1) 2DPSK信号原理 1.1 2DPSK信号原理 2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差,并规定:Φ=0表示0码,Φ=π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的,在接收端只能采用相干解调,它的时域波形图如图2.1所示。 图1.1 2DPSK信号 在这种绝对移相方式中,发送端是采用某一个相位作为基准,所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化,则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式,而采用相对移相方式。定义为本码元初相与前一码元初相之差,假设: →数字信息“0”; →数字信息“1”。 则数字信息序列与2DPSK信号的码元相位关系可举例表示如下: 数字信息: 1 0 1 1 0 1 1 1 0 1 DPSK信号相位:0

或 : 1.2 2DPSK 信号的调制原理 一般来说,2DPSK 信号有两种调试方法,即模拟调制法和键控法。2DPSK 信号的的模拟调制法框图如图1.2.1所示,其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。 图1.2.1 模拟调制法 2DPSK 信号的的键控调制法框图如图1.2.2所示,其中码变换的过程为将输入的基带信号差分,即变为它的相对码。选相开关作用为当输入为数字信息“0” 时接相位0,当输入数字信息为“1”时接pi 。 图1.2.2 键控法调制原理图 1.3 2DPSK 信号的解调原理 2DPSK 信号最常用的解调方法有两种,一种是极性比较和码变换法,另一种是差分相干解调法。 码变换 相乘 载波 s(t) e o (t)

信号与系统 MATLAB实验报告

《信号与系统》MATLAB实验报告 院系:专业: 年级:班号: :学号: 实验时间: 实验地点:

实验一 连续时间信号的表示及可视化 实验题目: )()(t t f δ=;)()(t t f ε=;at e t f =)((分别取00<>a a 及); )()(t R t f =;)()(t Sa t f ω=;)2()(ft Sin t f π=(分别画出不同周期个数的波形)。 解题分析: 以上各类连续函数,先运用t = t1: p:t2的命令定义时间围向量,然后调用对应的函数,建立f 与t 的关系,最后调用plot ()函数绘制图像,并用axis ()函数限制其坐标围。 实验程序: (1))()(t t f δ= t=-1:0.01:3 %设定时间变量t 的围及步长 f=dirac(t) %调用冲激函数dirac () plot(t,f) %用plot 函数绘制连续函数 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的围 (2))()(t t f ε= t=-1:0.01:3 %设定时间变量t 的围及步长 f=heaviside(t) %调用阶跃函数heaviside () plot(t,f) %用plot 函数绘制连续函数 title('f(t)=heaviside(t)') %用title 函数设置图形的名称 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的围 (3)at e t f =)( a=1时:

t=-5:0.01:5 %设定时间变量t 的围及步长 f=exp(t) %调用指数函数exp () plot(t,f) %用plot 函数绘制连续函数 title('f=exp(t)') %用title 函数设置图形的名称 axis([-5,5,-1,100]) %用axis 函数规定横纵坐标的围 a=2时: t=-5:0.01:5 f=exp(2*t) % 调用指数函数exp () plot(t,f) title('f=exp(2*t)') axis([-5,5,-1,100]) a=-2时: t=-5:0.01:5 f=exp(-2*t) plot(t,f) title('f=exp(-2*t)') axis([-5,5,-1,100]) (4))()(t R t f = t=-5:0.01:5 f=rectpuls(t,2) % 用rectpuls(t,a)表示门函数,默认以零点为中心,宽度为a plot(t,f) title('f=R(t)') axis([-5 5 -0.5 1.5]) (5))()(t Sa t f ω= ω=1时: t=-20:0.01:20 f=sin(t)./t % 调用正弦函数sin (),并用sin (t )./t 实现抽 样函数 plot(t,f)

2015产科护理案例分析题及答案解析

2015产科护理案例分析题及答案解析 天宇考王卫生资格高级职称考试题库包含:章节练习、综合复习题、模拟试卷、考前冲刺、历年真题等。 ------------------------第1题------------------------ 题型: 【单选题】 大问: 胃溃疡的发生原因 小问: 答案选项: A.保护因素消弱单选 B.胃壁血液供应丰富单选 C.胆汁反流单选 D.氧自由基单选 E.精神紧张单选 正确答案: A 答案解析:略 ---------------------------------------------------------- ------------------------第2题------------------------ 题型: 【单选题】 大问: 病人万某,因感染服磺胺药治疗,护士嘱其多饮水,其目的是 小问: 答案选项:

A.促进吸收单选 B.减少副作用单选 C.冲淡药味单选 D.防止在肾脏析出结晶单选 E.保护肝脏单选 正确答案: D 答案解析:略 ---------------------------------------------------------- ------------------------第3题------------------------ 题型: 【单选题】 大问: 洗发中,梳下的落发应置于( ) 小问: 答案选项: A.包好保存单选 B.垃圾袋中单选 C.盆中单选 D.污水桶中单选 E.焚烧单选 正确答案: B 答案解析:略 ---------------------------------------------------------- ------------------------第4题------------------------

数字信号处理课程设计报告

《数字信号处理》课程设计报告 设计题目: IIR滤波器的设计 专业: 班级: 姓名: 学号: 指导教师: 2010年月日

1、设计目的 1、掌握IIR 滤波器的参数选择及设计方法; 2、掌握IIR 滤波器的应用方法及应用效果; 3、提高Matlab 下的程序设计能力及综合应用能力。 4、了解语音信号的特点。 2、设计任务 1、学习并掌握课程设计实验平台的使用,了解实验平台的程序设计方法; 2、录制并观察一段语音信号的波形及频谱,确定滤波器的技术指标; 3、根据指标设计一个IIR 滤波器,得到该滤波器的系统响应和差分方程,并根据差分方程将所设计的滤波器应用于实验平台,编写相关的Matlab 程序; 4、使用实验平台处理语音信号,记录结果并进行分析。 3、设计内容 3.1设计步骤 1、学习使用实验平台,参见附录1。 2、使用录音机录制一段语音,保存为wav 格式,录音参数为:采样频率8000Hz、16bit、单声道、PCM 编码,如图1 所示。 图1 录音格式设置 在实验平台上打开此录音文件,观察并记录其波形及频谱(可以选择一段较为稳定的语音波形进行记录)。 3、根据信号的频谱确定滤波器的参数:通带截止频率Fp、通带衰减Rp、阻带截止频率Fs、阻带衰减Rs。 4、根据技术指标使用matlab 设计IIR 滤波器,得到系统函数及差分方程,并记录得到系统函数及差分方程,并记录其幅频响应图形和相频响应图形。要求设计 第 1页出的滤波器的阶数小于7,如果不能达到要求,需要调整技术指标。 5、记录滤波器的幅频响应和系统函数。在matlab 中,系统函数的表示公式为:

因此,必须记录系数向量a 和b。系数向量a 和b 的可以在Matlab 的工作空间(WorkSpace)中查看。 6、根据滤波器的系统函数推导出滤波器的差分方程。 7、将设计的滤波器应用到实验平台上。根据设计的滤波器的差分方程在实验平台下编写信号处理程序。根据运行结果记录处理前后的幅频响应的变化情况,并试听处理前后声音的变化,将结果记录,写入设计报告。 3.2实验程序 (1)Rs=40; Fs=1400; Rp=0.7; Fp=450; fs=8000; Wp=2*pi*Fp;Ws=2*pi*Fs; [N,Wn]=buttord(Wp,Ws,Rp,Rs,'s'); [b1,a1]=butter(N,Wn,'s'); [b,a]=bilinear(b1,a1,fs); [H,W]=freqz(b,a); figure; subplot(2,1,1);plot(W*fs/(2*pi),abs(H));grid on;title('频率响应'); xlabel('频率');ylabel('幅值');、 subplot(2,1,2); plot(W,angle(H));grid on;title('频率响应'); xlabel('相位(rad)');ylabel('相频特性'); 3.3实验结果(如图): N =5 Wn=6.2987e+003 第 2页

MATLAB课设报告

课程设计任务书 学生姓名:董航专业班级:电信1006班 指导教师:阙大顺,李景松工作单位:信息工程学院 课程设计名称:Matlab应用课程设计 课程设计题目:Matlab运算与应用设计5 初始条件: 1.Matlab6.5以上版本软件; 2.课程设计辅导资料:“Matlab语言基础及使用入门”、“Matlab及在电子信息课程中的应 用”、线性代数及相关书籍等; 3.先修课程:高等数学、线性代数、电路、Matlab应用实践及信号处理类相关课程等。 要求完成的主要任务:(包括课程设计工作量及其技术要求,以及说明书撰写等具体要求) 1.课程设计内容:根据指导老师给定的7套题目,按规定选择其中1套完成; 2.本课程设计统一技术要求:研读辅导资料对应章节,对选定的设计题目进行理论分析, 针对具体设计部分的原理分析、建模、必要的推导和可行性分析,画出程序设计框图,编写程序代码(含注释),上机调试运行程序,记录实验结果(含计算结果和图表),并对实验结果进行分析和总结。具体设计要求包括: ①初步了解Matlab、熟悉Matlab界面、进行简单操作; ②MATLAB的数值计算:创建矩阵矩阵运算、多项式运算、线性方程组、数值统计; ③基本绘图函数:plot, plot3, mesh, surf等,要求掌握以上绘图函数的用法、简单图形 标注、简单颜色设定等; ④使用文本编辑器编辑m文件,函数调用; ⑤能进行简单的信号处理Matlab编程; ⑥按要求参加课程设计实验演示和答辩等。 3.课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写,具体包括: ①目录; ②与设计题目相关的理论分析、归纳和总结; ③与设计内容相关的原理分析、建模、推导、可行性分析; ④程序设计框图、程序代码(含注释)、程序运行结果和图表、实验结果分析和总结; ⑤课程设计的心得体会(至少500字); ⑥参考文献(不少于5篇); ⑦其它必要内容等。 时间安排:1.5周(分散进行) 参考文献: [1](美)穆尔,高会生,刘童娜,李聪聪.MA TLAB实用教程(第二版) . 电子工业出版社,2010. [2]王正林,刘明.精通MATLAB(升级版) .电子工业出版社,2011. [3]陈杰. MA TLAB宝典(第3版) . 电子工业出版社,2011. [4]刘保柱,苏彦华,张宏林. MATLAB 7.0从入门到精通(修订版) . 人民邮电出版社,2010. 指导教师签名:年月日 系主任(或责任教师)签名:年月日

产科护理病例分析

王××,女,26岁,已婚,因停经55d,不规则阴道出血3d,剧烈右下腹痛2小时,病人既往月经规律,14岁初潮,5d/28~30d,末次月经2003年4月5日,婚育史:0-0-1-0。 体检:T37℃,P110次/分,R20次/分,Bp80/50mmHg,下腹压痛反跳痛明显,以右侧为甚,叩诊移动性浊音(+),妇检:阴道内有少量血液,子宫略大,软,后穹隆饱满,有触痛,宫颈举痛(+),子宫右侧可触及肿块,边界不清。 实验室检查:WBC7×109/L,Hb100g/L,尿妊娠试验(+) 1、根据病史资料,初步考虑该病人的疾病诊断是什么?有何依据?还需进行哪项检查? 2、病人主要存在哪些护理问题? 3、就病人现有的护理问题,你首先采取哪些护理措施? 4、此类妊娠并发症将如何处理? ● 2 某女,28岁,因“停经58日,腹痛.阴道出血一日”为主诉入院。患者平素月经 正常,末次月经2006年2月3日,停经后有恶心.呕吐等反应。停经40日,在我院检查尿妊娠试验阳性,同时B超检查提示“宫内妊娠5周”,昨日下午6时许同房后出现下腹痛,阵发性伴少许阴道流血,无血块.未见组织物排出,无畏寒.发热.头晕.乏力。 ●经带胎产史:初潮12岁,经期5天,周期30天,月经规则。25岁结婚,怀孕3次, 产育0次,2001年行药物流产1次,2003年行人工流产1次。白带无特殊。 ●专科检查:外阴正常,阴道通畅.内有中等量暗红色血污。宫颈光滑,宫口闭合, 未见组织物阻塞,抬举痛阴性。妇检子宫软,如孕8周大小,双附件未见明显包块。 ●实验室及辅助检查:血常规.尿常规正常。B超检查示宫内妊娠8周,活胎。双侧 附件未见异常。 1.该病诊断为何病? 2.写出该病的护理诊断及所应采取的护理措施 3.入院三天后,病人出现腹痛加剧,阴道出血量大如月经量,自诉有肉样物排出。妇科检查:外阴阴道大量血污,宫颈口见组织物阻塞,子宫小于孕周。 (1)此时应诊断为何病?(2)应采取何种手术处理。 26岁初产妇,胎儿娩出后无阴道流血,胎盘娩出后阴道流血不断,时多时少,1小时内出血量超过600ml,血压70/30mmhg,脉搏126次/分 ①此时出血的原因应考虑什么?应采取的紧急措施是什么? ②指出主要护理诊断二个,并指出相应的目标及措施。 某孕妇,36岁,G1P0,于规律宫缩前破水,现已临产17小时,宫口开大3cm,因精神过度紧张不能进食进水,不合作。宫缩每6-7分钟1次,持续20-30秒,宫口不能如期扩张,先露不能如期下降。 (1)、作出相应诊断及其依据为何? (2)、目前主要的护理诊断 (3)、针对护理诊断制定相应护理措施

信号分析课程设计报告书

信号分析课程设计 信号系统的时域分析 编程实现的卷积积分或卷积和 一、课程设计题目: 基于 MATLAB 的连续时间LTI 系统的时域分析 二、基本要求: ① 掌握连续时不变信号处理的基本概念、基本理论和基本方法; ② 学会 MATLAB 的使用,掌握 MATLAB 的程序设计方法; ③ 学会用 MATLAB 对信号进行分析和处理; ④ 编程实现卷积积分或卷积和,零输入响应,零状态响应; ⑤ 撰写课程设计论文,用信号处理基本理论分析结果。 三、设计方法与步骤: 一般的连续时间系统分析有以下几个步骤: ①求解系统的零输入响应; ②求解系统的零状态响应; ③求解系统的全响应; ④分析系统的卷积;⑤画出它们的图形. 下面以具体的微分方程为例说明利用MATLAB 软件分析系统的具体方法. 1.连续时间系统的零输入响应 描述n 阶线性时不变(LTI )连续系统的微分方程为: 已知y 及各阶导数的初始值为y(0),y (1)(0),… y (n-1)(0), 求系统的零输入响应。 建模 当LIT 系统的输入为零时,其零输入响应为微分方程的其次解(即令微分方程的等号右端为零),其形式为(设特征根均为单根) 其中p 1,p 2,…,p n 是特征方程a 1λ n +a 2λn-1+…+a n λ+a n =0的根,它们可以 用root(a)语句求得。各系数 由y 及其各阶导数的初始值来确定。对此有 1121111n n m n n m m n n m d y d y dy d u du a a a a y b b b u dt dt dt dt dt -++-++?????++=+????++1212()n p t p t p t n y t C e C e C e =++????+120n C C C y ++????+=11220 n n p C p C p C Dy ++????+=

matlab课程设计报告书

《计算机仿真及应用》课程设计报告书 学号:08057102,08057127 班级:自动化081 姓名陈婷,万嘉

目录 一、设计思想 二、设计步骤 三、调试过程 四、结果分析 五、心得体会 六、参考文献

选题一、 考虑如下图所示的电机拖动控制系统模型,该系统有双输入,给定输入)(t R 和负载输入)(t M 。 1、 编制MATLAB 程序推导出该系统的传递函数矩阵。 2、 若常系数增益为:C 1=Ka =Km =1,Kr =3,C2=0.8,Kb =1.5,时间常数T 1=5, T 2=0.5,绘制该系统的根轨迹、求出闭环零极点,分析系统的稳定性。若)(t R 和)(t M 分别为单位阶跃输入,绘制出该系统的阶跃响应图。(要求C 1,Ka ,Km ,Kr ,C2,Kb , T 1,T 2所有参数都是可调的) 一.设计思想 题目分析: 系统为双输入单输出系统,采用分开计算,再叠加。 要求参数均为可调,而matlb 中不能计算未赋值的函数,那么我们可以把参数设置为可输入变量,运行期间根据要求赋值。 设计思路: 使用append 命令连接系统框图。 选择‘参数=input('inputanumber:')’实现参数可调。 采用的方案: 将结构框图每条支路稍作简化,建立各条支路连接关系构造函数,运行得出相应的传递函数。 在得出传递函数的基础上,使用相应的指令求出系统闭环零极点、画出其根轨迹。 通过判断极点是否在左半平面来编程判断其系统是否稳定。 二.设计步骤 (1)将各模块的通路排序编号

(2)使用append命令实现各模块未连接的系统矩阵 (3)指定连接关系 (4)使用connect命令构造整个系统的模型 三.调试过程 出现问题分析及解决办法: 在调试过程出现很多平时不注意且不易寻找的问题,例如输入的逗号和分号在系统运行时不支持中文格式,这时需要将其全部换成英文格式,此类的程序错误需要细心。 在实现参数可调时初始是将其设为常量,再将其赋值进行系统运行,这样参数可调性差,后用‘参数=input('inputanumber:')’实现。 最后是在建立通路连接关系时需要细心。 四.结果分析 源代码: Syms C1 C2 Ka Kr Km Kb T1 T2 C1=input('inputanumber:') C2=input('inputanumber:') Ka=input('inputanumber:') Kr=input('inputanumber:') Km=input('inputanumber:') Kb=input('inputanumber:') T1=input('inputanumber:') T2=input('inputanumber:') G1=tf(C1,[0 1]); G2=tf(Ka*Kr,[0 1]); G3=tf(Km,[T1 1]); G4=tf(1,[T2 1]); G5=tf(1,[1 0]); G6=tf(-C2,1); G7=tf(-Kb,1); G8=tf(-1,1); Sys=append(G1,G2,G3,G4,G5,G6,G7,G8) Q=[1 0 0;2 1 6;3 2 7;4 3 8;5 4 0;6 5 0;7 4 0;8 0 0;]; INPUTS1=1; OUTPUTS=5; Ga=connect(Sys,Q,INPUTS1,OUTPUTS) INPUTS2=8; OUTPUTS=5; Gb=connect(Sys,Q,INPUTS2,OUTPUTS) rlocus(Ga)

MATLAB实验报告

实验一 名称:连续时间信号分析 姓名:王嘉琦 学号:201300800636 班级:通信二班 一、实验目的 (一)掌握使用Matlab 表示连续时间信号 1、学会运用Matlab 表示常用连续时间信号的方法 2、观察并熟悉常用信号的波形和特性 (二)掌握使用Matlab 进行连续时间信号的相关运算 1、学会运用Matlab 进行连续时间信号的时移、反褶和尺度变换 2、学会运用Matlab 进行连续时间信号微分、积分运算 3、学会运用Matlab 进行连续时间信号相加、相乘运算 4、学会运用Matlab 进行连续时间信号卷积运算 二、实验条件 Matlab 三、实验内容 1、利用Matlab 命令画出下列连续信号的波形图。 (1))4/3t (2cos π+ 代码: k=2;w=3;phi=pi/4; t=0:0.01:3; ft=k*cos(w*t+phi); plot(t,ft),grid on; axis([0,3,-2.2,2.2]) title('余弦信号')

(2) )t (u )e 2(t -- 代码: k=-1;a=-1; t=0:0.01:3; ft=2-k*exp(a*t); plot(t,ft),grid on axis([0,3,2,3]) title('指数信号')

(3))]2()(u )][t (cos 1[--+t u t π 代码: k=1;w=pi;phi=0; t=0:0.01:2; ft=1+k*cos(w*t+phi); plot(t,ft),grid on; axis([0,3,0,2]) title('余弦信号')

《妇产科学》病例分析题

《妇产科学》病例分析题 1、孕妇35 岁,G1P0 孕36 周,浮肿(),血压21.3/14.6kPa (160/110mmH蛋白尿( ),一天前出现头痛,恶心。入院诊断为重度妊高征,先兆子痫,( 1 ) 试写出此孕妇3 个主要的护理诊断。( 2)该孕妇目前应接受哪些护理措施。( 3) 目前用硫酸镁治疗,写出硫酸镁使用的注意事项。答案:( 1)知识缺乏:缺乏对妊高征处理的相关知识有母儿受伤的危险:与子痫发作时病人意识丧失有关焦虑:与担心高血压及其对母儿的影响有关( 2) 1 、卧床休息,以左侧卧位为宜2 、指 导病人摄入足够的水和富含纤维素的食品,及足够蛋白质的摄入3 、硫酸镁的用药护理4 、重度妊高征患者,保持病情稳定、预防子痫发生,为分娩做好准备①将病人安排在安静、光线较暗的病室,医护活动尽量集中;②准备下列物品:呼叫 器、放好床档防坠床、急救药物用物、产包。( 3)硫酸镁的治疗浓度和中毒浓度相近,故在使用过程中严密观察其毒性作用,控制入量,滴速以1g/h ,不超过2g/h。同时监测以下指标:1、膝腱反射必须存在;2、呼吸不少于16次/分钟; 3、尿量不少于25ml/h ;尿少提示肾脏排泄功能受到抑制,镁离子蓄积,故应随时备好10%葡萄糖酸钙,及时给予解毒。2、一孕妇,26岁,孕27W(夜晚无异常否认性生活)晨起时发一内裤及床上有鲜血,站立时也随着出血感,立即来院,入院后给予保守治疗,问病史曾有2 次怀孕,自然流产,强烈要求保住此胎,胎儿安好,给母亲对症,支持期待疗法,阴道间断、少量出血,孕32W 晚2Am突然阴道 大量出血,蹲痰盂,500mL仍有活动出血;请作出诊断(相应的检查)及处理? (12分)诊断:前置胎盘答案:进一步检查:(1)观察生面体征。(2)B超观察胎儿、胎盘、羊水。( 3)胎监监护胎儿处理:( 1)促进胎儿肺成熟;( 2)绝对卧床、左侧卧位、吸氧;( 3)静脉给予补液,立即准备剖宫产手术,结束妊娠。 3、某患者,LMP:2004年3月10日,停经50天时,阴道流血,量少,有早孕反应,于2004年5月3日入院。无痛或有轻微下腹疼痛,伴腰痛及下坠感。妇科检查子宫颈口闭,子宫大小与停经月份符合,B超示:宫内可见一妊娠囊,有胚芽和 胎心音波动。请给该孕妇作出诊断、治疗。( 8分)答案:诊断:先兆流产治疗: ( 1)卧床休息,吸氧;( 2)营养饮食;( 3)禁止性生活( 4)黄体酮20mg im QD;VitE20mg tid (5)观察胎儿宫内情况。4:女性,15岁,13岁初潮,近3个月月经紊乱,经量较多,肛诊:子宫正常大小,双侧附件区软,未触及包块。 (1)该患最可能的诊断是什么?(2)拟给该患调整月经周期,宜 采用何方法,具体如何应用?答案:( 1 )无排卵性功能失调性子宫出血( 5分)( 2)雌孕激素序贯疗法( 3分)具体方案( 2分)已稀雌酚,于出血第5日起,每晚1 次,连服20日至服药第11 日,每日加用黄体酮10mg肌注(或安宫黄体酮8~10mg口服),两药同时用完,停药后3~7日出血,于出血第5日重复用药,一般连用3个周期。5:患者,女,30岁。习惯性流产4 次,足月分娩1 次,畸形新生儿死亡。产后半年再孕入院保胎治疗,孕期

数字信号课程设计报告

一、语音信号去噪处理 1.设计要求: (1)在windows系统下的录音机录制一段1s左右的语音信号作为原声信号,在MATLAB软件平台下,利用函数wavread对语音信号进行采样,记住采样频率和采样点数; (2)画出语音信号的时域波形,对采样后的语音进行fft变换,得到信号的频谱特性;对语音信号分别加入正弦噪声和白噪声,画出加噪信号的时域波形和频谱图; (3)根据对加噪语音信号谱分析结果,确定滤除噪声滤波器的技术指标,设计合适的数字滤波器,并画出滤波器的频域响应; (4)用所设计的滤波器对加噪的信号进行滤波,在同一个窗口画出滤波前后信号的时域图和频谱图,对滤波前后的信号进行对比,分析信号变化; (5)利用sound(x)回放语音信号,验证设计效果。 2.设计步骤: (1)找到7s的语音信号,利用函数wavread对语音信号进行信号读取;(2)计算样本时刻和频谱图的频率,并进行N+1点FFT变换; (3)加噪声为5000Hz的正弦信号正弦噪声,采用awgn函数加信噪比为10的高斯白噪声; (4)设计滤波器; (5)绘出相应的时域、频域图; (6)利用sound函数进行原始信号的语音播放,加噪声音播放,以及滤波之后的语言播放。 3.设计实现: (1)时域图与频谱图(加正弦) 录入原始信号的时域图: 加入正弦信号后的时域图:

滤波后的时域图: 录入原始信号的频域图: 加入正弦信号后的频率图: 滤波后的频域图: 采用巴斯低通滤波器滤除正弦波:

(2)具体代码实现: [x,fs,bits]=wavread('E:\mcpass.wav');%原信号 n=size(x,1); %提取采样信号的长度 t=(0:length(x)-1)/fs; %计算样本时刻 f=fs*(0:(n+1)/2-1)/n+1; %计算频域图的频率 X=fft(x,n+1); %进行N+1点FFT变换 ts=0:1/fs:(size(x)-1)/fs; %将所加噪声信号的点数调整到与原始信号相同 s=x+0.05*sin(2*pi*5000*ts)'; %加噪声为5000Hz的正弦信号正弦噪声 S=fft(s,n+1); %加正弦噪声后的频域 %正弦滤波 wp=2000/fs*2*pi; %2000为通带截止频率 ws=3000/fs*2*pi; %3000为阻带下限截止频率 Rp=4; %通带波纹 Rs=25; %阻带波纹 T=1/fs;Fs=1/T; %定义采样间隔 Wp=2/T*tan(wp/2); %计算对应的数字频率 Ws=2/T*tan(ws/2); [N,wn]=buttord(Wp,Ws,Rp,Rs,'s'); %计算滤波器介数和截止频率 [c,d]=butter(N,wn,'s'); %计算滤波器系统函数分子分母系数 [B,A]=bilinear(c,d,Fs); %双线性变换得到数字滤波器系统函数分子分母系数[Hb,Wc]=freqz(B,A); sf=filter(B,A,s); %对加噪信号进行滤波 Sf=fft(sf,n+1); %对滤波后进行N+1点FFT变换 %绘图部分 figure(3); plot(fs*Wc/(2*pi),20*log10(abs(Hb)));title('巴斯低通滤波器频域响应图'); xlabel('频率(Hz)'); ylabel('幅度'); figure(1); subplot(3,1,1); plot(t,x);title('原信号时域') xlabel('时间(s)'); ylabel('幅度'); figure(2); subplot(3,1,1); plot(f,abs(X(1:(n+1)/2)));title('原信号频域') xlabel('频率(Hz)'); ylabel('幅度'); figure(1); subplot(3,1,2); plot(t,s);title('加正弦信号后的时域') xlabel('时间(s)'); ylabel('幅度');

matlab音频降噪课程设计报告

matlab音频降噪课程设计报告

燕山大学 医学软件课程设计说明书 题目:基于MATLAB巴特沃斯滤波器的音频去噪的GUI设计 学院(系):电气工程学院 年级专业: 13级生物医学工程 2 班 学号: 130103040041 学生姓名:魏鑫 指导教师:许全盛 1

院(系):电气工程学院基层教学单位:生物医学工程系 学号130103040041 学生 姓名 魏鑫 专业(班 级) 13级生 物医学 工程2 班 设计 题目 基于MATLAB音频去噪的GUI设计设 计 技术参数通带截止频率fp=2700;阻带截止频率fs=3000;采样频率FS=48000; 通带衰减不大于1dB;阻带衰减不小于10dB; 设计要求1.实现用MATLAB导入音频; 2.对音频进行频谱分析; 3.设计滤波器去噪并对含噪信号进行滤 2

波并进行功率谱分析; 4.设计能实现上述功能的GUI; 工作量1.完成音频录入及频谱分析相关程序的编写与调试; 2.设计滤波器去噪; 3.用MATLAB软件做GUI界面的设计; 工作计划11.21-11.24 MATLAB软件中GUIDE 工具箱的使用 11.25-11.29 各处理算法模块的编程实现 11.30-12.1 整体程序联调 12.2 撰写课程设计说明书,答辩 参考资料 1. 陈怀琛吴大正 MATLAB及在电子信息课程中的应用[M] 北京电子工业出版社 2006. 章节2.4; 2. 陈亚勇 MATLAB信号处理详解[M] 北京:人民邮电出版社 2000. 第十 3

章; 3.张康刘雅基于Matlab的巴特沃斯 数字低通滤波器的设计[J] 计算机与现代化 2007年 12期 98-100页 指导 教师签字许全盛 基层教学单 位主任签字 彭勇 目录 一、设计目的意义 (1) 1.1绪论 (1) 1.2设计目的 (1) 1.3意义 (1) 二、设计内容 (2) 2.1 设计原理 (2) 2.2 设计内容 (2) 三、设计过程及结果分析 (3) 3.1 设计步骤 (3) 4

相关主题