24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

我要飞

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 1515.3
散金: 144
红花: 5
帖子: 458
在线: 591.2小时
虫号: 671362
注册: 2008-12-09
专业: 计算数学与科学工程计算

[求助] 求助matlab里面画出x^3-1=0的Newton迭代的吸引域的程序

求助matlab里面画出x^3-1=0的Newton迭代的吸引域的程序

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有71人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

如何用matlab将这个方程代表的图画出来，最好能把画图的代码贴出来，谢谢啊！已经有11人回复
求用牛顿迭代法解非线性方程组的matlab程序已经有4人回复
请教matlab如何增加迭代次数已经有11人回复
求用牛顿迭代法求非线性方程组已经有4人回复
matlab数值计算时迭代项分母为零的时候怎么办？已经有15人回复
matlab菜鸟问3个问题~~~~~~~~~~~ 已经有2人回复
求助matlab画出螺旋线的程序已经有6人回复
【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】已经有15人回复
【求助】求一个用牛顿迭代法解非线性方程组的程序已经有10人回复

1楼 2011-10-15 21:34:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yf_2011

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 230.5
散金: 26
帖子: 48
在线: 19.8小时
虫号: 1453478
注册: 2011-10-21
专业: 电力电子学

你好，matlab初学者，要画吸引域，但是不知道怎么画啊，程序保存为m文件后运行没有图像出来啊，谢谢了

赞一下

回复此楼

6楼2012-06-07 21:43:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 6 个回答

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

CODE:

function AA

p=f(-1/2-sqrt(3)/2*i , -1/2+sqrt(3)/2*i, 1)

function p=f(A,B,C)

%>>p=f(-1/2-sqrt(3)/2*i , -1/2+sqrt(3)/2*i, 1)

h=0.01;

for a=-1:h:1

 for b=-1:h:1

 x0=a+b*i;

 y=newton(x0);

 if abs(y-A)<1.0e-6

 plot(a,b,'r.');

 hold on

 elseif abs(y-B)<1.0e-6

 plot(a,b,'g.');

 hold on

 elseif abs(y-C)<1.0e-6

 plot(a,b,'y.');

 hold on

 end 

 end

end

p=0;

function y=newton(z)

if (z==0)

 y=0;

 return;

end

for i=1:1:2000

 y=z-(z^3-1)/(3*z^2);

 if(abs(y-z)<1.0e-7)

 break;

 end

 z=y;

end