24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

shuangqingx

铁杆木虫 (著名写手)

宇宙超级无敌小菜鸟

应助: 10 (幼儿园)
金币: 5576.5
散金: 4056
红花: 32
帖子: 2006
在线: 228.5小时
虫号: 472069
注册: 2007-12-04
性别: GG
专业: 化学工程及工业化学

[求助] 各位帮忙处理一个信号解卷积的问题

我的实验中利用的是一个方波输入信号f(x)，假设系统的响应函数为g(x)，则出口处的输出信号h(x)=f(x)*g(x)。现在我测得了系统出口的输出信号h(x)，想利用解卷积操作反算响应函数g(x)。因为我学的化工专业，所以对这方面不大懂；附件里我上传了输入信号和输出信号，请各位帮忙处理一下，非常感谢。

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 小木虫求助解卷积.xls

2012-07-05 17:48:32, 715.5 K

» 猜你喜欢

低调做人

1楼 2012-07-05 17:50:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
应助: 367 (硕士)
贵宾: 0.457
金币: 7812
散金: 8897
红花: 20
沙发: 7
帖子: 2810
在线: 515.9小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

【答案】应助回帖

我试了下，算完之后再对u卷积得到v1，新得到的v1前一半和v几乎一样，
sum(v1(1:length(u))-v')的结果为 1.0842e-014

clc
clear
load a.mat
p(1)=v(1)/u(1);
for i=2:length(u)
k=p(1:i-1)*flipud(u(2:i));
p(i)=(v(i)-k)/u(1);
end
v1=conv(u,p);
sum(v1(1:length(u))-v')

plot(1:length(u),u,1:length(u),v,1:length(u),p,'r')
figure
plot(1:length(u),v,1:length(u),v1(1:length(u)),'r')