请问如何在matlab中求传递函数

fenfenyang

已有输入和输出的文本数据，怎么得到传递函数？恳求各位给予指点...

ch_j1985

原帖由 fenfenyang 于 2008-6-5 09:02 发表

                               
登录/注册后可看大图

已有输入和输出的文本数据，怎么得到传递函数？恳求各位给予指点...

LZ能否把背景知识也介绍一下？
是否可以用拟合？

nsin

期待解答，或者谁有这方面的实例

nsin

输入输出的数据均为原始的时间信号，找两者之间的传递函数，直接按定义做？？先做fft，然后相除，这在matlab里怎么实现

[ 本帖最后由 nsin 于 2008-7-14 16:52 编辑 ]

sogooda

这样？
X=fft(x);
Y=fft(y);
f=X./Y;

nsin

就是不知道啊，感觉楼上这么做不太好，同疑问

ChaChing

xf=fft(x); yf=fft(y);
Pxx=abs(xf).*abs(xf); Pyy=abs(yf).*abs(yf); Pxy=yf.*conj(xf); Pyx=xf.*conj(yf);
Txy1=Pxy./Pxx; Txy2 = Pyy ./ Pyx;

hadesloveln

这个行不行啊？希望有解答

ChaChing

xf=fft(x); yf=fft(y);
Pxx=abs(xf).*abs(xf); Pyy=abs(yf).*abs(yf); Pxy=yf.*conj(xf); Pyx=xf.*conj(yf);
Txy1=Pxy./Pxx; Txy2 = Pyy ./ Pyx;

新版的spectrum没再仔细看过, 但我可以确定matlab中的spectrum也是如此用(7F)!

shengshengchina

哈哈，线性系统中的问题，这种问题用专业命令即可，自己找找吧。

yuxinpeng

回复 7 # ChaChing 的帖子

xf=fft(x); yf=fft(y);
Pxx=abs(xf).*abs(xf);
Pyy=abs(yf).*abs(yf);
Pxy=yf.*conj(xf);
Pyx=xf.*conj(yf);
Txy1=Pxy./Pxx;
Txy2 = Pyy ./ Pyx;

请问，能够加上注释吗？最后的传递函数是哪个？先谢过！

ChaChing

回复 11 # yuxinpeng 的帖子

LS的问题已经并非matlab使用问题, 应该是DSP专业问题了:@)
建议有空去逛逛信号版块, 应该可以学习更多
个人水平/时间有限, 仅大约记得Pxx/Pyy是APSD(auto-spectrum), Pxy/Pyx是CPSD(cross-spectrum)
至於Txy1/Txy2两者, 若输出/输入皆无杂讯干扰的话, 理论上则无差异, 一般个人大都使用第一式
没再查证过, 希望没记错:@)

yuxinpeng

回复 12 # ChaChing 的帖子

我的问题是已知输入和传递函数，想求输入，但编程的时候遇到问题。
如果前辈有兴趣，还请移驾下面的帖子给与指点。不胜感激！
http://forum.vibunion.com/thread-108921-1-1.html

随心就动

ChaChing 发表于 2009-9-25 11:28

                               
登录/注册后可看大图

xf=fft(x); yf=fft(y);
Pxx=abs(xf).*abs(xf); Pyy=abs(yf).*abs(yf); Pxy=yf.*conj(xf); Pyx=xf.*conj(yf ...

这样求的话误差较大，并且把相位信息丢了，怎么求相位传递函数？

ChaChing

随心就动 发表于 2013-6-11 16:42

                               
登录/注册后可看大图

这样求的话误差较大，并且把相位信息丢了，怎么求相位传递函数？

Txy1=Pxy./Pxx; Txy2 = Pyy ./ Pyx;
得到的传递函数为复数, 其中不是就有幅值及相位信息!?
个人水平有限, 不清楚LS的"误差较大，并且把相位信息丢了"的含意?
可拨空说明下吗?