| 查看: 1739 | 回复: 0 | ||
[求助]
matlab求解一个非齐次偏微分方程
|
|
方程的具体形式在文件附件中（3）中，用C-N方法解的话我编写的程序为 p=2.1;%密度 c=1.17;%比热容 ke=0.0063;%热导率 g=0.93;%各向异性系数 u1=3;%吸收系数 u2=260;%散射系数 dx=0.015; x=0:dx:0.15; J1=length(x); dt=1e-4; N=3001; t=(0:N-1)*dt; u3=sqrt(3*u1*(u1+(1-g)*u2)); l=ke*dt/(p*c*(dx^2)); d=dt/(p*c); I0=180; S=zeros(N,J1); I=zeros(N,J1); for i=2:J1-1 for m=1:9 I(m,1)=I0; I(m,i)=I0; I(m,J1)=I0; S(m,1)=u1*I(1,1); S(m,i)=u1*I(1,i)*exp(-u3*(i-1)*dx); S(m,J1)=u1*I(1,11)*exp(-u3*10*dx); S(10:N,  =0;end end u=zeros(N,J1); ubc_l=S(:,1)/(p*c); ubc_r=S(:,11)/(p*c); u_initial=310*ones(1,J1); u(1, =u_initial;a=ones(J1,1)*(l+1); b=-l/2*ones(J1-1,1); c=-l/2*ones(J1-1,1); f=zeros(J1,1); for n1=2:N temp1=l/2*ubc_l(n1-1)+(1-l)*u(n1-1,1)+l/2*u(n1-1,2); f(1)=temp1+l/2*ubc_l(n1)+1000*d*S(n1-1,1); f(2:J1-1)=l/2*u(n1-1,1:J1-2)+(1-l)*u(n1-1,2:J1-1)+l/2*u(n1-1,3:J1)+1000*d*S(n1-1,2:J1-1); f(J1)=l/2*u(n1-1,J1-1)+(1-l)*u(n1-1,J1)+l/2*ubc_r(n1-1)+l/2*ubc_r(n1)+1000*d*S(n1-1,J1); u(n1, =my_chasing(a,b,c,f);end plot(t,u(:,1)-273) 其中S相当于Q,u相当于T,为什么我画出的图形衰减太快，my_chasing是追赶法 另外对于这类方程我也不清楚它的边界调节如何确定。 请专家们帮帮我万分感谢！！！！！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 新建WPS文字文档(3).wps
2012-07-02 17:42:22, 12.5 K
» 猜你喜欢
304求调剂
已经有4人回复
-
326求调剂
已经有4人回复
-
一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项
已经有4人回复
-
环境工程调剂
已经有5人回复
-
本人考085602 化学工程 专硕
已经有14人回复
-
一志愿211 0703方向310分求调剂
已经有3人回复
-
0703化学336分求调剂
已经有3人回复
-
318求调剂
已经有3人回复
-
化学工程321分求调剂
已经有9人回复
-
0703化学调剂
已经有9人回复
找到一些相关的精华帖子，希望有用哦~
这样一个偏微分方程有没有解析解？
已经有4人回复
- Matlab求解偏微分方程组
已经有22人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 【求助】matlab求解一个简单的积分（急！！！）
已经有7人回复
- 【求助】用pdepe求解偏微分方程
已经有10人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】如何求解下列偏微分方程
已经有13人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】用MATLAB求解偏微分方程的数值解 悬赏20
已经有8人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
- 【分享】PDECOL-经典数值求解非线性偏微分方程组源程序
已经有142人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
