Matlab或者EXCEL高手进。。学术问题。。

hpPlayer

回家作业有道编程题，可以使用任何的程序编辑。。但是我之前从来没有编过程。。。这不是故意难为我么。。我电脑上只有EXCEL和MATLAB，大家想想办法吧。。。。废话不多说，上题：
一条河流的初始浓度公式已经知道：
C=1/((2*3.14*0.5^2)^0.5)*EXP((-(x-50)^2)/(2*0.5^2))   ‘ X为与原点的距离，范围为0:到100
t+△t时间后，浓度的分布公式已经知道：
C(i)(t+△t)=△t*D/((△x)^2)*(C(i+1)(t)+C(i-1)(t)-2C(i)(t))+C(i)(t)
i=任一点，i+1, 任一点后面一个单位的点，i-1, 任一点前面一个单位的点.
t, 任意时间；△t， 时间的差值，可以任意假设(但是要小于0.08,。。。经验公式所得)，我在做题时，为了方便计算，假设△t=0.05；D=0.5；
boundary conditions are zero on both sides.
题目让求10,50,100,200,500s 后的浓度分布公式，用程序形式表达。
我的解法是，算出各个点的初始浓度C，
然后代入到公式2，就可以求出某一点△t=0.5时间后的浓度，
同理求得△t=0.5时间后该点前后点的浓度
如此反复代入，最后就可以求出题目要求时间内的浓度公式。
很明显，这是个循环语句。。。很高深。。我不会编程。。。
哪位高手能帮我。。。谢谢大家了。。悬赏分神马的都是浮云，要多少给多少。。。只求做出这道题。。。

MartiniBlanco

1000.........

tensor

程序已经调试过，可以用。程序必须是英文字母，注释%后面可以写汉字。
matlab是矩阵运算，所以变量乘、除、平方前面加字符“.”，例如*应该写成.*
回复不能加.m的附件，我压缩成了。rar的文件，在附件中。
例子中，取时间t=10，你也可以取别的值（50,100,200,500s），我有事，所以简单的这样编了。
x分点我取了100个，不应该取得太大，否则间距dx太小，第二个方程dx的平方在分母上，会得到很大的值，应该避免。

t=10; % 输入时间t，这里取t=10，t也可以取50,100,200,500s
x=linspace(0,100,99);   % x 范围[0,100],平均分为100个点,改变第三个输入参数99，可以改变分的点数，x如果分得太大，dx太小，浓度计算有问题，自己分析
C=1/((2*3.14*0.5^2)^0.5)*exp((-(x-50).^2)/(2*0.5^2)); %计算x对应的C值，共1000个
C(1)=0;c(1000)=0;  %边界为零
%计算时间t的浓度
dt=0.05;
dx=x(2)-x(1)  %计算x的间距dx
D=0.5;
n=t/dt;
for i=1:n
    for j=2:length(x)-1    %边界为零，只改变除边界的其他值
        C(j)=dt*D/(dx^2)*(C(j+1)+C(j-1)-2*C(j))+C(j);
    end
end
% 输出t时间的浓度值
C,

anomaly

楼上真有干劲， 不愧是学 tensor 的

xjzf

楼上是高手

tensor

把程序中的C(1000)=0，改成C(100)=0，我没注意看，只有100个分点！下面附件是修改过的

代表太阳惩罚你

niuB...

菊菊思密达

我草，果断进来膜拜程序帝！

菊菊思密达

看得我等彩笔一愣一愣的，
犹如当年第一眼看到包哥第一视角时的震惊

baggiozerg

tensor

tensor同学你那个ij二重循环会不会有点小问题啊，比如你那个C(j)=dt*D/(dx^2)*(C(j+1)+C(j-1)-2*C(j))+C(j)中间，C(j+1)还是t时刻的不是t+dt的。
建议用二维的C(i,j)=C(t+idt,Xj)来表示。
lofi 发表于 2010-11-3 10:16

我当时有急事，没仔细看，谢谢。是有问题，我改改。

tensor

改正后的程序
%程序开始。要求输入时间 t 和 x 的分的点数 N 。
t=input('请输入时间 t ：t = ');
N=input('请输入[0 x]区间均分的段数 N ：N = ');
x=linspace(0,100,N);    % x 范围[0,100],平均分为N个点。x如果分得太多，dx太小，浓度计算有问题，自己分析
C=1/((2*3.14*0.5^2)^0.5)*exp((-(x-50).^2)/(2*0.5^2));   %计算x对应的C的初始值
C(1)=0;  c(length(x))=0;  %边界为零
%计算时间t的浓度
dt=0.05;       %时间变化dt
dx=x(2)-x(1);  %计算x的间距dx
D=0.5;
n=t/dt;       %时间的段数
for i=1:n
    temp=C;    %每次计算前，保存C的值，否则影响计算
    for j=2:length(x)-1    %边界为零，只改变除边界的其他值
        C(j)=dt*D/(dx^2)*(temp(j+1)+temp(j-1)-2*temp(j))+temp(j);
    end
end
% 输出t时间的浓度值
fprintf('时间 t = %d 对应的浓度为：',t),
C,

kingcobra

进来观摩编程帝，真威武
将来有问题还是得来8DA

arcturus

stone2618

万能8达啊，看水晶

BeRush

哪个编程语言啊，%竟然是注释用的。。