| 查看: 431 | 回复: 0 | ||
小木虫zb木虫 (正式写手)
|
[求助]
matlab程序的问题
|
|
小弟用matlab计算一个偏微分方程，程序运行不出来，求大神相助，程序如下 function DBT_MNA%DBT用1表示,MNA用2表示 clear all;clc A=importdata('DBT+MNA-RS-03.txt');%数据源文件中横坐标为g-MDF/g-AC t2=A(:,1)/7.03;y2=A(:,2);%横坐标由g-MDF/g-AC化为h-1 q10=zeros(1,20); c10=zeros(1,20); q20=zeros(1,20); c20=zeros(1,20); y0=[c10 q10 c20 q20]; c0=7.03; options=odeset('relTol',1e-6); tspan=[0:0.02:t2(end)]; [t,y]=ode15s(@fangcheng,tspan,y0,options); plot(t*7.03,y(:,40)/c0,'ro',t*7.03,y(:,80)/c0,'bo','LineWidth',3);hold on plot(A(:,1),A(:,2),A(:,1),A(:,3),'LineWidth',3) xlabel('Amount of treated MDF(g-MDF/g-AC)'),ylabel('c/c0'),axis([-0.05 120 -0.05 1.05]), legend('拟合DBT','拟合MNA','实验曲线DBT','实验曲线MNA','Location','best') ,title('DBT+MNA-RS') grid,hold off %--------------------------------------------------------- function dydt=fangcheng(t,y) kf=200; a=0.363;%空隙率 u=0.361032/a;%床层颗粒间隙中流体流速 c0=7.03; L=8.31e-2; N=20; dz=L/N; p=436.8/(1-0.363);%颗粒表观密度对曲线影响较大从434.2该为436.8 q1=y(1:N);c1=y(N+1:2*N); q2=y(2*N+1:3*N);c2=y(3*N+1:4*N); [ce1,ce2]=equilibrium(q1,q2); %-----------定义dq/dt----------------------------- for j=1:N dq1dt(j)=kf/p*(c1(j)-ce1(j)); dq2dt(j)=kf/p*(c2(j)-ce2(j)); end %---------定义dC/dt----------------------------------- dc1dt(1)=-u*(c1(1)-c0)/dz-(1-a)/a*p*dq1dt(1); dc2dt(1)=-u*(c2(1)-c0)/dz-(1-a)/a*p*dq2dt(1); for j=2:N dc1dt(j)=-u*(c1(j)-c1(j-1))/dz-(1-a)/a*p*dq1dt(j); dc2dt(j)=-u*(c2(j)-c2(j-1))/dz-(1-a)/a*p*dq2dt(j); end dydt=[dq1dt dc1dt dq2dt dc2dt]'; function [ce1,ce2]=equilibrium(q1,q2) [ce1,ce2]=solve('q1=0.45967ce1^0.3839/(ce1^0.0457+0.3*ce2^0.0003206',... 'q2-0.17356*ce2^(-1.55201)/(ce2^(-2)+0.1*ce1^(-1))');[ Last edited by 小木虫zb on 2013-4-26 at 09:30 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : DBT+MNA-RS-03.txt
2013-04-26 09:30:48, 990 bytes
» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
面上意见出来了
已经有18人回复
-
找导师
已经有7人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
找到一些相关的精华帖子，希望有用哦~
求助MATLAB做三维图问题
已经有4人回复
- matlab语句问题
已经有5人回复
- matlab中遇到的画图小问题，程序没问题。标题要长@@@@@@@@@@@@@@@@@@@@@@
已经有8人回复
- 求助一个matlab关于分类讨论的问题
已经有9人回复
- matlab编程中遇到的问题
已经有4人回复
- Matlab中save问题
已经有5人回复
- 问下这个问题怎么解决.....真心啊
已经有4人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- 请教关于matlab矩阵转换问题
已经有6人回复
- 关于Matlab作图的问题
已经有14人回复
- matlab的函数调用问题
已经有3人回复
- 问个关于matlab里find函数的问题！
已经有4人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求助使用MATLAB计算绿色叶片的覆盖度的问题
已经有24人回复
- matlab非线性参数拟合问题
已经有7人回复
- MATLAB 画图小问题
已经有11人回复
- 同样的MATLAB程序，为何在别人的电脑上运行不了
已经有21人回复
- 关于matlab 2010b 的一个小问题
已经有7人回复
- matlab做函数图象的问题
已经有6人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
