Question

Matlab 求值　fzero　第二个值怎么也求不出来　　请大师帮忙

clear
>> x=[0:0.5:5];
>> y=2.0095*sqrt(x).*tan(0.7265*sqrt(x))-sqrt(5.6976-0.3205*x);
>> plot(x,real(y),'r','LineWidth',1.5);
>> h=line([4.6748,4.6749],[0,0]);
>> h=line([0,5],[0,0]);
>> hold on
>> plot(x,imag(y),'bo');
>> grid
f=@(x)real(2.0095*sqrt(x).*tan(0.7265*sqrt(x))-sqrt(5.6976-0.3205*x));
[x1,f1,exitflag1]=fzero(f,1);
[x2,f2,exitflag2]=fzero(f,4);
x=[x1,x2];
f=[f1,f2];
>> x
x =
1.2228 4.6749
>> f
f =
1.0e+016 *
0.0000 1.5338
>> x1
x1 =
1.2228
>> f1
f1 =
4.4409e-016
>> x2
x2 =
4.6749
>> f2
f2 =
1.5338e+016
>>
第二个值明显错误，请大师帮忙
两个人答的都很好，奖励谁？

Answer 1

你的错误在于，x0的初值选择不正确。因初值离真值较远，所以出现错误。

从图形看，有无数个解。

第一个解：

>> [x,fval,exitflag]=fzero(f,1)
x =
    1.2228
fval =
  4.4409e-016
exitflag =
     1

第二个解：

>> [x,fval,exitflag]=fzero(f,20)
x =
   18.6994
fval =
  2.7948e-015
exitflag =
     1

第三个解：

>> [x,fval,exitflag]=fzero(f,80)
x =
   74.7977
fval =
  1.1179e-014
exitflag =
     1

。。。

追问

你这图怎么画出来的？发shi_w_y@我163的，谢谢，

追答

用 ezplot（）函数画出来。
>> ezplot('y=2.0095*sqrt(x).*tan(0.7265*sqrt(x))-sqrt(5.6976-0.3205*x)',[-10 50 -100 100]),grid on

Answer 2

程序选的x太少了。

 x=[0:0.05:5];

>> y=2.0095*sqrt(x).*tan(0.7265*sqrt(x))-sqrt(5.6976-0.3205*x);

>> plot(x,real(y),'r','LineWidth',1.5);

>>  h=line([0,5],[0,0]);

>> grid on

 f=@(x)real(2.0095*sqrt(x).*tan(0.7265*sqrt(x))-sqrt(5.6976-0.3205*x));f=@(x)real(2.0095*sqrt(x).*tan(0.7265*sqrt(x))-sqrt(5.6976-0.3205*x));

>> [x2,f2,exitflag2]=fzero(f,4);

>> x2

x2 =

    4.6749

结果正确

追问

不对，f2= 1.5338e+016;也太不趋近于0了，能对呢？

追答

之前你画图只选用了10个点，也太少了。
从图片的数值上看，x2求得是正确的。
可能是fzero求解时，f2的计算是不对的。

追问

有这种情况？那也不正常，把x2代入求解也一样，怎么能不对了？

追答

你的函数中有tan，有奇点啊。失误，之前没有想到这点。
只是楼主，要求的x范围是[0 5]的话，便不用考虑后续的解了。
画图就是描点画图啊。分开画就可以，显示出这个分开的效果。
x=[0:0.0005:4.6749];
>> y=2.0095*sqrt(x).*tan(0.7265*sqrt(x))-sqrt(5.6976-0.3205*x);
plot(x,y)
>> x=[4.6749:0.0005:5];
>> hold on
>> y=2.0095*sqrt(x).*tan(0.7265*sqrt(x))-sqrt(5.6976-0.3205*x);
plot(x,y)

根据0.7265*sqrt(x)=pi/2+k*pi

可以解出来奇值点的坐标

追问

怎么解出奇值点？？？？

追答

添加命令控制坐标系 axis([4 5 -1000 1000])

>> k=0:10;

>> x=(pi/2+k*pi).^2./0.7265^2;x=x';

便是奇点

x =

  1.0e+003 *

    0.0047

    0.0421

    0.1169

    0.2291

    0.3787

    0.5657

    0.7901

    1.0518

    1.3510

    1.6876

    2.0616

解释：

第一个奇点：

x(1)

4.6749

第二个：

x(2)

42.0737

第三个

x(3)

116.8715

追问

还是楼上专业，不过通过你我也学了不少东西，这奖励给谁？

追答

奖励随意，哈哈，楼主可以抛硬币决定。答题的过程也是学习的过程。