求解电路传输响应函数

青蓝忘情水

那本书我在网上找到了，要定购，不过我已经找到其他相关得网站了，多谢提供宝贵意见跟信息

熊猫羊

可以直接下载,pdf版

青蓝忘情水

我用MATLAB计算这个阶跃响应得输出函数
SYS=tf([1],[L*C C*R 1]); step(tf);为什么我得OUTPUTSIGNAL出不来呀:mad: 不知道哪里出问题了

青蓝忘情水

还是这个问题，现在要求我计算输出端电流
我的具体程序如下：
R=1; L=1e-3; C=1e-6;
t=0:.001:1;
lamda1=-R/(2*L)+sqrt(R/(2*L)-1/(L*C));
lamda2=-R/(2*L)-sqrt(R/(2*L)-1/(L*C));
ia=1/(L*(lamda1-lamda2))*(exp(lamda1*t)-exp(lamda2*t));
subplot(3,1,1)
plot(1000*t(1:500),abs(ia(1:500)))
title('Ausgangsstrom in Zeitbereich')
xlabel('Zeit(milisekonde)')
ylabel('Ausgangsstrom(Ampêr)')
grid on

Ia=fft(ia,512);
f = 1e+6*(0:255)/512;
subplot(3,1,2)
plot(f,abs(Ia(1:256)))
title('Frequenzbebeich')
xlabel('Frequency(Hz)')
ylabel('|Y(f)|')
grid on

y=ifft(fft(ia),512);
subplot(3,1,3)
plot(abs(y))
title('Zeitbebeich')
xlabel('Zeit')
ylabel('Strom');
grid on
请高手看看，为什么我的结果跟我用别的软件模拟电路出的结果不一致呢？:(
应该在频域内幅值是1.6^(-3)，我的是0.025，请高手帮我看看，急

eisenstange

f = 1e+6*(0:255)/512;

你Sweep的区间是 从0 到 500K赫兹么？你电路的谐振频点在多少频率？

青蓝忘情水

是的，因为用另一个软件模拟只能从10kHz开始，我选择了10K-1MHz,所以我这里选到500K。谐振频率？我没有考虑呀，内组数值也是随便带入的难道还要考虑谐振频率？
我开始时算输出电压
clc;
clf;
R=100; L=10^(-3); C=10^(-6);
s=tf('s');
H=1/(1+s*C*R+s^2*C*L);
h=tf([0 1],[C*L R*C 1]);
ua=step(h);
subplot(3,1,1)
plot(ua)
grid on
title('Sprungantwort')
xlabel('time')
ylabel('Spannung')

Ua=fft(ua,512);
Y= Ua.* conj(Ua) / 512;
f = 1000*(0:255)/512;
subplot(3,1,2)
plot(f,Y(1:256)), grid on
title('Frequenzrebeich')
xlabel('frequency (Hz)')
y=ifft(Ua,512);
subplot(3,1,3)
plot(abs(y)),grid on
title('IFFT')
xlabel('time')
ylabel('Spannung')
经过ifft为什么在大于100时突然到零？
因为输出电压变化不大，所以教授要去我plot输出电流，我就做了上边的内些改动，不知道是否正确

[ 本帖最后由 青蓝忘情水 于 2006-12-7 16:46 编辑 ]

青蓝忘情水

Y= Ua.* conj(Ua) / 512;
这个是算功率谱的吧？我想问问可不可以省略这一行命令？

eisenstange

你两个软件 Sweep的范围要保持一致，至少上限要一致，因为S=jw，也就是说你的传输函数和频率相关，谐振可以先不考虑。你把Matlab的上限调到1M试试看。

青蓝忘情水

谢谢楼上的及时回复,1M(1e+6)?那我的抽样点要多大呀?

eisenstange

这个你要看FFT函数的说明了，应该对于数字信号计算，是取主周期，你需要自己看帮助。