我怎么能解决下面的问题?

3视图(30天)

显示旧的评论

PULAK咕 2020年12月26日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/702947-how-can-i-solve-following-problems

评论道: 沃尔特·罗伯森 2020年12月28日

答:接受沃尔特·罗伯森

的代码是:

                          %的程序代码在MATLAB由Pulak找到根
                         
                          clc,清晰所有
                         
                          信谊x
                         
                          一个=输入(的输入变量x的函数的形式:“);
                         
                          呃= 10;
                         
                          disp (“你在你的数学有初始近似值吗?”)
                         
                          b =输入(如果是的,请按1,如果没有按2:“);
                         
                          如果(b ~ = 1)
                         
                          x (1) = 1;
                         
                          其他的
                         
                          x(1) =输入(“输入初始猜测:”);
                         
                          结束
                         
                          n =输入(进入小数位的);
                         
                          托尔= 1 / (10 ^ (n - 1))
                         
                          f =内联(a)
                         
                          dif = diff(信谊(a));
                         
                          d =内联(dif);
                         
                          为k = 2:1000
                         
                          x (k) = x (k - 1) - ((f (x (k - 1)) / d (x (k - 1))));
                         
                          呃= abs ((x (k) - x (k - 1)) / x (k));
                         
                          如果犯错<托尔
                         
                          打破
                         
                          结束
                         
                          结束
                         
                          k
                         
                          流(“% * f '。n x (k))

但在运行这段代码之后,我明白了:

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

输入变量x的函数形式:x ^ 2 + (4 * sin (x))

你有初始近似值你的数学吗?

如果是的,请按1,如果没有按2:2

输入十进制place4

托尔=

1.0000 e 03

f =

内联函数:

f (x) = sin (x) * 4.0 + x ^ 2

转换从符号逻辑是不可能的。

错误final_newton_raphson(22)行

如果犯错<托尔

我怎样才能解决这个?&我也可以不给输入复杂的值x(1)我也不能得到复杂根先生。我能解决这个问题吗?

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

接受的答案

沃尔特·罗伯森 2020年12月26日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/702947-how-can-i-solve-following-problems answer_584752

                              如果(b ~ = 1)
                             
                              x (1) = 1;
                             
                              其他的
                             
                              x(1) =输入(“输入初始猜测:”);
                             
                              结束

改变,

                              如果(b ~ = 1)
                             
                              x0 = 1;
                             
                              其他的
                             
                              x0 =输入(“输入初始猜测:”);
                             
                              结束

之前,

为k = 2:1000

插入

x = x0;

我也可以不给输入复杂的值x (1)

只要输入提示

> > x0 =输入(输入初始猜测:)

我输入初始猜测:3 + 5

x0 =

3 + 5

11日评论
显示十年长的评论隐藏10年长的评论

沃尔特·罗伯森 2020年12月27日

%的程序代码在MATLAB由Pulak找到根

信谊x

disp (的输入变量x的函数的形式:“);

输入变量x的函数形式:

% {

=输入();

= x ^ 4 + x ^ 3 + 5 * x ^ 2 + 4 * x + 4%的演示

一个=

呃= 10;

disp (“你在你的数学有初始近似值吗?”)

你有初始近似值你的数学吗?

disp (如果是的,请按1,如果没有按2:“)

如果是的,请按1,如果没有按2:

% {

b =输入();

b = 1%的演示

b = 1

如果(b ~ = 1)

x0 = 1;

其他的

disp (“输入初始猜测:”)

% {

x0 =输入();

x0 = 1我%的演示

结束

输入初始猜测:

我x0 = 0.0000 + 1.0000

disp (进入小数位的)

进入小数位

% {

n =输入();

n = 5%的演示

n = 5

托尔= 1 / (10 ^ (n - 1))

托尔= 1.0000 e-04

f =内联(a)

=内联函数:f (x) = x ^ 2 * 4.0 + x。* 5.0 + x ^ 3 + x。^ 4 + 4.0

dif = diff(信谊(a));

d =内联(dif);

x = x0;

为k = 2:1000

x (k) = x (k - 1) - ((f (x (k - 1)) / d (x (k - 1))));

呃= abs ((x (k) - x (k - 1)) / x (k));

如果犯错<托尔

打破

结束

k = 5

流(' %。* f % + * fi。\ n ',n,真正的(x (k)) n,图像放大(x (k)));

-0.50000 + 0.86603我

PULAK咕 2020年12月27日

当我们运行以下代码:

                                  %在MATLAB程序代码的牛顿迭代方法
                                 
                                  clc,清晰所有
                                 
                                  流(“* * * * * * * * <强>代码建立了由Pulak茶室(181201)< / >强* * * * * * * * \ n”)
                                 
                                  流(“* * * * * * * * <强>代码使用Newton-Rapshon虚根的方法< / >强* * * * * * * * \ n”)
                                 
                                  信谊x
                                 
                                  一个=输入(的输入变量x的函数的形式:“);
                                 
                                  呃= 10;
                                 
                                  disp (“你在你的数学有初始近似值吗?”)
                                 
                                  b =输入(如果是的,请按1,如果没有按2:“);
                                 
                                  如果(b ~ = 1)
                                 
                                  x0 = 1;
                                 
                                  其他的
                                 
                                  x0 =输入(“输入初始近似值:”);
                                 
                                  结束
                                 
                                  n =输入(“输入小数地方:”);
                                 
                                  托尔= 1 / (10 ^ (n - 1));
                                 
                                  disp (”“你的帐号多项式方程是:)
                                 
                                  一个
                                 
                                  f =内联(a);
                                 
                                  dif = diff(信谊(a));
                                 
                                  d =内联(dif);
                                 
                                  x = x0;
                                 
                                  为k = 2:1000
                                 
                                  x (k) = x (k - 1) - ((f (x (k - 1)) / d (x (k - 1))));
                                 
                                  呃= abs ((x (k) - x (k - 1)) / x (k));
                                 
                                  如果犯错<托尔
                                 
                                  打破
                                 
                                  结束
                                 
                                  结束
                                 
                                  流(的<强>使用Bewton-Raphson方法,我们得到了< / >强”)
                                 
                                  流(' \ n \ t \ t \ \ t <强>根是:%。* f + %。* fi \ n < / >强”,n,真正的(x (k)) n,图像放大(x (k)));
                                 
                                  流(“* * * * * * * * <强>代码建立了由Pulak茶室(181201)< / >强* * * * * * * * \ n”)

我们可以看到下面的输出:

* * * * * * * *代码建立了由Pulak茶室(181201)* * * * * * * *

* * * * * * * *代码使用Newton-Rapshon虚根的方法* * * * * * * *

输入变量x的函数形式:

x ^ 4 - (5 * x ^ 3) + (20 * x ^ 2) - (40 * x) + 60

你有初始近似值你的数学吗?

如果是的,请按1,如果没有按2:

输入小数:

你的帐号多项式方程是:

一个=

x ^ 4 - 5 * x ^ 3 + 20 * x ^ 2 - 40 * x + 60

通过使用Bewton-Raphson方法,我们得到了

我根是:-9.5195 + 0.0000

* * * * * * * *代码建立了由Pulak茶室(181201)* * * * * * * *

但我真正的答案是1.915 + 1.908。

现在我能做些什么呢?

沃尔特·罗伯森 2020年12月28日

在您的代码中,当用户表示他们没有一个初始值,您可以使用一个默认的初始值为1。

牛顿迭代法使用的比例函数及其导数。给定函数总是产生真正的值时输入,和真正的差异值总是真实的,所以我们可以看到后片刻的思考导数的定义限制斜坡上的差异,在一个区域,一个函数产生真正的输出为输入,它遵循导数也必须实值。然后把两个真实值的比值,结果将是实值。

因此,当你使用牛顿迭代函数产生真正的输出为输入,它遵循牛顿迭代方法将只项目到一个新的实值位置和永远不会项目复值的位置。因此当你的函数产生真正的输出为输入和使用实值默认输入,牛顿迭代方法永远找不到一个复杂的根源。

函数停止运行,因为你达到迭代极限,不是因为它到达错误宽容。

沃尔特·罗伯森 2020年12月28日

改变

                                   f =内联(a)
                                  
                                   dif = diff(信谊(a));
                                  
                                   d =内联(dif);

来

                                   f = matlabFunction (,“var”,x);
                                  
                                   d = matlabFunction (diff),“var”,x);

登录置评。

类别

MATLAB 编程功能函数创建

找到更多的在函数创建在帮助中心和文件交换

社区寻宝

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

开始狩猎!

我怎么能解决下面的问题?

0评论
显示1年长的评论藏1年长的评论

接受的答案

11日评论
显示十年长的评论隐藏10年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

我怎么能解决下面的问题?

0评论 显示1年长的评论藏1年长的评论

接受的答案

11日评论 显示十年长的评论隐藏10年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

0评论
显示1年长的评论藏1年长的评论

11日评论
显示十年长的评论隐藏10年长的评论