| 查看: 714 | 回复: 3 | |||
[交流]
【求助】MATLAB的符号运算已有3人参与
|
|
谁能帮我把下面的程序中 yy变成容易看的形式： syms er y=[0.162936345 0.139659921 0.0658347 0.00551692 0.00062984 9.36509E-06 7.83179E-06 7.95103E-06 7.41275E-06]; %实验测定值 h=0.03; %床层高度 ri=[0 0.003 0.006 0.009 0.012 0.015 0.02 0.025 0.03]; %径向上测验点的距离 ai=[3.832 7.016 10.173 13.324 16.471 19.616 22.760 25.904 29.047 32.1897]; % 1阶bessel函数的正根 r0=0.036; %床层半径 uz=0.001474; %床层内液速 zc=1/180; %示踪剂与水流量之比 x=ri/r0; %er=2.22856425*1.0e-8; %径向扩散系数Er，给定初值 b1=besselj(0,ai*x(1)); %求b b2=besselj(0,ai*x(2)); b3=besselj(0,ai*x(3)); b4=besselj(0,ai*x(4)); b5=besselj(0,ai*x(5)); b6=besselj(0,ai*x(6)); b7=besselj(0,ai*x(7)); b8=besselj(0,ai*x(8)); b9=besselj(0,ai*x(9)); j0=besselj(0,ai); % 0阶bessel函数 w=zc./j0.^2; e=(uz*r0^2-(uz^2*r0^4+4*er^2*ai.^2).^0.5)*h/2/er; g1=w.*b1.*exp(e); g2=w.*b2.*exp(e); g3=w.*b3.*exp(e); g4=w.*b4.*exp(e); g5=w.*b5.*exp(e); g6=w.*b6.*exp(e); g7=w.*b7.*exp(e); g8=w.*b8.*exp(e); g9=w.*b9.*exp(e); y1=zc+sum(g1); %求y y2=zc+sum(g2); y3=zc+sum(g3); y4=zc+sum(g4); y5=zc+sum(g5); y6=zc+sum(g6); y7=zc+sum(g7); y8=zc+sum(g8); y9=zc+sum(g9); yn=[y1 y2 y3 y4 y5 y6 y7 y8 y9] yn1=diff(yn); yy=sum((yn-y).*yn1) |
回复此楼» 猜你喜欢
工材口青年基金大概什么样能上会？
已经有9人回复
-
怎么带研究生？
已经有4人回复
-
这个博士要读吗
已经有29人回复
-
求助
已经有4人回复
-
累并快乐着
已经有21人回复
-
豫北虫友互识
已经有12人回复
-
试剂价格
已经有6人回复
-
34岁读博士晚吗
已经有41人回复
-
2024国社科通讯评审
已经有45人回复
-
刚刚收到科研之友邮件
已经有25人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 关于Matlab算符号积分
已经有6人回复
- matlab 曲线拟合的问题，急急急！！拜托了
已经有27人回复
- 【分享】Matlab经典教程——从入门到精通【已搜索无重复】
已经有794人回复
- 求助四重积分的matlab求解程序
已经有5人回复
- MATLAB解方程出现??? Subscripted assignment dimension mismatch.
已经有3人回复
- 【求助】MATLAB、Mathematica和Maple这三大数学软件的优缺点
已经有41人回复
- 【讨论】进行复杂公式赋值计算时请用Maple，慎用MathCAD
已经有43人回复
- 【求助】请问怎么求函数方程的函数表达式？
已经有11人回复
2楼2010-05-24 18:47:01
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
3楼2010-05-24 18:58:28
yznian
至尊木虫 (著名写手)
- 应助: 25 (小学生)
- 金币: 10304.6
- 散金: 69
- 红花: 39
- 帖子: 1653
- 在线: 417.4小时
- 虫号: 959310
- 注册: 2010-03-01
- 性别: GG
- 专业: 零件加工制造
4楼2010-05-25 11:41:01
10