Help Plotting Diffusion Equation in MatLab

Question

Anthony Koning 2021-9-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1459509-help-plotting-diffusion-equation-in-matlab

回答： Star Strider 2021-9-23

I am currently working on a diffusion distribution project for an engineering class. The class requires me to use MatLab for many of the problems, a program that I have minimal experience with. Currently, I am trying to figure out how to plot Fick's second law (equation C(x,t)=N*e^(-x^2/(4Dt))/sqrt(4piDt), where D and N are given constants.). I was wondering if someone would be able to teach me how to plot this function in MatLab or provide a template script to use. The final result should look similar like this:

I would appreciate any help or hints that someone could provide. Thank you very much, and have a nice day.

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2021-9-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1459509-help-plotting-diffusion-equation-in-matlab#answer_793769

在 MATLAB Online 中打开

One option is to stack them in the same axes —

x = linspace(-200, 200);

v = logspace(-2, 2, 10);

y = exp(-(x./v(:)).^2);

y = 10×100

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0005 0.0009 0.0017 0.0029 0.0049 0.0183 0.0215 0.0251 0.0293 0.0341 0.0394 0.0455 0.0524 0.0601 0.0687 0.0783 0.0889 0.1007 0.1136 0.1278 0.1433 0.1601 0.1783 0.1979 0.2190 0.2416 0.2655 0.2910 0.3178 0.3459 0.3753 0.4059 0.4376 0.4702 0.5036

figure

plot(x, y)

grid

figure

plot(x, y+log(v(:)))

grid

Other options are subplot or stackedplot.

.