| 查看: 2530 | 回复: 15 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
各位朋友，如何用matlab编写这个方程，并求解？
|
|||||
|
各位朋友，最近遇到这么一个方程，不知道如何用matlab解决，方程如下： 公式中的未知数是t，其他的都是已知的。怎么用matlab编写这个方程，并求解？请各位朋友帮忙，谢谢！  方程 |
回复此楼» 收录本帖的淘帖专辑推荐
 学习知识 |
» 猜你喜欢
人气不行了
已经有10人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有6人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab求解方程中的参数
已经有21人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab求解非线性方程组
已经有16人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】如何用matlab求解这个微分方程？【已解决】
已经有3人回复
- 【资源】用matlab求解超越方程
已经有13人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★
xiegangmai: 金币+2, 专家考核, 辛苦了！ 2012-07-30 21:50:33
xiegangmai: 金币+2, 专家考核, 辛苦了！ 2012-07-30 21:50:33
|
或者这么写也可以 function yan global s phi lam alpha s=3.2; phi=1; lam=3; alpha=0.6; [X,FVAL,EXITFLAG,OUTPUT]=fsolve(@sumt,10) function y=sumt(t) global s phi lam alpha y=0; % for N=0:floor(s*t-phi) % y=y+exp(-lam*t)*(lam*t)^N/gamma(N+1); % end N=0:floor(s*t-phi); yN=exp(-lam*t)*(lam*t).^N./gamma(N+1); y=sum(yN); y=y-alpha; |
4楼2012-07-30 21:15:00
2楼2012-07-30 20:39:08
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
xiegangmai: 金币+3, 谢谢应助！ 2012-07-30 21:50:21
yan_fei: 金币+130, 感谢这位朋友的回复，金币全部给你！不过不知道为什么初值不同，结果差别很大，我想只改变alph的值求相应的t，随便一个初值都能求出准确的解，如知道望告知！感谢您的帮助~~~ 2012-07-30 22:39:56
感谢参与，应助指数 +1
xiegangmai: 金币+3, 谢谢应助！ 2012-07-30 21:50:21
yan_fei: 金币+130, 感谢这位朋友的回复，金币全部给你！不过不知道为什么初值不同，结果差别很大，我想只改变alph的值求相应的t，随便一个初值都能求出准确的解，如知道望告知！感谢您的帮助~~~ 2012-07-30 22:39:56
|
Equation solved. fsolve completed because the vector of function values is near zero as measured by the default value of the function tolerance, and the problem appears regular as measured by the gradient. X = 10.0881 FVAL = -9.6716e-010 EXITFLAG = 1 OUTPUT = iterations: 2 funcCount: 6 algorithm: 'trust-region dogleg' firstorderopt: 2.0551e-010 message: [1x695 char] |
3楼2012-07-30 21:09:48
|
首先感谢您的回复，我刚试了下您的程序，好像对初始值的依赖性很强，没次给不同的初始值求的结果都不一样，而且结果几乎和给的初始值差不多。我是想改变不同的alpha的值求出对应的t. 程序运行后的提示为： No solution found. fsolve stopped because the relative size of the current step is less than the default value of the step size tolerance squared, but the vector of function values is not near zero as measured by the default value of the function tolerance. fsolve stopped because the relative norm of the current step, 4.694166e-013, is less than max(options.TolX^2,eps) = 1.000000e-012. However, the sum of squared function values, r = 9.759502e-003, exceeds sqrt(options.TolFun) = 1.000000e-003. Optimization Metric Options relative norm(step) = 4.69e-013 max(TolX^2,eps) = 1e-012 (default) r = 9.76e-003 sqrt(TolFun) = 1.0e-003 (default) fsolve stopped because the relative norm of the current step, 4.694166e-013, is less than max(options.TolX^2,eps) = 1.000000e-012. However, the sum of squared function values, r = 9.759502e-003, exceeds sqrt(options.TolFun) = 1.000000e-003. Optimization Metric Options relative norm(step) = 4.69e-013 max(TolX^2,eps) = 1e-012 (default) r = 9.76e-003 sqrt(TolFun) = 1.0e-003 (default) |
5楼2012-07-30 21:53:23
20