Matlab 二阶差分图
我正在做一个项目,我需要计算数据系列的二阶差分图。数据向量是xn,我需要计算以下值来绘制图形:Matlab 二阶差分图,matlab,Matlab,我正在做一个项目,我需要计算数据系列的二阶差分图。数据向量是xn,我需要计算以下值来绘制图形: X(n) = x(n+1)- x(n)-----Equation(1) Y(n) = x(n+2) - x(n-1)-------Equation(2) 那么Xn对Yn的情节将是我的草料。 如何在MATLAB中计算第一个方程的xn+1和第二个方程的xn+2。我尝试了以下代码来计算xn+1的值: X = zeros(1,N); % Assume the first element is A1 X(1)
X(n) = x(n+1)- x(n)-----Equation(1)
Y(n) = x(n+2) - x(n-1)-------Equation(2)
那么Xn对Yn的情节将是我的草料。
如何在MATLAB中计算第一个方程的xn+1和第二个方程的xn+2。我尝试了以下代码来计算xn+1的值:
X = zeros(1,N);
% Assume the first element is A1
X(1) = A1;
for n = 1:(N-1)
if rem(n,2) == 0
% even case
An = A0;
else
% odd case
An = A1;
end
% step to increment X(n) to X(n+1)
X(n+1) = An*X(n)
end
但我得到的是与输出相同的xn数据。
我感谢你对这个问题的任何意见。
提前谢谢大家
x=rand(1E4,1); %example vector
X=x(2:end)-x(1:end-1); %x(n+1)-x(n)
Y=[nan;x(3:end)-x(1:end-2)]; %x(n+2)-x(n-1)
plot(X,Y,'.')
Y的计算如下:
Y1=nan无法计算x1-1,x0不存在
Y2=x3-x1
Y3=x4-x2
等等。您的阵列x和y是什么样子的?从你的算法来看,你的x似乎是[A1,A1*A1,A0*A1^2,A0*A1^3…]。这是正确的吗?亲爱的Anthony,“x”是一个TXT文件,由4096个ASCII码的EEG时间序列样本组成。我以xn的形式读取这个文件,需要计算xn+1和xn+2。希望这能给你一个清晰的想法。谢谢你的回答。在信号和系统透视图中,信号xn在时间标度上移动以获得xn+1。我不明白你为什么在Y=[NaN;x3:end-x1:end-2]中使用'NaN';请给出任何建议,移动信号能给我所需的结果吗?亲爱的Dattaprasad,我之所以选择nan是因为我无法计算x1-1,x0不存在。所以Y的第一个值无法计算。是的,Gelliant,我知道了。这就是我想要的结果。谢谢你的回答。