| 查看: 2405 | 回复: 0 | |||
mayaya26木虫 (小有名气)
|
[求助]
请教大家一个MATLAB求解四元一次方程组的问题
|
|
编程中里面存在方程组fun_2，方程组中的未知量为 x y theta theta1 ，我尝试用简单的牛顿法求解，好像不行，求各位高手相助。谢谢 程序如下： function f=fun_2(x) d=0.4;epsilong=9; xt=-0.5;xr1=0.5;L=5; syms x y theta theta1 phi=asin(sin(theta)./sqrt(epsilong)); phi1=asin(sin(theta1)./sqrt(epsilong)); Lt=d/cos(phi); L1=d/cos(phi1); Ltt=(y-d)/cos(theta); L11=(y-d)/cos(theta1); f1=(x-xt)^2+y^2-(Lt^2+Ltt^2 +2*Lt*Ltt*cos(theta - phi)); f2=(x-xr1)^2+y^2-(L1^2+L11^2 +2*L1*L11*cos(theta1 - phi1)); f3=L-sqrt(epsilong)*Lt+Ltt+sqrt(epsilong)*L1+L11; f=[f1 f2 f3]; end function df=dfun_2(x) syms x y theta theta1 f=fun_2(x); df=[diff(f,'x');diff(f,'y');diff(f,'theta');diff(f,'theta1')]; df=conj(df'); end function x=newton_2(x0,eps,N) con=0; % syms x y theta theta1 %其中x0为迭代初值eps为精度要求N为最大迭代步数con用来记录结果是否收敛 for i=1:N; x0=[0 1 0.5 0.4];eps=1.0e-3;N=10; for i=1:N; f=subs(fun_1(x0),{'x' 'y' 'theta' 'theta1'},{x0(1) x0(2) x0(3) x0(4)}); df=subs(dfun_1(x0),{'x' 'y' 'theta' 'theta1'},{x0(1) x0(2) x0(3) x0(4)}); x=x0-f/df; for j=1:length(x0) il(i,j)=x(j); end if norm(x-x0)<eps con=1; break; end x0=x; end 运行结果： newton_2 Error using mupadmex Error in MuPAD command: Out of memory. Error in sym/privBinaryOp (line 1506) Csym = mupadmex(op,args{1}.s, args{2}.s, varargin{:}); Error in sym/mrdivide (line 237) X = privBinaryOp(A, B, 'symobj::mrdivide'); Error in newton_2 (line 9) x=x0-f/df; 是哪里出错了？  大神们帮帮忙，比较急 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
matlab 如何求解多元二阶微分方程组
已经有7人回复
- Maxwell方程如何求解
已经有10人回复
- 求解一个方程组，急！！！
已经有19人回复
- 请教高手，如何解非线性方程组！！！
已经有13人回复
- 请教matlab如何增加迭代次数
已经有11人回复
- matlab如何求解多元多次方程，包括绘图
已经有7人回复
- matlab求三元4次符号代数方程组，出现warning且答案不正确，该怎么办？
已经有12人回复
- matlab解方程组求助攻
已经有4人回复
- 简单四阶导数求解！！
已经有7人回复
- matlab解四次方程的问题
已经有6人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- MATLAB求解非线性方程组
已经有5人回复
- Mathematica 简单 四元一次方程组 ，求解问题~
已经有9人回复
- 初学者，请教大家关于matlab解多元二次方程组
已经有4人回复
- 求助matlab求解三角函数方程
已经有9人回复
- Matlab求解偏微分方程组
已经有22人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人