应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 重金求Matlab计算一个简单的问题
北京石油化工学院2026年研究生招生接收调剂公告
61/1返回列表
查看: 1050  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

likeweid

金虫 (正式写手)

[求助] 重金求Matlab计算一个简单的问题

希望会的高手帮下忙，希望能提供一下计算结果及编好的matlab程序，感激不尽
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-09-22 20:46:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmc_1979

金虫 (小有名气)

【答案】应助回帖

likeweid(金币+200): 2011-09-25 18:55:52
clc;
x0=0.01;
x1=0.09;
[kk,kk1]=fun(x0);
while abs(kk)>0.000000001
    x1=x0-kk/kk1;
    f=x0;
    x0=x1;
    num=num+1;  
   [kk,kk1]=fun(x0);
end
x0
num
编写M文件
function [y,d]=fun(x)
   y=(5.7182*x^2-9.4364*x+5.577)*log((1-x)/x)+(1-2*x)*(3.7182*log(x)-6.799);
   d=(28591/2500*x-23591/2500)*log((1-x)/x)+(28591/5000*x^2-23591/2500*x+5577/1000)*(-1/x-(1-x)/x^2)/(1-x)*x-18591/2500*log(x)+6799/500+18591/5000*(1-2*x)/x;
关于多少次迭代的问题,主要看你的初值与精度要求.比如x0=0.01;精度0.000000001(函数值)得到根 0.02633273389612迭代214次
x0=0.4精度0.000000001得到根0.5迭代222次;
注意牛顿法的收敛性与初值是有关的,有可能的不到要求的解.
应该该题就两个解
赞一下 回复此楼
2楼2011-09-24 22:36:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likeweid

金虫 (正式写手)
引用回帖:
2楼: Originally posted by wmc_1979 at 2011-09-24 22:36:08:
clc;
x0=0.01;
x1=0.09;
[kk,kk1]=fun(x0);
while abs(kk)>0.000000001
    x1=x0-kk/kk1;
    f=x0;
    x0=x1;
    num=num+1;  
   [kk,kk1]=fun(x0);
end
x0
num
编写M文件
function [y,d] ...

大哥，你的导数是自己求的吗？感觉好像有点问题哦
赞一下 回复此楼
3楼2011-09-25 18:39:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likeweid

金虫 (正式写手)
引用回帖:
2楼: Originally posted by wmc_1979 at 2011-09-24 22:36:08:
clc;
x0=0.01;
x1=0.09;
[kk,kk1]=fun(x0);
while abs(kk)>0.000000001
    x1=x0-kk/kk1;
    f=x0;
    x0=x1;
    num=num+1;  
   [kk,kk1]=fun(x0);
end
x0
num
编写M文件
function [y,d] ...

我不是太懂，你里面的d是什么啊？导数还是那个迭代的函数？
赞一下 回复此楼
4楼2011-09-25 18:42:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likeweid

金虫 (正式写手)
引用回帖:
2楼: Originally posted by wmc_1979 at 2011-09-24 22:36:08:
clc;
x0=0.01;
x1=0.09;
[kk,kk1]=fun(x0);
while abs(kk)>0.000000001
    x1=x0-kk/kk1;
    f=x0;
    x0=x1;
    num=num+1;  
   [kk,kk1]=fun(x0);
end
x0
num
编写M文件
function [y,d] ...

能QQ交流一下吗？感激不尽，672468529
赞一下 回复此楼
5楼2011-09-25 18:47:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmc_1979

金虫 (小有名气)

【答案】应助回帖

引用回帖:
4楼: Originally posted by likeweid at 2011-09-25 18:42:17:
我不是太懂，你里面的d是什么啊？导数还是那个迭代的函数？

导数是matlab求的,我把结果拿过来,d就是导数表达式,对每一个x值返回一个导数值
赞一下 回复此楼
6楼2011-09-25 18:48:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 likeweid 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +9 一样YWY 2026-04-06 9/450 2026-04-06 17:17 by lijunpoly
[考研] 0854求调剂 +6 亨氏番茄沙司 2026-04-06 7/350 2026-04-06 16:38 by 蓝云思雨
[考研] 0703总分331求调剂 +10 ZY-05 2026-04-04 14/700 2026-04-06 11:51 by xhai2011
[考研] 一志愿武汉理工大学080200机械工程308分，求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
[考研] 085410人工智能 初试316分 求调剂 +7 残星拂曙 2026-03-31 7/350 2026-04-06 10:09 by 蓝云思雨
[考研] 320分人工智能调剂 +8 振—TZ 2026-04-03 8/400 2026-04-05 22:33 by 范式思维
[考研] 288求调剂，一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 413求调剂 +4 柯某某 2026-03-31 4/200 2026-04-04 22:18 by 学员6BFVa3
[考研] 0854求调剂 +4 assdll 2026-04-03 4/200 2026-04-04 22:17 by hemengdong
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 一志愿华北电力大学（北京），材料科学与工程学硕265，求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 一志愿沪9，求生物学调剂，326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
[考研] 0703求调剂 +6 zizimo 2026-03-31 6/300 2026-04-04 14:16 by 无际的草原
[考研] 278求调剂 +6 Yy7400 2026-04-03 6/300 2026-04-04 09:53 by zhangdingwa
[考研] 生物学硕341求调剂 +4 你笑起来像云朵 2026-04-03 4/200 2026-04-03 10:32 by macy2011
[考研] 312 化工或制药调剂 +8 小小墨123 2026-04-02 9/450 2026-04-03 09:12 by zhouxiaoyu
[考研] 初试301，代码085701环境工程，本硕一致，四六级已过，有二区一作，共发表5篇论文 +6 axibli 2026-04-01 6/300 2026-04-02 13:42 by Ecowxq666！
[考研] 一志愿北交大材料工程总分358 +8 cs0106 2026-04-01 9/450 2026-04-02 10:36 by 不吃魚的貓
[考研] 307分求调剂 +14 (o~o) 2026-03-31 15/750 2026-04-01 20:43 by longlotian
[考研] 313求调剂 +6 卖个关子吧 2026-03-31 6/300 2026-03-31 10:58 by Jaylen.
信息提示
关闭
请填处理意见
关闭