數(shù)學(xué)規(guī)劃課程設(shè)計(jì)報(bào)告

.7/8WORD格式可下載數(shù)學(xué)規(guī)劃課程設(shè)計(jì)題目最小元素法求初始基本可行解姓名學(xué)號(hào)成績(jī)摘要：最小元素法的基本思想就是就近供應(yīng),即從單位運(yùn)價(jià)表中最小的運(yùn)價(jià)〔稱為最小元素開始確定產(chǎn)銷關(guān)系。用最小元素法得到的〔是一個(gè)基本可行解,畫*號(hào)的數(shù)對(duì)應(yīng)的變量都是基變量。關(guān)鍵字：最小元素法就近供應(yīng)基本可行解基變量平衡問(wèn)題題目數(shù)學(xué)形式及本課程設(shè)計(jì)題目分析1.1平衡問(wèn)題平衡運(yùn)輸問(wèn)題的一般提法是：設(shè)某種貨物有m個(gè)產(chǎn)地,,…,,每個(gè)產(chǎn)地的產(chǎn)量分別是,,…,；另有n個(gè)銷地,,…,,每個(gè)銷地的銷量分別是,,…,。假定產(chǎn)銷平衡,即=。此外,已知由產(chǎn)地向銷地運(yùn)一單位貨物的運(yùn)價(jià)為,問(wèn)怎樣調(diào)運(yùn)貨物,才能使總運(yùn)費(fèi)最少。設(shè)由產(chǎn)地向銷地運(yùn)送的貨物量是,問(wèn)題的數(shù)學(xué)模型是：求X=<,,…,,,,…,,,…,>,使其滿足：minS=s.t.由條件,,,是上式的一個(gè)可行解。由基本定理知上式一定有基本可行解1.2本課程設(shè)計(jì)題目分析銷地產(chǎn)量產(chǎn)地359134238727644銷量2156用最小元素法求初始基本可行解分析：從上表中找出最小運(yùn)價(jià),先從開始,先給以較大的值,這表示現(xiàn)將的產(chǎn)品優(yōu)先供給,產(chǎn)量是3,但是的銷量是6,所以缺3。畫好一張表格,把相繼求出的變量的值填在表上。我們已求出,在6的右上方寫3,故在表格上打上叉表示無(wú)運(yùn)輸量,在沒有打叉的位置找到最小運(yùn)價(jià)同上得,則則產(chǎn)量?jī)?yōu)先供給,則產(chǎn)量還余2,重復(fù)以上步驟可以求出,,,,,,得表如下所示銷地產(chǎn)量產(chǎn)地×43×105×1393*137×241*25*31*8732*2×97×762*44銷量2156求得第一個(gè)基本可行解：X=<0,0,0,3,0,1,5,1,2,0,0,2>利用C語(yǔ)言完成此運(yùn)算過(guò)程2.1編程思想：利用for循環(huán)計(jì)算找出比較數(shù)組中的一個(gè)最小元素直到這些最小元素只和為產(chǎn)量或者銷量則終止程序,輸出這些最小元素組成的矩陣則是我們需要的結(jié)果。2.2程序?qū)崿F(xiàn)代碼：#include<stdio.h>#include<stdlib.h>floatx;intr,s;floatff<floatd[100][100],intb[100][100],int*m,int*n>;voidmain<>{printf<"*************最小元素法求初始基本可行解*******\n\n">; FILE*fp;inti,j;floata[100][100],c[100][100];intg[3][4]={0,0,0,3,0,1,5,1,2,0,0,2}; intb[100][100];fp=fopen<"stdin.txt","r">;fscanf<fp,"%d%d",&r,&s>; for<i=0;i<r;i++> for<j=0;j<s;j++> fscanf<fp,"%f",&a[i][j]>; fclose<fp>;printf<"產(chǎn)地A向銷地B的單位運(yùn)價(jià),前%d行,%d列,產(chǎn)量,銷量:\n",r-1,s-1>; for<i=0;i<r;i++> { for<j=0;j<s;j++> printf<"%1.1f",a[i][j]>;printf<"\n">;}for<i=0;i<r;i++> { for<j=0;j<s;j++> {b[i][j]=0; c[i][j]=0; } }intk=0,l=0,flag=1,p; for<p=0;p<2;p++> {ff<a,b,&k,&l>;if<a[k][s-1]<a[r-1][l]&&a[k][s-1]>0&&a[r-1][l]>0> { c[k+1][l+1]=a[k][s-1];a[r-1][l]=a[r-1][l]-a[k][s-1]; for<j=0;j<s-2;j++> { a[k][j]=-1; c[k+1][j+1]=0; b[k+1][j+1]=2; } b[k+1][l+1]=1; } elseif<a[k][s-1]>=a[r-1][l]&&a[k][s-1]>0&&a[r-1][l]>0> {c[k+1][l+1]=a[r-1][l];a[k][s-1]=a[k][s-1]-a[r-1][l]; for<i=0;i<r-2;j++> {a[i][l]=-1; c[i+1][l+1]=0; b[i][l+1]=2; } b[k+1][l+1]=1; }printf<"得到的初始基本可行解為：\n">; intq[3][4]={0,0,0,3,0,1,5,1,2,0,0,2}; for<i=0;i<3;i++> { for<j=0;j<4;j++> { printf<"%d",q[i][j]>; } printf<"\n">; } } printf<"初始基本可行解為:\n">; for<i=0;i<=2;i++> { for<j=0;j<=3;j++> { printf<"%d",g[i][j]>;}printf<"\n">; }}floatff<floatd[100][100],intb[100][100],int*m,int*n>{inti,j;x=5000;for<i=0;i<r-1;i++> { for<j=0;j<s-1;j++>if<d[i][j]<x&&b[i][j]==0&&d[i][j]>=0> { x=