如何解决以下方程,利用Matlab

2视图(30天)
嗨,所有
当我们解决了这个方程\开始{对齐}
f (z) = (0.10 - -0.3) z ^ {1} + (0.2121 - 0.0008) z ^ {2} + (0.9 + 0.001) z ^ {3} = 0
通过使用Matlab,我们得到以下结果,
通过Matlab数值解,
1 / z_1 = (0.0000 + 0.0000)
1 / z_2 = (-0.4716 - 0.4706)
1 / z_3 = (0.2359 + 0.4717)
当我解这个方程byhand我们了;
1 / z_1 = 0.0000 + 0.0000
1 / z_2 = (-0.4929597603 - -0.4990725263)
1 / z_3 = (0.2561831187 + 0.5002232656)
我的段,
通过使用Matlab,我怎么能重复的解决方案(例如100次)检查每次数值解仍然接近用手或成为不同的解决方案,我有吗?
通过其他方式我可以使用任何循环重复利用Matlab解决方案100年时间?
我很感谢任何帮助
谢谢

接受的答案

卢卡铁
卢卡铁 2023年1月19日
试着用这个函数:
函数checksol (N)% N是它能解决方程的次数(不使用一个庞大的数字,10是足够的)
信谊z;
索尔= {};
diffcount = 0;
情商= (0.10 - -0.3)* z ^ 1 + (0.2121 - 0.0008) * z ^ 2 + (0.9 + 0.001) * z ^ 3 = = 0;
lastsol =解决(eq, z);
jj = 1: N
currsol =解决(eq, z);
如果currsol ~ = lastsol
diffcount = diffcount + 1;%增加计数器,每次发现一个不同的解决方案
结束
lastsol = currsol;
结束
如果diffcount ~ = 0
disp (发现不同的结果的)
其他的
disp (的结果总是相同的)
结束
结束
在我的机器上,它解决了持之以恒,我会仔细检查你的手工解决方案而不是matlab。
1评论
阿以莎默罕默德
阿以莎默罕默德 2023年1月20日
谢谢 卢卡铁
我想运行这个函数,但我没有发现任何resopnd Matlab。
我很感谢任何帮助

登录置评。

答案(1)

约翰D 'Errico
约翰D 'Errico 2023年1月20日
你的问题是:
f = @ (z) (0.10 - -0.3) * z。^ (1) + (0.2121 - 0.0008) * z。^ (2) + (0.9 + 0.001) * z。^ (3);
请不要在一个方程中使用花括号。这不是有效的MATLAB语法。,MATLAB使用*相乘。没有隐含乘以一个常数。MATLAB找到了解决方案:万博 尤文图斯
信谊Z
zsol =解决(f (Z) Z,“maxdegree”3)
zsol =
使用MATLAB解决你声称已经发现是假,因为z = = 0除以零。所以z = 1/0不是一个方程的解,即使你发现它用z = 1 / u,然后解决根使用u。对不起。只有两种解决方案。万博 尤文图斯
vpa (zsol)
ans =
逆这些数字是什么,所以我们可以把它比作声称,他们已经发现使用MATLAB ?
1. / vpa (zsol)
ans =
这是相同的MATLAB的回报。
我们可以验证这些确实是真正的解决你的问题。万博 尤文图斯
f (vpa (zsol))
ans =
所以数字的1 e-40。因为这是vpa的默认数量的数字,这只是浮点垃圾。
有什么有意义的原因这个解决方案执行一百万次,甚至100 ?不。这些解决方案,随着时间的推移万博 尤文图斯,他们不会改变。
不管怎样,你声称已经发现的解决方万博 尤文图斯案是用手真的解决问题你展示吗?
z2invclaimed = -0.4929597603 - -0.4990725263我;
f (1 / z2invclaimed)
我答= 0.0240 - 0.0130
所以实际上不是一个解决方案。这是接近。你使近似,删除结果只有几位有效数字,所以你自称是手工解决方案,你发现只是近似的解决方案。万博 尤文图斯有效,它们大约解决的问题写下来。万博 尤文图斯
再一次,多次重复的解决方案如您希望将无关紧要。结果不会改变。
2的评论
Torsten
Torsten 2023年1月22日
问题是为什么解的数值解精度比我们有手动吗?这是有意义的吗?
没有约翰D 'Errico试图与他最后回答这个句子吗?
你使近似,删除结果只有几位有效数字,所以你自称是手工解决方案,你发现只是近似的解决方案。万博 尤文图斯有效,它们大约解决的问题写下来。万博 尤文图斯
然后呢?你把数字从多项式系数的数值解?

登录置评。

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!