| 查看: 717 | 回复: 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) |
回复此楼» 猜你喜欢
请问有知道博士后面上专家评审什么时候结束的么？
已经有10人回复
-
夏至，要不要硕博联谊
已经有10人回复
-
今年什么时候会评啊
已经有9人回复
-
大佬们，打的核磁氢谱与chemdraw预测的有偏差
已经有7人回复
-
为何全世界都用JCR分区，咱们要弄个中科院分区？意义何在？
已经有4人回复
-
offer选择
已经有10人回复
-
论文被硕导卡
已经有6人回复
-
三颈瓶加热
已经有5人回复
-
请问三个字的中文名字，在文中到底怎么表示？
已经有7人回复
-
高校辞职，要求赔偿，这到底合不合理
已经有36人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 关于Matlab算符号积分
已经有6人回复
- matlab 曲线拟合的问题，急急急！！拜托了
已经有27人回复
- 【分享】Matlab经典教程——从入门到精通【已搜索无重复】
已经有794人回复
- 求助四重积分的matlab求解程序
已经有5人回复
- MATLAB解方程出现??? Subscripted assignment dimension mismatch.
已经有3人回复
- 【求助】MATLAB、Mathematica和Maple这三大数学软件的优缺点
已经有41人回复
- 【讨论】进行复杂公式赋值计算时请用Maple，慎用MathCAD
已经有43人回复
- 【求助】请问怎么求函数方程的函数表达式？
已经有11人回复
5