第十章 MATLAB在电子信息课程中的应用

matlab

第十章 MTALBA在电子息课 信程中的应用10.1 MTALAB在路电中应的 1用02 MA.TAB在L数字信处理中的应用号

matlab

0. 1AMTABL在电中路应的例10.用 电1阻路电的求 如图1解.10所的电路，示 已：知R1=2,R2 =4, R3=1 2R4,=4 ,R5 =1 ,26=R4 ,7=R2

。图10. 110例.的电路图1

()1如知us已1=V0,求3i,4uu,;7( 2如)已知u46V,=求s,u3,ui。7

matlab

matlab

模：建用孔法网，照图按 建模用网孔：法，按照10图1.列出可网方孔 可列程网孔出方程( 1R+ R 2 +R 3)ia Ri3b =s u R3ai (+ 3 + RR +4 5 )ibR 5Ric= 0 R ib5 + (R + 56 + RR )ic7= 0

方程组写把成阵形矩式：A为X=B把 方程组写矩成阵式为： R形 1 R2 ++ R3 R 3 0R 3R +3 R4 R5 R5+ i a 1 R 5 bi = 0 usR 5 +6 R+R7 c i 0 0 (

matlab

)1us=1令V,0i3由i=-ai,b4=R4u*biu,=R7*7i即c可 到问得题1)的解。 (2(由)路电的性性线，质令可 i3:k=*1s,uu=k42us,u*7k=3u*。s可 由问以题(1)结的得到比例系数果 k=i13u/s,2=k4u/usk3=u,/u7s suu=4k/2i3=k,*us1k1*=4uk/,u7=2k*3s=uk*u4/k23

matlab

AMTLA程B序R =21R2=4;R;=13;24=R4;R=52;R164=;R7=2 ;dislpy(a'解题问()'1) 解问题%1 a11R=1+2+RR3a1;2=-R;313a0=;a21= R3;a-22R=3R+4+R5;23=aR5-; 31=0;a32=-Ra5;3a=R53+R+R76 A=[a1;1,1a,2a31;a1,2a2,a22;a313,a23a33,; ]b11;b2=0;b==3; B0[b=1b;;b23] u;s=npuit('su')= %,入输(解)的已知1件条I =AB*u\; si=Ia()1i;=Ib2)(;ci=(I)3; 3i=aiib,u4-R4=*bi,7=Ru*i7c

matlab

> ans>= 问解题 ( 解)题(问1 u)s =u = s1 i0 =3 03.047 4 u= .2222 27 u =.7004

matlab

7 7u=k2*3u42 /k

2ans= 解题( ) 解问问题() 给2u42定6 给定=dis lpa(y解'题(2)'问)us2 = 2 7.000 u420i=punt'给(u定4=2') 32 i =1=k3iu/sk;2u=/4usk3=u7/;s; u.0001 0s2=uu42k2/,3i2=1*u42k/2 ku2 = 72

matlab

例 01.2阶低一通的率响频，应以U为c响，应 求频率响应函，画数其幅频响应(幅频出特性 )H|(ω)|j和频的响应(相频相特性)θ ω(。)US

/1jωcUc

matlab

频 响率函数(即应弦稳态正络函数)定 网义为应响输()出量向与激励(输Y)向入 量F比，即之(Hωj=|) H(ωje)jθω)(=Y/|F 。用MA TLBA中ab的s(H)an和gelH()句可语以直 接计算幅频应响相频响应和

。

matlab

建 模分用公式可求压得频响率函数 应 1 & Ujc ω1C ( Hj ω ) = == &ω U SU s R j 1 1 j+ ωCω 1c ω = 为截止c频角率 C

1/jRcωUc 设无量纲频率ω=ωwωc/0=0.,,2,4,…出画频响幅应 相及频应响

。

matlab

leca,frormt comaacp tww0:=0.:2; 4H=1.(/1j+w*); wfigru(e)1 sbplout(2,,11)p,ol(ww,atbs())H grdixl,bal('wew',y)abel(la'sb()'H s)bpuol(2t1,2),p,lot(w,wanlg(He)) grdi,laxeb(lw''w)yla,bl('enalegH()')f igure()2subpl t(2,o11),se,ilmgx(wow,0*2log01(bs(a)H) )grid,laxbel'(ww)'ylab,e(l分贝'') sbuplt(o2,,21,p)ot(lwwa,nlegH()) rig,dxlabl('wwe' ylabel)'(agle(nH'))

matlab

matlab

0.2 M1ALATB在号信系与统中的用 应例1 .30 连续号信的连续信号的M ATLB描A 描述 述(1)单冲位函激数 ： )1 x 1 (t) = δ( tt 1 ) = 0其 余

t 1< <t 1t

+2)(单位跃阶数函：)单 位跃阶函：数1 x 2t ) (= (ut t ) 1=

t01< t t1 + < t 0<

%

matlab

1)单(冲位击应 响%在t(t1≤01t≤ft)有处一续持时为间td,积为1面 脉冲的信号，其时余间均零。为cl ae,r0t=0tf;=5;t1= 1 dt;0.=05 t=[;t0:td:tf;]st=legnhtt(; n)1f=olor((1-t0)t/td;) %求t1对的应本序样 号x1=erosz(,1ts) x1(n;1)=1/td sta;ir(t,xs1,)gidr no,xasi(0,[50,,2]2

)

matlab

(%2)单阶跃信号位% 信从t号0到tf在t1(t,0≤1ttf)前为≤,0t1处到有一 跃，变后以1为c lar,et=0;tf05; =t11; d==0t05; .=tt0[:dt:ft;]tslen=tg(ht;)n1=floo r((1-t0)t/td; x)=2ze[ors(,1n1,)nose1,st-n()] ; s1atrsi(,tx2,)rgi odn,aixs[(05,0,,1.]1)