应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 各位帮忙处理一个信号解卷积的问题
142/2返回列表上一页12
查看: 1757  |  回复: 13
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

我试了下，算完之后再对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')
赞一下 回复此楼
showmethemoney
11楼2012-07-06 15:01:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuangqingx

铁杆木虫 (著名写手)

宇宙超级无敌小菜鸟
引用回帖:
11楼: Originally posted by csgt0 at 2012-07-06 15:01:08
我试了下，算完之后再对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( ...

算完应该是和你最初贴的那个图相似，但是是正峰，不是这样的周期性的峰
赞一下 回复此楼
低调做人
12楼2012-07-06 15:58:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

引用回帖:
12楼: Originally posted by shuangqingx at 2012-07-06 15:58:48
算完应该是和你最初贴的那个图相似，但是是正峰，不是这样的周期性的峰...

等大牛解决吧。
我也不知道对不对，只是用这个方法计算得到的变换系数，再与初始输入进行卷积后得到的结果与初始输出是一样的。
赞一下 回复此楼
showmethemoney
13楼2012-07-06 16:03:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuh668

新虫 (初入文坛)
引用回帖:
5楼: Originally posted by shuangqingx at 2012-07-06 11:29:40
如果用卷积定理来做呢？就是将两个信号先傅里叶变换，相除，之后再反变换回来能行吗？
我刚刚请教别人给我试了一下，结果还是出现震荡的趋势。。。
请教您一下，您知道问题出在哪里了吗？...

按照你说的样子，估计是病态的问题。我也碰到这个问题，还没解决
赞一下 回复此楼
14楼2015-12-31 10:28:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shuangqingx 的主题更新
142/2返回列表上一页12
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭