單純形法課程論文-資料下載頁
2024-10-29 02:25本頁面
【正文】 [z,i]。endz z(y)=[]。z Xtfor i=1:length(z)。yz=A(z(i),s(k))A(z(i),:)=A(z(i),:)+A(y,:).*yz B(z(i))B(y)yzB(z(i))=B(z(i))+B(y).*yz endfor i=1:length(Xt)if Xt(i)==y1 Xt(i)=s(k)。breakend end Xtdisp(39。170。187。187。186。243。39。)A=A B=B AB=[A,B]。for i=1:length(C)D=0。for j=1:length(Xt)D=D+AB(j,i)*C(Xt(j))。endxi(i)=C(i)D。end xi s=[]。for i=1:length(xi)1if xi(i)endend svpa([A,B。C])。f=length(s)。h=h+1。if h==5breakend endxi(length(xi))第五篇：單純形法C語言程序代碼長 春 工 業(yè) 大 學課程設(shè)計程序代碼課程設(shè)計名稱 運籌學課程設(shè)計 專 業(yè) 信息管理與信息系統(tǒng) 班 級 130506班 學 生 姓 名 于松南、張鑫蕊、趙改玲、趙海潮指 導 教 師王亞君、王忠吉2015年7月3日include includeint m。//記錄約束條件方程組的個數(shù) int n。//記錄未知量的個數(shù) float M=。float A[100][100]。//用于記錄方程組的數(shù)目和系數(shù)float C[100]。//用于存儲目標函數(shù)中各個變量的系數(shù) float b[100]。//用于存儲常約束條件中的常數(shù) float CB[100]。//用于存儲基變量的系數(shù) float seta[100]。//存放出基與入基的變化情況 float [100]。//存儲檢驗數(shù)矩陣 float x[100]。int num[100]。//用于存放出基與進基變量的情況 float Z=0。//記錄目標函數(shù)值void shuru()。void print()。int mincz()。int find_line(int a)。void exchange(int a,int b)。int main(){int i,j=0。int p,q,temp。//q：換入，p：換出shuru()。printf(“nn”)。printf(“ tCBtXBtbt”)。for(i=0。iprintf(“ X(%d)t”,i+1)。for(i=0。ix[i]=0。printf(“n”)。while(1){q=mincz()。if(q==1){print()。printf(“n所得解已經(jīng)是最優(yōu)解!n”)。printf(“n最優(yōu)解為:n”)。for(j=0。jtemp=num[j]1。x[temp]=b[j]。}for(i=0。iprintf(“x%d=%.2f ”,i+1,x[i])。Z=Z+x[i]*C[i]。}printf(“Z=%.2f”,Z)。break。}print()。p=find_line(q)。printf(“np=%d,q=%d”,p,q)。if(q==1)break。exchange(p,q)。}return 0。} int mincz(){int i,k=0。int flag=0。//檢驗數(shù)標記float min=0。for(i=0。iif([i]=0)flag=1。else {flag=0。break。}if(flag==1)return1。//進行到此處，說明存在//找到最小的檢驗數(shù)，作為換入變量for(i=0。iif(min[i]){min=[i]。k=i。}}return k。} int find_line(int a){int i,k,j。int flag=0。float min。k=a。for(i=0。iif(A[i][k]flag=1。else {flag=0。break。}if(flag==1){printf(“n該線性規(guī)劃無最優(yōu)解!n”)。return1。}for(i=0。iif(A[i][k]0)seta[i]=b[i]/A[i][k]。else seta[i]=M。}min=M。for(i=0。iif(min=seta[i]){min=seta[i]。j=i。}}num[j]=k+1。CB[j]=C[k]。return j。} void exchange(int p,int q){int i,j,c,l。float temp1,temp2,temp3。c=p。//行號，換出l=q。//列號，換入temp1=A[c][l]。//A[c][l]主元b[c]=b[c]/temp1。for(j=0。jA[c][j]=A[c][j]/temp1。//主元化為1for(i=0。iif(i!=c)if(A[i][l]!=0){temp2=A[i][l]。b[i]=b[i]b[c]*temp2。//主元所在列，其余元素化為0for(j=0。jA[i][j]=A[i][j]A[c][j]*temp2。}}temp3=[l]。for(i=0。i[i]=[i]A[c][i]*temp3。} void print(){int i,j=0。printf(“nn”)。for(i=0。iprintf(“%(%d)% ”,CB[i],num[i],b[i])。for(j=0。jprintf(“% ”,A[i][j])。printf(“n”)。}printf(“nn”)。printf(“ttt”)。for(i=0。iprintf(“ %”,[i])。printf(“nn”)。} void shuru(){int i,j。//循環(huán)變量int k。printf(“請輸入線性規(guī)劃問題的約束條件個數(shù)M：”)。scanf(“%d”,amp。m)。printf(“請輸入線性規(guī)劃問題的決策變量個數(shù)N：”)。scanf(“%d”,amp。n)。printf(“n請輸入方程組的系數(shù)矩陣A(%d行%d列):n”,m,n)。for(i=0。ifor(j=0。jscanf(“%f”,amp。A[i][j])。printf(“n請輸入初始基變量的數(shù)字代碼矩陣:n”)。for(i=0。iscanf(“%d”,amp。num[i])。printf(“n請輸入方程組右邊的值矩陣b:n”)。for(i=0。iscanf(“%f”,amp。b[i])。printf(“n請輸入目標函數(shù)各個變量的系數(shù)所構(gòu)成的系數(shù)陣C:n”)。for(i=0。iscanf(“%f”,amp。C[i])。for(i=0。i[i]=C[i]。for(i=0。ik=num[i]1。CB[i]=C[k]。}}
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號-1