matlab无限长序列卷积,怎样求未知长度序列的卷积

程序是用fortran写的，但这里没有单独的板块，所以发到这里了，其实这不是一个程序语言的问题，用matlab编的话也会出现这个问题。

关于卷积，怎样求两个未知长度序列的卷积，就是不声明数组，只保存当前步数据

我们一般求卷积(非fft方法)的程序是：

这里求的是两组已知长度的序列x(n)和h(m)的卷积，卷积结果由y表示

dimension x(0-1),h(0-1),y(0-1)

l=m+n-1

Ierror=0

if(l.ge.(m+n-1))goto 1

Ierror=1

return

1 do 2 i=n,L-1 !由于y长度是m+n-1，所以x和h分别在n和m之外的部分用零

!补足

2 x(i)=0.0

do 3 i=m,L-1

3 h(i)=0.0

do 11 k=0,L-1

sum=0.0

do 12 i=0,k

sum=sum+x(i)*h(k-i)

12 continue

y(k)=sum

write(111,*) y(k)

11 continue

return

end

我们都知道，卷级的核心就是翻转，成绩，叠加，也就是这一步：

sum=sum+x(i)*h(k-i)

我现在想知道的是，可不可以不定义数组dimension x(0-1),h(0-1),y(0-1)，

也就是x或h的长度不知道，只知道要计算到某一时刻t就停止(这里就当成计算两组时间序列的卷积)，只保存当前步骤的卷积数值(也就是这里的sum)，把当前步骤的卷积值写到文本文件里就行了，比如下面这样(但这个结果不对，想问一下大家这个应该怎么写呢？)

program test_conv

!计算两组数据 ht和ft的卷积，ht和ft的具体取值如下面程序所写：

!ht=exp(-time)等等

time=0.0

do i=1,600

ht=exp(-time)

if(i>=1 .and.i<=101) then

ft=1.

else

ft=0.

end if

ans=ans+ht*ft*0.01 !现在计算的是相关，不是卷积，因为没有翻转

! 我感觉错误还是出现在这一步，但不知道应该怎么写

!才能跟上一个程序的结果一样呢？

! 然后在这里把这一步得到的ans，也就是卷积结果写入文本

time=time+0.01

end do

end

我的意思就是不把ht和ft定义成数组，怎样来求ht和ft的卷积？ (这里ft和ht的例子好像给的不好) 好像类似于叫动态卷积？

我写的有点多，但想尽量写得详细一写大家才能弄明白我的问题是什么意思吧？如果还有写的不清楚的地方，大家尽管问我。我被这个问题困扰很久了非常感谢

matlab无限长序列卷积,怎样求未知长度序列的卷积相关推荐

matlab无限长一维原子链,固体物理 03-03一维双原子链
<固体物理 03-03一维双原子链>由会员分享,可在线阅读,更多相关<固体物理 03-03一维双原子链(35页珍藏版)>请在人人文库网上搜索. 1.固体物理Solid Stat ...
【算法】未知长度序列等概率采样
问题是这样的,给一个序列,长度未知,要求从中随机等概率采样出m个元素. 最直观的解法是首先遍历序列,得到序列长度n,然后以m/n的概率从中随机取出m个元素.先不论从长度为n的序列中如何高效率的随机取出 ...
【算法】【树】已知先序中序序列求后序序列（详细解释）
题目描述如题所示,已知先序中序序列建树与求后序序列算法原理利用递归和分制的思想,找到当前树先序序列的根节点,然后找到对应中序序列的位置,然后根据根节点在中序序列中的位置来判断左右子树分别的位置, ...
求在整数序列A中连续权值最大的子序列的权值
题目描述给定整数序列A. 求在整数序列A中连续权值最大的子序列的权值. 输入描述: 第一行输入整数n.(1<=n<=1000) 第二行输入n整数a.(-1000<=a<=10 ...
有限长序列的z变换收敛域_求出以下序列的Z变换及收敛域
[单选题]已知N点有限长序列X(k)=DFT[x(n)],0≤n,k [判断题]目前我国的证券投资基金均为契约型基金 [简答题]求出以下序列的Z变换及收敛域 [u(n)-u(n-10)] [简答题]用 ...
matlab用卷积积分求单自由度阻尼系统响应（零极点模型）
%spring1.m-用卷积积分求单自由度阻尼系统的强迫振动(先求系统的零极点模型) m=1;c=4;K=100;dt=0.015; %dt为卷积积分时的采样间隔 w0=sqrt(K/m); %固 ...
matlab 卷积求图像梯度,图像的卷积(滤波)运算(一)——图像梯度
1. 卷积/滤波原理首先要明确的一点是图像的卷积/滤波运算,是针对原图像每一个像素进行处理,得到一个新的图像的过程.那么进行怎么样的运算呢?要知道图像能够被人所识别,是因为图像中每个像素并不完全是离 ...
【 MATLAB 】【 MATLAB 】DFT的性质讨论（三）序列的循环卷积及其 MATLAB 实现
这篇博文要将的是循环卷积,循环卷积和线性卷积还是有很大区别的,我们都知道,两个N点序列之间的线性卷积会得出一个更长的序列,不得不再一次要将区间限制在 0 <= n <= N -1.因此代替 ...
求助：MATLAB中实现卷积运算和理论分析中的卷积运算有什么区别？
MATLAB中实现卷积运算和理论分析中的卷积运算有什么区别. 欢迎使用Markdown编辑器你好! 这是你第一次使用 Markdown编辑器所展示的欢迎页.如果你想学习如何使用Markdown编辑 ...

matlab无限长序列卷积,怎样求未知长度序列的卷积

matlab无限长序列卷积,怎样求未知长度序列的卷积相关推荐

最新文章

热门文章