| 查看: 602 | 回复: 3 | |||
chuzhenhui金虫 (小有名气)
|
[交流]
【求助】一个方程，怎么解？？ 已有3人参与
|
方程如下： 已知X ̅ 和t，求Deff，其他都是常数。 n好像到8项以后值很小了，误差就可以忽略 想请教下各位，这个方程怎么解的？用循环语句+solve（）可以吗 非常感谢！ [ Last edited by chuzhenhui on 2010-5-10 at 11:05 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这样一个偏微分方程有没有解析解？
已经有4人回复
- 椭圆形偏微分方程标准化后该如何解？
已经有6人回复
- 一个偏微分方程组（包括边界条件）如何求解得到自振频率和振型？
已经有4人回复
- 【求助】如何解这样的方程组【有效期至2011年4月30日】
已经有9人回复
- 【求助】请教一个关于解方程的问题
已经有8人回复
- 【求助】求教 多元指数方程如何解 有软件否？
已经有6人回复
- 【求助】100金币求助高人一个程序解方程
已经有15人回复
- 【求助】一个常微分方程的通解
已经有6人回复
- 【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】
已经有15人回复
- 一个方程组的解法
已经有8人回复
2楼2010-05-10 13:04:39
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12258.2
- 散金: 200
- 帖子: 642
- 在线: 680.6小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
★ ★ ★
robert2020(金币+2):谢谢提供解决问题程序，欢迎常来仿真模拟版交流。 2010-05-14 22:37:45
zzuwangshilei(金币+1):帮助他人解决问题 2010-05-15 11:52:17
chuzhenhui(金币+6):非常感谢，程序代码我回去好好研究研究，努力学习！ 2010-05-16 12:52:03
robert2020(金币+2):谢谢提供解决问题程序，欢迎常来仿真模拟版交流。 2010-05-14 22:37:45
zzuwangshilei(金币+1):帮助他人解决问题 2010-05-15 11:52:17
chuzhenhui(金币+6):非常感谢，程序代码我回去好好研究研究，努力学习！ 2010-05-16 12:52:03
|
我觉得用2分法就可以了，我不知道方程已知系数取值。这里随便给出 作为例子，程序和结果如下： clc;clear;close all; Xp=0.5; % Xp是X ̅ Xe=1; Xo=8; L=1; t=1; Deff=linspace(-8,0); % Deff的区间可能与已知系数有关，原则是该区间内包含一个零点 F=(Xp-Xe)/(Xo-Xe); for n=0:10; F=F+8*exp([2*n+1]^2*pi^2*Deff*t/[4*L^2])/[pi*(2*n+1)^2]; end plot(Deff,F); % 作图法确定Deff两个端点值，使之F(Deff1)*F(Deff2)<0; Deff1=-8; Deff2=0; F1=(Xp-Xe)/(Xo-Xe); for n=0:10; F1=F1+8*exp([2*n+1]^2*pi^2*Deff1*t/[4*L^2])/[pi*(2*n+1)^2]; end F2=(Xp-Xe)/(Xo-Xe); for n=0:10; F2=F2+8*exp([2*n+1]^2*pi^2*Deff2*t/[4*L^2])/[pi*(2*n+1)^2]; end while abs(F1*F2)>1e-10; % 二分法 Deffm=[Deff1+Deff2]/2; Fm=(Xp-Xe)/(Xo-Xe); for n=0:10; Fm=Fm+8*exp([2*n+1]^2*pi^2*Deffm*t/[4*L^2])/[pi*(2*n+1)^2]; end if Fm*F2>=0; Deff2=Deffm; F2=Fm; else Deff1=Deffm; F1=Fm; end end Deffr=[Deff1+Deff2]/2; Fr=[F1+F2]/2; hold on;plot(Deffr,Fr,'ro'); % 画出解所在点 plot(xlim,[0,0],'r'); 运行上述程序可得到下图：  |
3楼2010-05-14 19:02:07
chuzhenhui
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 600.7
- 散金: 35
- 帖子: 129
- 在线: 10.2小时
- 虫号: 763256
- 注册: 2009-05-05
- 性别: GG
- 专业: 食品&包装
4楼2010-05-16 12:53:01
20