How can I use one matrix to assign values to a second matrix using a for loop?

Question

Harold Sandahl 2019-11-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/493920-how-can-i-use-one-matrix-to-assign-values-to-a-second-matrix-using-a-for-loop

评论： Harold Sandahl 2019-11-29

I need to create a diffusion simulation for class and I am stuck on the linear algebra. The first point needs to remain constant and the second point begins to be changed. I need the next point to be calculated through the equation Ci + phi*(Ci-1 - 2*Ci + Ci+1). It appears that I am doing something wrong with the loop itself. Any help is much appreciated. Thank you.

matrix = [1 1 1 1 1 1 1 1 1 1];
vector = [0 0 0 0 0 0 0 0 0 0];
phi = 0.5;
n=1;
for n=1:1:10
    if n = 1
        vector(n)=matrix(1)
    else
        vector(n)=matrix(n) + phi*(matrix(n-1)+matrix(n+1)-2*matrix(n));
    end
end

2 个评论
显示无隐藏无

KALYAN ACHARJYA 2019-11-29

编辑：KALYAN ACHARJYA 2019-11-29

No issue with Ci-1 (previous data), but how you get the next value of Ci+1?

Harold Sandahl 2019-11-29

C is a n by 1 matrix and the i represented the row of the matrix. So the i-1 means the previous row, the i means the current row, and the i+1 means the next row.

请先登录，再进行评论。

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

Answer 1

ME 2019-11-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/493920-how-can-i-use-one-matrix-to-assign-values-to-a-second-matrix-using-a-for-loop#answer_403939

Based on your definition there is nothing to diffuse. If the values is the same all along your domain then there won't be any change because you will always have vector(n) = matrix(n) as the rest of your expression will cancel out.

I would also have thought you'd get an error of index exceeds the number of array elements. This would be because you are trying to take Ci+1 as the point just outside of your arrays. To fix that issue you'll have to consider and implement some kind of boundary condition on the n=10 end.

I hope that helps! If you still have problems then you may need to explain exactly what issues you are having.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Harold Sandahl 2019-11-29

This was just an attempt to get the code to run before I implement actual values into it. The boundary conditions make sense, but as it stands I am getting Incorrect use of '=' operator. To assign a value to avariable, use '='. To compare values for equality, use '=='. for every equal sign in the for loop.

I have not gotten an index error, but I will define both boundaries and see if that fixes things.

请先登录，再进行评论。