版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

mayaya26

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2243.2
散金: 104
紅花: 2
帖子: 132
在線: 115.5小時
蟲號: 1683262
注冊: 2012-03-11
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 請教大家一個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為最大迭代步數(shù)con用來記錄結(jié)果是否收斂 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

運行結(jié)果： 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;

是哪里出錯了？

大神們幫幫忙，比較急

回復(fù)此樓

» 猜你喜歡

2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！已經(jīng)有5人回復(fù)
298求調(diào)劑已經(jīng)有3人回復(fù)
數(shù)一英一271專碩（085401）求調(diào)劑，可跨已經(jīng)有8人回復(fù)
295求調(diào)劑已經(jīng)有4人回復(fù)
0856求調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有10人回復(fù)
299求調(diào)劑已經(jīng)有6人回復(fù)
一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑已經(jīng)有6人回復(fù)
289求調(diào)劑已經(jīng)有12人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有6人回復(fù)

1樓 2015-05-25 18:54:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

matlab 如何求解多元二階微分方程組已經(jīng)有7人回復(fù)
Maxwell方程如何求解已經(jīng)有10人回復(fù)
求解一個方程組，急�。。� 已經(jīng)有19人回復(fù)
請教高手，如何解非線性方程組�。�！已經(jīng)有13人回復(fù)
請教matlab如何增加迭代次數(shù) 已經(jīng)有11人回復(fù)
matlab如何求解多元多次方程，包括繪圖已經(jīng)有7人回復(fù)
matlab求三元4次符號代數(shù)方程組，出現(xiàn)warning且答案不正確，該怎么辦？已經(jīng)有12人回復(fù)
matlab解方程組求助攻已經(jīng)有4人回復(fù)
簡單四階導(dǎo)數(shù)求解！！已經(jīng)有7人回復(fù)
matlab解四次方程的問題已經(jīng)有6人回復(fù)
matlab如何求解一個非線性微分方程組已經(jīng)有8人回復(fù)
請教MATLAB中的LMI求解問題，急�。。� 已經(jīng)有7人回復(fù)
MATLAB求解非線性方程組已經(jīng)有5人回復(fù)
Mathematica 簡單四元一次方程組，求解問題~ 已經(jīng)有9人回復(fù)
初學(xué)者，請教大家關(guān)于matlab解多元二次方程組已經(jīng)有4人回復(fù)
求助matlab求解三角函數(shù)方程已經(jīng)有9人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【求助】求MAtlab求解程序！�。� 已經(jīng)有20人回復(fù)

點擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 mayaya26 的主題更新