ÐźÅÓëϵͳʵÑéÖ¸µ¼Êé(2008528) ÏÂÔØ±¾ÎÄ

ͼ4£­1 Ë¥¼õµç·

¼ÆËãÈçÓÒ£º

H(?)?U0(?)?Ui(?)R2j?C21R2?j?C2R1R2j?C1j?C2?11R1?R2?j?C1j?C2R21?j?R2C2R1R2?1?j?R1C11?j?R2C2

£½ £¨6-5£©

Èç¹û R1C1?R2C2 Ôò H(?)?R2Êdz£Êý£¬?(?)?0 £¨6-6£©

R2?R1 ʽ£¨6-6£©Âú×ãÎÞÊ§Õæ´«ÊäÌõ¼þ¡£

ËÄ¡¢ÊµÑéÒªÇó

£¨1£©»æÖƸ÷ÖÖÊäÈëÐźÅÊ§ÕæÌõ¼þϵÄÊäÈëÊä³öÐźţ¨ÖÁÉÙÈýÖÖ£©¡£ £¨2£©»æÖƸ÷ÖÖÊäÈëÐźÅÎÞÊ§ÕæÌõ¼þϵÄÊäÈëÊä³öÐźţ¨ÖÁÉÙÈýÖÖ£©¡£ £¨3£©±àÖÆ³öÍêÕûµÄʵÑé³ÌÐò£¬½øÐÐÑéÖ¤£¬»æÖÆÂ˲¨Æ÷µÄƵÂÊÏìÓ¦ÇúÏߣ¬ÐγÉʵÑ鱨¸æ¡£ ½â£º £¨1£©

R1=input('µç×èR1=') R2=input('µç×èR2=') C1=input('µçÈÝC1=') C2=input('µçÈÝC2=') syms t W;

x1=cos(2*pi*t); x2=exp(-2*abs(t));

x3=2*cos(2*pi*t)+3*sin(2*pi*1.5*t); F1=fourier(x1); F2=fourier(x2); F3=fourier(x3);

H1=R2/(1+i*W*R2*C2); H2=R1/(1+i*W*R1*C1); H=H1/(H2+H1); R1=H*F1; R2=H*F2;

13

R3=H*F3; f1=ifourier(R1) f2=ifourier(R2) f3=ifourier(R3)

subplot(321);ezplot(x1); subplot(322);ezplot(f1); subplot(323);ezplot(x2); subplot(324);ezplot(f2); subplot(325);ezplot(x3); subplot(326);ezplot(f3);

Ö´ÐкóÊäÈë²ÎÊý

µç×èR1=2000 R1 =

2000 µç×èR2=1000 R2 =

1000 µçÈÝC1=0.01 C1 =

0.0100 µçÈÝC2=0.01 C2 =

0.0100

µÃµ½Èçͼ6.2ËùʾµÄÈýÖÖÐźÅÔÚÊ§ÕæµÄÇé¿öϵÄÊäÈëÊä³öÐźźÍÊä³öÐźŵıí´ïʽ

f1 =

cos(2*pi*x)*(1+20*i*W)/(3+40*i*W) f2 =

(1+20*i*W)*(exp(2*x)*heaviside(-x)+exp(-2*x)*heaviside(x))/(3+40*i*W) f3 =

(3*sin(3*pi*x)+2*cos(2*pi*x))*(1+20*i*W)/(3+40*i*W)

14

ͼ4£­2 ÈýÖÖÐźÅÔÚÊ§ÕæµÄÇé¿öϵÄÊäÈëÊä³öÐźÅ

£¨2£©

R1=input('µç×èR1=') R2=input('µç×èR2=') C1=input('µçÈÝC1=') C2=input('µçÈÝC2=') syms t W;

x1=cos(2*pi*t); x2=exp(-2*abs(t));

x3=2*cos(2*pi*t)+3*sin(2*pi*1.5*t); F1=fourier(x1); F2=fourier(x2); F3=fourier(x3);

H1=R2/(1+i*W*R2*C2); H2=R1/(1+i*W*R1*C1); H=H1/(H2+H1); R1=H*F1; R2=H*F2; R3=H*F3; f1=ifourier(R1) f2=ifourier(R2) f3=ifourier(R3)

subplot(321);ezplot(x1); subplot(322);ezplot(f1); subplot(323);ezplot(x2);

15

subplot(324);ezplot(f2); subplot(325);ezplot(x3); subplot(326);ezplot(f3);

Ö´ÐкóÊäÈë²ÎÊý

µç×èR1=2000 R1 =

2000 µç×èR2=1000 R2 =

1000 µçÈÝC1=0.01 C1 =

0.0100 µçÈÝC2=0.02 C2 =

0.0200

µÃµ½Èçͼ6.3ËùʾµÄÈýÖÖÐźÅÔÚÊ§ÕæµÄÇé¿öϵÄÊäÈëÊä³öÐźźÍÊä³öÐźŵıí´ïʽf1 =

1/3*cos(2*pi*x) f2 =

1/3*exp(2*x)*heaviside(-x)+1/3*exp(-2*x)*heaviside(x) f3 =

sin(3*pi*x)+2/3*cos(2*pi*x)

ͼ4£­3 ÈýÖÖÐźÅÔÚÎÞÊ§ÕæµÄÇé¿öϵÄÊäÈëÊä³öÐźÅ

£¨3£©

clf;

R1=input('µç×èR1=') R2=input('µç×èR2=') C1=input('µçÈÝC1=') C2=input('µçÈÝC2=')

16