1樓：
可以用fsolve直接求出。計算方法如下：
1、%p=x(1);q=x(2);
%例如資料x隨機給出
x=round(100*rand(10,1))
n=length(x);
f=@(x)[sum(x.^x(1).*log(x))./sum(x.^x(1))-1/x(1)-sum(log(x))/n;
x(2)-(sum(x.^x(1))/n).^(1/x(1))];
options=optimset('maxfunevals',2000,'maxiter',1000);
x = fsolve(f,[0.1;5],options)
p=x(1)
q=x(2)
2、由於你的方程組的特殊性，兩個變數之間沒有耦合，q沒有出現在第一個方程中。於是更加方便的，能每次求出的結果的解法是：
clear
%p=x
x=round(100*rand(10,1))
n=length(x);
f=@(x)sum(x.^x.*log(x))./sum(x.^x)-1/x-sum(log(x))/n;
options=optimset('maxfunevals',2000,'maxiter',1000);
[x fval]= fsolve(f,0.1,options)
p=xq=(sum(x.^x)/n).^(1/x)
個人推薦用第二種方法。因為解方程組，有事fsolve會出現求解失敗。
希望對你有幫助，如若解決，望採納
2樓：穗子和子一
非線性方程組數值解法 - 正文
n個變數n個方程(n >1)的方程組表示為
(1)式中ƒi(x1,x2,…,xn)是定義在n維歐氏空間rn 的開域d上的實函式。若ƒi中至少有一個非線性函式，則稱(1)為非線性方程組。在rn 中記 ƒ＝ 則(1)簡寫為ƒ(尣)=0。
若存在尣*∈d，使ƒ(尣*)＝0，則稱尣*為非線性方程組的解。方程組(1)可能有一個解或多個解，也可能有無窮多解或無解。對非線性方程組解的存在性的研究遠不如線性方程組那樣成熟，現有的解法也不象線性方程組那樣有效。
除極特殊的方程外，一般不能用直接方法求得精確解，目前主要採用迭代法求近似解。根據不同思想構造收斂於解尣*的迭代序列(k=0，1，…)，即可得到求解非線性方程組的各種迭代法，其中最著名的是牛頓法。
matlab中如何求解符號變數的非線性方程組的解，求大神賜教！！！**如下： 65
3樓：匿名使用者
matlab中如何求解符號變數的非線性方程組的解，對於簡單明瞭的方程組可以用solve()求得。solve()求解格式為
eq1=f1(x,y,z),eq2=f2(x,y,z),eq3=f3(x,y,z)
solve(eq1,eq2,eq3)
但分析了你給出的**，我覺得只能用數值的方法求解，得到其數值解。你想要解析值是有一定的難度。
常用的數值方法，有二分法，牛頓法等等。
用matlab解多元非線性方程組，求大神
4樓：僅僅是追憶
首先定義函式:
function f=fx(x)
f(1)=3*x(1)-cos(x(2)*x(3))-0.5;
f(2)=x(1)^2-81*(x(2)+0.1)^2+sin(x(3))+1.06;
f(3)=exp(-x(1)*x(2))+20*x(3)+(10*pi-3)/3;
>> x=fsolve(@fx,[1,1,1])最後求得
x =0.5000 0.0000 -0.5236
5樓：匿名使用者
很高興為您解答，首先定義函式myfun
function f = myfun(x)f = [3*x(1)-cos(x(2)*x(3))-0.5;
x(1)^2-81*(x(2)+0.1)^2+sin(x(3))+1.06;
exp(-x(1)*x(2))+20*x(3)+(10*pi-3)/3];
前面這段函式定義在myfun.m檔案中
主程式如下：
[x,fval]=fsolve(@myfun,[0,0,0]);
x(1)
x(2)
x(3)
最後求得
x(1)=0.5000
x(2)=0.0000
x(3)=-0.5236
滿意請採納回答，謝謝！
6樓：匿名使用者
建立 myfun.m 檔案
function f = myfun(x,a)
e = a(1);
i = a(2);
r0 = a(3);
r1 = a(4);
t = a(5);
a = a(6);
v = a(7);
rho = a(8);
f = [ (t - rho * a * v^2) *
sin(x(3)) * x(1) - (t * cos(x(3)) + rho * a * v^2 - rho * a * v^2 *
cos(x(3))) * x(2) - e*i/(r0 + r1);
(1/3) * (t - rho * a * v^2) *
sin(x(3)) * x(1)^3 - (1/2) * (t * cos(x(3)) + rho * a * v^2 - rho *
a * v^2 * cos(x(3))) * x(2) * x(1)^2 - e* i * x(2);
(t - rho * a * v^2) * sin(x(3)) *
x(1)^2 - (t * cos(x(3)) + rho * a * v^2 - rho * a * v^2 *
cos(x(3))) * x(2) * x(1) - e* i * x(3)];
建立一個執行檔案
clcclear
a = zeros(8);
display('# pls input the known
parameters: #')
a(1) = input('e = ');
a(2) = input('i = ');
a(3) = input('r0 = ');
a(4) = input('r1 = ');
a(5) = input('t = ');
a(6) = input('a = ');
a(7) = input('v = ');
a(8) = input('rho = ');
display('# pls input the initial
point: #')
x0 = zeros(3);
% make a
starting guess at the solution
x0(1) = input('x1 = ');
x0(2) = input('y1 = ');
x0(3) = input('phi = ');
options =
optimset('display','iter');
% option to display
output
[x,fval] = fsolve(@(x)
myfun(x,a),x0,options)
% call solver
執行，輸入已知的幾個引數，再輸入初始搜尋點，即可！
7樓：匿名使用者
fx.m
function f=fx(x)
f(1)=3*x(1)-cos(x(2)*x(3))-0.5;
f(2)=x(1)^2-81*(x(2)+0.1)^2+sin(x(3))+1.06;
f(3)=exp(-x(1)*x(2))+20*x(3)+(10*pi-3)/3;
>> x=fsolve(@fx,[1,1,1])x =0.5000 0.0000 -0.5236
8樓：瑞秋老大
1.首先定義函式myfun
function f = myfun(x)
f = [3*x(1)-cos(x(2)*x(3))-0.5;
x(1)^2-81*(x(2)+0.1)^2+sin(x(3))+1.06;
exp(-x(1)*x(2))+20*x(3)+(10*pi-3)/3];
2.前面這段函式定義在myfun.m檔案中。
主程式如下：
[x,fval]=fsolve(@myfun,[0,0,0]);
x(1)
x(2)
x(3)
3.最後求得：
x(1)=0.5000
x(2)=0.0000
x(3)=-0.5236
擴充套件資料
matlab的簡介
matlab是美國mathworks公司出品的商業數學軟體，用於演算法開發、資料視覺化、資料分析以及數值計算的高階技術計算語言和互動式環境，主要包括matlab和simulink兩大部分。
非線性方程組數值解法
20世紀60年代中期以後，發展了兩種求解非線性方程組(1)的新方法。一種稱為區間迭代法或稱區間牛頓法，它用區間變數代替點變數進行區間迭代，每迭代一步都可判斷在所給區間解的存在惟一性或者是無解。
求大神幫我看一下這個matlab程式怎麼改
yyy名 quadgk比quad好用一些 clear all clc fi 5.pi.180 f 25400 e 5 vx 10 w1 370.2.pi.60 w2 100.2.pi.60 r0 50.5 h 0.0005 a 0.00005 振幅alb 0.00001 abc 1 25.4 f1 ...
matlab急求大神幫忙跪求大神幫忙做個多元線性規劃的matlab模型
x 128453,102398,28774,90302,128.6,64.1,64.5 129227,116694,28559,90976,124.1,64.0,60.1 129988,136515,27947,92184,122.9,64.2,58.7 130756,182321,26504,94...
求日語大神給下面這首歌歌詞的漢字標上假名不要羅馬音謝謝
透 通 水槽 世界 墨 掬 魚 口 浸 金 買 事実 並 模型 馴 合 踏 臺 痛 溺 魚 水槽 中 鳥 鳥 窮屈 空 僕 気 愛 憎 消 消 暗 映 綺麗 僕 目 悪 殘 溫 値 消 消 暗 溺 魚 水槽 中 上 積 偉 英雄 名乗 飛行機 落 死 數字 明日 忘 幼馴染 遊 子 自殺 名前 思 ...