matlab 组合结果有个题,想把exact solution里面循环的出来得结果组合到一起,最好形成一个类似于N*1的矩阵形式,以下是我的代码:N=6 % n cellsnsquare=25Ta=293,Tb=373a=0,b=1L=b-adeltax=L/NA=(-2-nsquare*deltax*delt
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/28 11:00:13
![matlab 组合结果有个题,想把exact solution里面循环的出来得结果组合到一起,最好形成一个类似于N*1的矩阵形式,以下是我的代码:N=6 % n cellsnsquare=25Ta=293,Tb=373a=0,b=1L=b-adeltax=L/NA=(-2-nsquare*deltax*delt](/uploads/image/z/4532885-53-5.jpg?t=matlab+%E7%BB%84%E5%90%88%E7%BB%93%E6%9E%9C%E6%9C%89%E4%B8%AA%E9%A2%98%2C%E6%83%B3%E6%8A%8Aexact+solution%E9%87%8C%E9%9D%A2%E5%BE%AA%E7%8E%AF%E7%9A%84%E5%87%BA%E6%9D%A5%E5%BE%97%E7%BB%93%E6%9E%9C%E7%BB%84%E5%90%88%E5%88%B0%E4%B8%80%E8%B5%B7%2C%E6%9C%80%E5%A5%BD%E5%BD%A2%E6%88%90%E4%B8%80%E4%B8%AA%E7%B1%BB%E4%BC%BC%E4%BA%8EN%EF%BC%8A1%E7%9A%84%E7%9F%A9%E9%98%B5%E5%BD%A2%E5%BC%8F%2C%E4%BB%A5%E4%B8%8B%E6%98%AF%E6%88%91%E7%9A%84%E4%BB%A3%E7%A0%81%EF%BC%9AN%3D6+%25+n+cellsnsquare%3D25Ta%3D293%2CTb%3D373a%3D0%2Cb%3D1L%3Db-adeltax%3DL%2FNA%3D%28-2-nsquare%2Adeltax%2Adelt)
matlab 组合结果有个题,想把exact solution里面循环的出来得结果组合到一起,最好形成一个类似于N*1的矩阵形式,以下是我的代码:N=6 % n cellsnsquare=25Ta=293,Tb=373a=0,b=1L=b-adeltax=L/NA=(-2-nsquare*deltax*delt
matlab 组合结果
有个题,想把exact solution里面循环的出来得结果组合到一起,最好形成一个类似于N*1的矩阵形式,以下是我的代码:
N=6 % n cells
nsquare=25
Ta=293,Tb=373
a=0,b=1
L=b-a
deltax=L/N
A=(-2-nsquare*deltax*deltax)*eye(N);
for i=1:N-1
A(i,i+1)=1;
A(i+1,i)=1;
end
A(1,1)=-3-nsquare*deltax*deltax;
A(N,N)=-1-nsquare*deltax*deltax
B=-Ta*nsquare*deltax*deltax*ones(N,1);
B(1,1)=-Ta*nsquare*deltax*deltax-2*Tb
T=A\B
% exact solution
for x=deltax/2:deltax:L
t=sqrt(nsquare)*(L-x)
cosht=(exp(t)+exp(-t))/2
Texact=cosht*(Tb-Ta)/cosh(sqrt(nsquare)*L)+Ta
end
%
matlab 组合结果有个题,想把exact solution里面循环的出来得结果组合到一起,最好形成一个类似于N*1的矩阵形式,以下是我的代码:N=6 % n cellsnsquare=25Ta=293,Tb=373a=0,b=1L=b-adeltax=L/NA=(-2-nsquare*deltax*delt
没必要用循环,直接用向量运算就可以了:
x=deltax/2:deltax:L;t=sqrt(nsquare)*(L-x);
cosht=(exp(t)+exp(-t))/2;
Texact=cosht*(Tb-Ta)/cosh(sqrt(nsquare)*L)+Ta
得到的结果与用循环一致:
Texact =345.7424 315.9320 302.9899 297.3960 295.0358 294.1730