How to apply an equation to each element in the array?

k = 10^-7 % filtration of soil
M = 19000
q = 80 % external force Q acting on soil
DeltaZ = 0.1 % depth of soil subsidization
GammaW = 10  % water weight [Kn/m^3]
Cv = k*M/GammaW % consolidation factor
DeltaT = DeltaZ^2/(2*Cv) % max iteration step
alpha = Cv*DeltaT/DeltaZ^2 
% NUMERICAL PARAMETERS
u = q*ones(10:10)  % matrix of u nodes
%symbolic variables
% SOLUTION
% discrete form of the governing equation:
soilSub = conv( u, [0 1 0]+alpha*[1 -2 1] );

-----------------------------------------------------

You know I'm pretty new to MATLAB using it for 6 weeks or so, I probably do not have enough knowledge for advanced models. I do not expect folks here will do my homework for me :) but I'm eager to discuss methods in which you can: 1) build a model of rock layer 2) implement finite difference equation to that model so it goes through it and computes external pressure.

Matt J 2021-11-28

在 MATLAB Online 中打开

Use conv2 instead

soilSub = conv2( u, [0 1 0]+alpha*[1 -2 1] )

KarolN 2021-11-28

IT WORKS! Thank you very much for your valuable advice!

请先登录，再进行评论。

Answer 2

Image Analyst 2021-11-28

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1597244-how-to-apply-an-equation-to-each-element-in-the-array#answer_842164

在 MATLAB Online 中打开

Try conv2():

k = 10^-7 % filtration of soil
M = 19000
q = 80 % external force Q acting on soil
DeltaZ = 0.1 % depth of soil subsidization
GammaW = 10  % water weight [Kn/m^3]
Cv = k*M/GammaW % consolidation factor
DeltaT = DeltaZ^2/(2*Cv) % max iteration step
alpha = Cv*DeltaT/DeltaZ^2 
% NUMERICAL PARAMETERS
u = q*ones(10:10)  % matrix of u nodes
% SOLUTION
kernel = [0, 1, 0] + alpha * [1, -2, 1]
% Discrete form of the governing equation:
soilSub = conv2(u, kernel) % Use conv2(), not conv()

4 个评论
显示 2更早的评论隐藏 2更早的评论

Image Analyst 2021-11-28

No problem. I'm sure using conv() instead of conv2() was just an oversight by Matt. I know he knows better - he's a bright guy. You can vote for multiple answers though you can only Accept one. You can even vote for one you Accepted to award double the reputation points to that one.

KarolN 2021-11-28

And so additionally I learned about that vote feat :) Voted for both of you naturally. Thanks guys

请先登录，再进行评论。

How to apply an equation to each element in the array?

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

采纳的回答

5 个评论
显示 3更早的评论隐藏 3更早的评论

更多回答（1 个）

4 个评论
显示 2更早的评论隐藏 2更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

How to apply an equation to each element in the array?

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

采纳的回答

5 个评论 显示 3更早的评论隐藏 3更早的评论

更多回答（1 个）

4 个评论 显示 2更早的评论隐藏 2更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

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

5 个评论
显示 3更早的评论隐藏 3更早的评论

4 个评论
显示 2更早的评论隐藏 2更早的评论