規劃問題(一)

2022-06-24 21:39:09 字數 2429 閱讀 4227

[x,y]=linprog(f,a,b,aeq,beq,lb,ub)

其中,x是最優解時自變數的取值,y是最優解,f是價值向量,就是目標函式的係數,a是線性不等式的約束的係數矩陣,b是不等式的的另一側的係數矩陣,aeq是等式的約束的係數矩陣,beq是等式約束的另一側的係數矩陣,lb是x上界,ub是下界。

其中,求的是最小值。具體使用辦法:

整數規劃

min∑cij·xij

s.t.:第j列只有乙個人是1,第i行只有乙個是1.

用法和linprog幾乎一樣,也是求最小值。

指派問題2.7,

如何把二維決策變數變成一維決策變數?  c=c(:)

其他的看不懂了!這個a到底是什麼鬼啊!!!

例2.8注意混合整數規劃中,lb ub怎麼寫。

非線性規劃

約束極值問題

例題:飛行管理問題

function [output 1,....] = name(input1,...)

呼叫的話,函式一定要放在當前目錄下,要不然就報錯說「未定義函式或變數」。

v是變數,n是階數。

對於一些簡單的運算:

;simplify:

各種化簡↑;

pretty 分數線居中顯示。

對函式求導=0即可:

solve函式

syms x

y=balabalaba;

dy=diff(y);

dy_zero=solve(dy);

dy_zero_num=double(dy_zero);%變成數值型別

ezplot(y);%畫圖

此處求得的是x軸座標,就算求出來的點也是駐點,不一定是極值點。再利用高數的知識:

hessian矩陣的特徵值就是形容其在該點附近特徵向量方向的凹凸性,特徵值越大,凸性越強;

表示將符號表示式s中的符號變數old替換為新的值new。

syms x

syms y

f=x^3-y^3+3*x^2+3*y^2-9*x;

df=jacobian(f); %求一階導

d2f=jacobian(df); %求hessian陣

[xx,yy]=solve(df) %求駐點

xx=double(xx);yy=double(yy);

for i=1:length(xx)

a=subs(d2f,,);

b=eig(a); %求特徵值

f=subs(f,,);f=double(f);

if all(b>0)

fprintf('(%.1f,%.1f)是極小值點,對應的極小值為%.1f\n',xx(i),yy(i),f);

elseif all(b<0)

fprintf('(%.1f,%.1f)是極大值點,對應的極大值為%.1f\n',xx(i),yy(i),f);

elseif any(b>0)&any(b<0)

fprintf('不是極值點\n');

else

fprintf('無法判斷\n');

endend

以上是符號解,以下是數值解

另一種 fminunc函式和fminsearch函式等等等等一堆函式;

還有梯度法。

還能用optimtool,這個感覺就是上面**的另一種圖形化形式而已,不看了(其實是不會哈哈哈)。

不定積分:int(expr,v)

定積分:int(expr,v,a,b)

symsum(expr,v,a,b)

見台大郭老師ppt。

和c很像。

初值問題是指在自變數的某值給出適當個數的附加條件,用來確定微分方程的特解的這類問題

syms y(x)

dsolve(x.^2+y+(x-2*y)*diff(y)==0)

求高階微分方程:

後話:參考書是《數學建模演算法與應用》,以上隨手記錄,加強印象~

[chū zhí wèntí shì zhǐ zài zì biànliàng de mǒu zhí gěi chū shìdàng gè shǔ de fùjiā tiáojiàn, yòng lái quèdìng wéifēn fāngchéng de tè jiě de zhè lèi wèntí]

動態規劃之揹包問題(一)

一切都要從一則故事說起。話說有一哥們去森林裡玩發現了一堆寶石,他數了數,一共有n個。但他身上能裝寶石的就只有乙個揹包,揹包的容量為c。這哥們把n個寶石排成一排並編上號 0,1,2,n 1。第i個寶石對應的體積和價值分別為v i 和w i 擦,目測這哥們是一苦逼程式設計師 排好後這哥們開始思考 揹包總...

乙個有趣的規劃問題

最近看到了乙個有趣的小問題,整理了一些自己的思考,對此做乙個簡單的總結。問題如下 我在某寶買了100袋樂事大波浪薯片,薯片到了菜鳥驛站,菜鳥驛站距離我家500公尺,我要把薯片全部拿回家。但是,我每次最多只能拿50袋,並且我嘴饞,只要手上有薯片,每走10公尺就要吃掉一袋。請問,我最多能把幾袋薯片拿回家...

分享非線性動態規劃問題(一)

摘要 前面的都是簡單的線性動規,本篇開始樹形動規問題的入門.1.正文 樹是我最喜歡的資料結構 沒有之一 它既不是像線性表那樣過於單調的線性關係,又沒有像圖那樣有著複雜的網狀關係,規律的非線性關係常常是我們現實生活中諸多問題的基本模型,比如家族關係 企業架構等。當它與動態規劃相結合時,似乎很高大上,但...