尋求幫助以優化 matlab 程式，新手問題。

你很困惑。

函式 kutta % 最好有輸入和輸出，例如函式 [y1，y2]=kutta（x1，x2）。

sym tt=

以這種方式將 t 定義為符號變數基本上失去了意義。

u1=inv(kd+kt)*f

如果使用 u1，u1=inv（kd+kt）*（f+kt*u1） 幾乎毫無意義。

點選這個迴圈...

我基本上不明白該按什麼週期。

為什麼我在MATLAB中更改了目標函式，優化結果沒有變化，出了什麼問題？

1.題主給出的**不能得到最優解，方法錯誤。

2.分析題主給出的具體問題，實際的非線性函式優化問題應使用fmincon（）函式。 求解方法的思路是：

1. 自定義目標函式 myfun（x）。 即。

x1=x(1);y1=x(2);z1=x(3);

f=(1100962500*x1*y1*pi*(500+z1)+2652318750000*z1*pi);

2. 自定義約束函式 mycon（x）。 即。

eq = symsum(16*y1/(500+m*y1)*(500+z1+m*y1)^2/((500+z1+m*y1)^2-500^2),m,1,n);求和。

約束。 3.確定初始值，x0=[6 8 8]。

4.確定x的上下限，即lb=[6,8,8]; ub=[12,20,20];

5. 使用 fmincon（） 函式查詢 x、y 和 z 的值。

6. 驗證約束是否接近 0

按照上述方法程式設計，執行即可得到以下結果。

th1=zeros(36,2);

th1(:,2) =1;

為了判斷優化效果，我們每輪執行100萬次兩種方法，執行10輪，比較執行時間，結果如下：

T1 是優化前的時間，T2 是優化後的執行時間。

可以看出，優化前每一輪的平均執行時間約為秒，優化後每輪的平均執行時間約為秒。

主題的優化問題可以用 fmincon（） 函式解決。 解決方法如下：

x0=[ 0 0]

a=;b=;aeq=[1,1,1,1];beq=[1];

lb=[0,0,0,0];ub=[1,1,1,1];

x,fval,exitflag]=fmincon(@func,x0,a,b,aeq,beq,lb,ub,@fcon)

執行結果。

這個優化問題可以用 fmincon（） 獲得，結果 x =

fval =

目標函式：函式 f = myfun（x）。

f=;結束約束函式：

function [c,ceq] = mycon(x)c(1)=;

c(2)=;

c(3)=750*1728-pi*x(3)^2*x(4)-4*pi*x(3)^3/3;

c(4)=x(4)-240;

c(5)=;

c(6)=x(2);

c(7)=10-x(3);

c(7)=x(4)-200;

ceq=;end

錯誤的主要原因是所有乘法和除法運算都沒有更改為點乘法和除法。

建議的更改如下：

clearclc

q1=1;q2=-4;d=2;

x,y]=meshgrid(:5,:5);

z=(q1*(y-d/2)./(sqrt(x.^2+(y-d/2).^2)))q2*(y+d/2)./(sqrt(x.^2+(y+d/2).^2)))

k=:5;s,h]=contour(x,y,z,k);

結果如下圖所示：

您好，答案如下。

x = [200 250 300 400 500];

y = [

plot(x,y,'o')

hold on

z = polyfit(x,y,1)

y2 = polyval(z,x);

plot(x,y2,'r')z =