查看: 144|回复: 2|关注: 0

[已解决] 求赐教,Fsolve 要求函数返回的所有值都是数据类型的 double

[复制链接]

新手

12 麦片

财富积分


050


3

主题

4

帖子

0

最佳答案
本帖最后由 GRotta 于 2020-10-12 11:19 编辑

我的程序:


  1. clc;
  2. clear;
  3. % r1=0.256;
  4. % r2=0.332;
  5. F=22500;
  6. E=6.312467.*10^9;
  7. g=38/73.5/pi;
  8. h=((g.^3+1)/2/(g.^2)).^0.5;
  9. e=(1-g.^2).^0.5;
  10. syms t1;
  11. y1=1/((1-(e*sin(t1)).^2).^0.5);
  12. Ke=int(y1,t1,0,pi/2);
  13. syms t2;
  14. y2=((1-(e*sin(t2)).^2).^0.5);
  15. Ee=int(y2,t2,0,pi/2);
  16. f1=Ee/(((1-e.^2).^0.5)*Ke);
  17. f2=pi*(Ee.^0.5)/(2*((1-e.^2).^0.5)*(Ke.^1.5));
  18. m=(4*E*h*f2/3/F).^2;
  19. global x;
  20. A=0:0.0001:0.005;
  21. n=length(A);
  22. for i=1:n
  23.     x=A(i);
  24.     [y(i),fval(i)]=fsolve(@(y) (y.^3)*m+acos((y-x)/x)-pi/2,0.0001);
  25. end
  26. [A;y];
  27. figure(1);
  28. subplot(1,1,1);
  29. plot(A,y);
复制代码



错误信息:


错误使用 fsolve (line 269)
FSOLVE requires all values returned by functions to be of
data type double.

出错 GIJIN1 (line 25)
    [y(i),fval(i)]=fsolve(@(y)
    (y.^3)*m+acos((y-x)/x)-pi/2,0.0001);


是公式里面的“m”设置的有问题吗








回复主题 已获打赏: 0 积分

举报

新手

12 麦片

财富积分


050


3

主题

4

帖子

0

最佳答案
 楼主| 发表于 2020-10-12 11:14:54 | 显示全部楼层
是我的m 处理的有问题吗
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

9

主题

1650

帖子

353

最佳答案
  • 关注者: 81
发表于 2020-10-12 21:12:27 | 显示全部楼层 |此回复为最佳答案
仅供参考
  1. clc;
  2. clear;
  3. % r1=0.256;
  4. % r2=0.332;
  5. F=22500;
  6. E=6.312467.*10^9;
  7. g=38/73.5/pi;
  8. h=((g.^3+1)/2/(g.^2)).^0.5;
  9. e=(1-g.^2).^0.5;
  10. syms t1;
  11. y1=1/((1-(e*sin(t1)).^2).^0.5);
  12. Ke=double(int(y1,t1,0,pi/2));
  13. syms t2;
  14. y2=((1-(e*sin(t2)).^2).^0.5);
  15. Ee=double(int(y2,t2,0,pi/2));
  16. f1=Ee/(((1-e.^2).^0.5)*Ke);
  17. f2=pi*(Ee.^0.5)/(2*((1-e.^2).^0.5)*(Ke.^1.5));
  18. m=(4*E*h*f2/3/F).^2;
  19. global x;
  20. A=0:0.0001:0.005;
  21. n=length(A);
  22. for i=1:n
  23.     x=A(i);
  24.     [y1(i),fval(i)]=fsolve(@(y) (y.^3)*m+acos((y-x)/(x+eps))-pi/2,0.0001);
  25. end
  26. [A;y1];
  27. figure(1);
  28. subplot(1,1,1);
  29. plot(A,y1);
复制代码
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表