How can I curve fit this diffusion equation?

3 次查看（过去 30 天）

Joaquin 2014-7-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/144201-how-can-i-curve-fit-this-diffusion-equation

编辑： Star Strider 2014-7-31

I am plotting the following points

%t is the time points
t=[0 300 600 1800 3600 7200 10800 14400 18000 21600 25200 36000 43200 86400 108000 129600];
%C1 is the X-RAY attenuation
C1=[0 3.4118 5.015 8.3483 14.5637 23.4199 30.2665 35.5957 40.3528 43.7849 47.6186 54.3443 58.3553 69.2554 70.6798 72.0403];
To curve fit those points I am using the following function
clearfunction C1=Conc(t,D)
C1=0;
z=2; %Contact interface cartilage-bath
h=2; %Cartilage height
%D is the diffusion coefficient that we are looking for
%t corresponds to the different time points
for i=1:1000
C=(((-1)^i)/(2*i-1))*(cos((i-0.5)pi(-0.5)))*(exp(-((i-0.5)^2)pi^2)(D/4)*t);
C1=C1+C;
end
C1=1+(4/pi)*C1
end