Solution approaching wrong value

2 次查看（过去 30 天）

Abigail Grein 2018-12-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/435051-solution-approaching-wrong-value

评论： Torsten 2018-12-11

Hello, I have written the following code as an analytical solution to excess pore pressure in 1D soil consolidation:

z=(1:-delZ:0)';
tt=[0.01,0.1,0.2,0.5,1,2];
a=4/pi; pa=pi/2; pp=pa*pa;
Wan=zeros(21,1);
    for k=1:5
        jj(k)=(2*k-1);
        for ii=2:N+1
        Wan(ii)=Wan(ii)+(((-1)^(k-1))/jj(k))*cos(jj(k)*pa*z(ii))*exp(-jj(k)*jj(k)*pp*tt(1));
        end
    end      

My problem is it results in this graph for t=0.01, with the values approaching 0.8 when they should be approaching 1.0. Can anybody spot where I've gone wrong with this? Thank you!