sum of a series

2 次查看（过去 30 天）

显示更早的评论

dav 2013-3-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/67716-sum-of-a-series

采纳的回答： Azzi Abdelmalek

在 MATLAB Online 中打开

Hi,

Can someone please help me to write the following code using a summation

for i = 1:1000
v(i) =y(i) - a0- a1*y(i-1)-a2*y(i-2);
end

Thanks.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Azzi Abdelmalek 2013-3-19

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/67716-sum-of-a-series#answer_79151

在 MATLAB Online 中打开

a0=1;
a1=2;
a2=1;
y=rand(1,1000) % your data
coeff=[-a2 -a1 1];
for ii = 3:1000
  v(ii) =sum(y(ii-2:ii).* coeff)-a0;
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

更多回答（1 个）

Youssef Khmou 2013-3-18

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/67716-sum-of-a-series#answer_79149

编辑：Youssef Khmou 2013-3-18

在 MATLAB Online 中打开

hi, what is Y[nT]? random variable?

N=1000;
a0=2.36;
a1=4.56;
a2=5.57;
y=randn(N,1);
v=zeros(N,1);
v(1:2)=y(1:2); % Optional, take is as initial condition 
 % Using Loops
 for n=3:N
          v(n)=y(n)-a0-(a1*y(n-1))-(a2*y(n-2));
 end
 figure, plot(y), hold on, plot(v,'r'), hold off
 % Vectorization
 Y1=zeros(size(y));
 Y1(2:end)=y(1:end-1);
 Y2=zeros(size(y));
 Y2(3:end)=y(1:end-2);
 V=y-a0-a1*Y1-a2*Y2;
 % Comparaison
 figure, plot(v), hold , plot(V,'r'), hold off