MATLAB实验报告第二次第三次.doc

MATLAB数学实验报告姓名： 李犁学号： 1304120724指导老师：易昆南专业班级：统计1203目录1学生管理程序 32自测题 83动画函数的生成.104 MATLAB在高等数学中的应用.125利用MATLAB求的值.146零件的参数设计储.167油罐的变位识别与罐容表标定.21一.学生管理程序 学号1304120724班级统计1203姓名李犁指导教师易昆南实验题目学生成绩管理程序评 分1.序言： 这个学期，我们学习了一种全新的计算机软件-MATLAB,MATLAB 是一种与数学密切相关的算法语言，是目前在工程界流 行最为广泛的计算机语言，也是当今最优秀的科技应用软件之一。它可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。在上过几次MATLAB实验后，我们已经初步掌握MATLAB的知识，而在这次作业部署中，我们要独立完成的是制作一个学生成绩管理系统，以方便老师查询学生成绩，这是更是对我们上课听课效果的考验。2.实验内容： 利用MATLAB编写一个学生管理程序，包括学号、姓名、性别、年龄、成绩，及成绩的分类，及格率，各分数段所占比例，并绘制相关的图像。3.设计（实习）目的：1.了解MATLAB在实际问题中的应用2.通过实践加深对这门语言中M文件的了解3.熟悉简单程序结构，如循环结构（for循环、while循环）选择结构（if-else-if）、分支语句（switch-case-otherwise）。4.实验构思： 运用MATLAB语言，进行程序设计，运用MATLAB在带函数实现对成绩的分类，计算等。在这个成绩管理系统，我们要将学生的各个成绩进行分类。将不同成绩段的分为“满分”、“优秀”、“较好”、“及格”、“不及格”、“缺考”，主要是运用循环语句，十七每个分数都被标记。然后在从成绩中选出最高分，最低分，平均分，还有方差。求各分数段人数，及格率。在这个环节中主要运用函数的调用，因为MATLAB含有丰富的函数资源，这为我们编程者提供了很多便利。而在最后我又调用bar函数，绘制直方图，更直观的观察学生的成绩情况。调用legend函数，绘制饼形图。5：具体操作：clear;for i=1:10;ai=89+i;%把90 分及以上的学生成绩赋值给abi=79+i;%把80 分及以上的学生成绩赋值给bci=69+i;%把70 分及以上的学生成绩赋值给cdi=59+i;%把60 分及以上的学生成绩赋值给dend;c=d,c;StudentNo=1,2,3,4,5,6,7,8,9,10,11,12;%学生学号Name=张三,赵四,王五,马六,聪聪,宝宝,刘谦,二蛋,二胖,鸣笛,泡泡,吉吉;%学生姓名Sex=男,女,男,男,男,男,男,女,女,女,男,男;%学生性别Age=17,16,17,18,16,15,18,17,19,15,17,16;%学生年龄Mark=91,73,66,0,80,100,76,82,92,77,86,59;%学生成绩PostScript=cell(1,12);S=struct(StudentNo,StudentNo,Name,Name,Sex,Sex,Age,Age,Marks,Mark,PostScript,PostScript);%for i=1:6 num(i)=0;endfor i=1:12%循环12次 switch S(i).Marks%分支条件 case 100 S(i).PostScript=满分 ;num(1)=num(1)+1; %当100 时备注显示满分 case 0 S(i).PostScript= 缺考;num(6)=num(6)+1; %当0 时备注显示缺考 case a S(i).PostScript=优秀; num(2)=num(2)+1; %当满;足a 情况时备注显示优秀 case b S(i).PostScript=比较好;num(3)=num(3)+1; %当满足b 情况时备注显示比较好 case c S(i).PostScript=已及格;num(4)=num(4)+1; %当满足c 情况时备注显示已及格 otherwise S(i).PostScript=不及格;num(5)=num(5)+1; %其他情况显示不及格 endend%获取成绩中最大值，最小值，平均值，标准差s=91,73,66,0,80,100,76,82,92,77,86,59;N=0;for i=1:4 N=N+num(i);endmi=min(s);ma=max(s);mea=mean(s);fc=std(s);jg=N/12;disp(学号 , 学生姓名 , 性别, 年龄, 成绩 , 备注 );disp()%声明显示for i=1:12; disp(num2str(S(i).StudentNo),blanks(6),S(i).Name,blanks(6),S(i).Sex,blanks(6),num2str(S(i).Age),blanks(6),num2str(S(i).Marks),blanks(6),S(i).PostScript);%按循环依次显示姓名 成绩 备注end;disp(各个分数段人数统计:) disp(满分:,num2str(num(1) disp(缺考:,num2str(num(6) disp(优秀:,num2str(num(2) disp(比较好:,num2str(num(3) disp(已及格:,num2str(num(4) disp(不及格:,num2str(num(5) disp(及格率:,num2str(jg)disp(最小值:,num2str(mi)disp(最大值:,num2str(ma)disp(平均值:,num2str(mea)disp(标准差:,num2str(fc)h=figure(1);q=num(1),num(6),num(2),num(3),num(4),num(5); pie(q,0 0.1 0.2 0.3 0.4 0.5 ), legend(满分,缺考,优秀,比较好,已及格,不及格);h=figure(2);x=1,2,3,4,5,6,7,8,9,10,11,12;bar(x,s)set(gca,XTick,1:1:12)set(gca,XTickLabel,1,2,3,4,5,6,7,8,9,10,11,12)title(Sine Function);xlabel(Radians); ylabel(Function Value);colormap(color);6.实验结果：学号 学生姓名 性别 年龄 成绩 备注 1 张三 男 17 91 优秀2 赵四 女 16 73 已及格3 王五 男 17 66 已及格4 马六 男 18 0 缺考5 聪聪 男 16 80 比较好6 宝宝 男 15 100 满分 7 刘谦 男 18 76 已及格8 二蛋 女 17 82 比较好9 二胖 女 19 92 优秀10 鸣笛 女 15 77 已及格11 泡泡 男 17 86 比较好12 吉吉 男 16 59 不及格各个分数段人数统计:满分:1缺考:1优秀:2比较好:3已及格:4不及格:1及格率:0.83333最小值:0最大值:100平均值:73.5标准差:25.7777：实验总结 在学习MATLAB过程中我们要多动脑，多动手写程序、调试，看程序时候多想几个为什么，理解别人程序中精华的东西，多自己调调程序，可以改改程序，把知识转成自己的，弄懂一两个经典例子比粗粗的看一大书要好，值得一提是不是一下子就能看懂相对比较长，比较复杂的程序的，要做的是自己慢慢去分析，利用已经有的资料和帮助，等慢慢搞懂它，要多看几遍，相信多看一遍你就会有一些新的认识，一直到你已经掌握里例子的精华。 MATLAB是一款功能强大的软件，在这次实践中，我也是收获颇丰，刚开始面对这个问题，不知道该如何设计，但是努力着回忆老师每节课讲的内容，再加上对课本的查阅，渐渐有了思路，设计出了大概的框架。但是在具体编程过程中的确遇到了好多麻烦，自身原因还是因为对MATLAB的不熟悉，经过一番努力，终于设计出了一项相对简单的学生成绩管理系统。今后还有继续听老师认真讲课，在下次课题实践中设计更好的程序。 最后多读matlab高手写的程序，帖子，找到一个高手，在几个大的论坛可以搜索出一大堆的帖子，然后慢慢去看吧，从中可以学到很多东西,也许你会发现很多你都看不懂，这很正常，别人写程序考虑的是能实现这样的功能，运行效率等问题，而不是考虑是不是很容易能让别人理解，你就看你现在能看懂，再看些通过你断点调试，查help能看明白的，看不懂的就先放下，遇到好的和看不懂，你又很挺感兴趣的话题，你先收藏起来，过段时间你再回头看看，你会发现，每看一遍，你就会多懂一些东西，当你第三遍甚至第四遍看同一个程序的时候，也许你就会有豁然开朗的感觉。一个比较长的程序，首先不要怕长，把它分解成一句一句，一句一句看懂它。要明白一口吃不成胖子，好东西要慢慢体会，才能领会。 参考MATLAB语言实用教程 马莉编著 姓名：李犁 2013年10月20日2 自测题1.缺省的Current Directory窗口是什么？当前目录浏览器。2.命令窗口中输入x时，MATLAB的搜索顺序是什么？ 如果在MATLAB命令窗口中输入字符串xMATLAB编译器会按顺序尝试解释它 1 把x作为一个参数或变量去搜索2 检查x是否为一个内置的函数3 在默认的目录即当前工作目录Current Directory下搜索文件名为x.m的文件 4 在搜索路径的所有目录中搜索文件名为x.m的文件。 如果多个函数具有相同的名称那么只有搜索到的第一个函数会执行其他函数会被认为不存在而不被执行。3.如何扩展MATLAB的搜索路径？把程序所在目录扩展成MATLAB的搜索路径的方法：1.搜索路径对话框：利用file-set path菜单，进入设置搜索路径对话框2.path ：如path(path,e:myfiles); 把新目录扩展到搜索路径3.Addpath：扩展目录,如addpath e:myfiles end把新目录加到末尾；addpath e:myfiles begin把新目录加到开头。3.Genpath：得到所有搜索路径连接而成的长字符串。4.使用editpath或pathtool扩展目录。4.输入矩阵A和B的操作是： A=1 1 0;-1 1 1;1 -1 3B=1 0 1;0 1 1;-1 0 1A = 1 1 0 -1 1 1 1 -1 3B = 1 0 1 0 1 1 -1 0 15.执行A*B的结果是：A*Bans = 1 1 2 -2 1 1-2 -1 36.执行A.*B的结果是：A.*Bans = 1 0 0 0 1 1-1 0 37.执行A.B的结果是：A.Bans = 1 1 0 1 1 1 1 1 38.生成由A的主对角元素构成的对角阵的操作是：C=diag(diag(A)C = 1 0 0 0 1 0 0 0 39.赋值表达式A=A;1, 3, 5的结果是：A=A;1, 3, 5A = 1 1 0 -1 1 1 1 -1 3 1 3 510.从A中得到第二行第三列元素的操作是：D=A(2,:)E=A(3,:)D = -1 1 1E = 1 -1 311.从A中得到A的第一列元素构成的列向量的操作是：F=transpose(A(1,:)F = 1 1 012.A(:,1,3)的结果是：A(:,1,3)ans = 1 0 -1 1 1 3 1 5 13.把A还原成第1题中的A的操作是： A(4,:)=;a=A(:,1);b=A(:,2);A=a,1,1,-1,bA = 1 1 0 -1 1 1 1 -1 314.执行W1=zeros(2,3)的结果是：W1=zeros(2,3)W1 = 0 0 0 0 0 015.执行W2=ones(3)的结果是：W2=ones(3)W2 = 1 1 1 1 1 1 1 1 116.执行W1=eye(3,4)的结果是：W1=eye(3,4)W1 = 1 0 0 0 0 1 0 0 0 0 1 03 动画函数的生成1.编写语言的电影动画M=moviein(16);x=9*pi;for j=1:600 x(j+1)=x(j)-pi/50; y1=(x-sin(x)./(x+sin(x); y2=1+0.1; y3=1-0.1; plot(x,y1,b*,x,y2,r,x,y3,r,x,0,r); axis(0 9*pi -0.2 2.5) M(j)=getframe; endmovie(M,3)2.编写语言的电影动画。xlim(0,10);ylim(-1,1);x0=0;a=0.05;y1=-1:0.001:1;y2=0+0.05;y3=0-0.05;for n=0:0.01:10;hold on;y=sin(n)/n;plot(n,y,b.);if abs(n-x0)=aplot(x0-a,y1,r,x0+a,y1,r)endif abs(y)=aplot(n,y2,r,n,y3,r)endpause(0.1);end4. MATLAB在高等数学中的应用1.符号导数求函数的一阶导数和二阶导数输入：syms xdiff(xn,1)diff(xn,2)执行后得的一阶导数和二阶导数分别为：ans=xn*n/xans=xn*n2/x2-xn*n/x22.符号微分设，求。syms x yz=(1+x*y)y;diff(z,x)diff(z,y)则输出结果：ans=（1+x*y）y*y2/(1+x*y)ans=（1+x*y）y*(log(1+x*y)+y*x/(1+x*y)3. 符号积分求输入：syms xint（sin（x）/x,x）则输出为：ans=sinint（x）4.常微分方程求微分方程y+2xy=的通解。clear；dsolve（Dy+2*x*y=x*exp(-x2),x）便得到微分方程的通解：（1/2*x2+C1）*exp(-x2)5. 无穷级数求级数的和输入：syms xs2=symsum(x(3*k),k,1,inf)得到和函数s2=-x3/（x3-1）6. 线性方程组求解线性方程组：输入clear；A=1,1,-2,-1;3,-2,-1,2;0,5,7,3;2,-3,-5,-1;D=det(A)X=null(A)输出为：D = 0X = -0.4714 0.2357 -0.4714 0.70717. 利用积分求面积求由圆和双纽线，所围成的面积。syms thr1=3*cos(th);r2=1+cos(th);s1=int(1/2*r22,th,0,pi/3);s2=int(1/2*r12,th,pi/3,pi/2);S=2*(s1+s2)输出为：5/4*p5 利用MATLAB求的值方案1（分析方法）函数的幂级数展开（精确到）Arctanx的麦克劳林级数展开式为：取x=1得：=4arctan1=4（1-+-L）取前n项作为近似值，产生的误差为，要精确到，则要计算前二十万位，收敛速度太慢。在1时远比x=1收敛速度快，利用正切的倍角公式可得：这时可简单的利用来估计判断误差。具体代码如下：syms xx=0;n=10200;for i=1:n x=x+4*(-1)(i+1)/(2*i-1);endp=vpa(x);方案2（迭代加速算法）利用函数arctanx的泰勒展开式取x=1，则得到记 便可利用埃特金加速收敛方法来加快已知序列的收敛速度。方案3 （数值积分法）具体代码如下：syms xp=vpa(4*int(1/(1+x2),x,0,1);方案4 概率方法（蒙特卡洛随机模拟）具体代码如下：syms xx=0;n=500;for i=1:n a=rand(1,2); if a(1)2+a(2)2=1 x=x+1; endendx=4*x/n;方案5 （公式法） 具体代码如下：syms xf1=(-1)x/(2*x+1)*2(2*x+1);f2=(-1)x/(2*x+1)*3(2*x+1);p=vpa(4*(symsum(f1,x,0,inf)+symsum(f2,x,0,inf);方案6 （符号积分法）具体代码如下syms xf=4/(1+x2);y=int(f,x,0,1);numefic(y)方案7（Taylor展开法）具体代码如下:syms xf=(-1)x/(2*x+1);p=vpa(4*symsum(f,0,inf);方案8（刘徽割圆法）具体代码如下：syms xn=61;for i=1:nx=sqrt(2-sqrt(4-x2);endx=subs(x,x,1);p=3*2n*vpa(x,n+5);零件的参数设计61.问题重申：粒子分离器某参数（记作y）由7个零件的参数决定, 其经验公式为：y的目标值（记作yo）为1.50, 当y偏离yo0.1时，产品为次品，质量损失1000元。当y偏离yo0.3时，产品为废品，质量损失9000元。零件容差分别为A(1%)、B(5%)、C(10%)三个等级，容差为均方差的三倍，成本如表3.1.1所示，求零件参数标定值及零件容差等级所使总费用最小。 表3.1.1标定值容许范围C等B等A等0.075，0.125/25/0.225，0.3752050/0.075，0.12520502000.075，0.125501005001.125，1.87550/12，2010251000.5625，0.935/25100目标函数 （3.1.12）其中：达最小，即：用蒙特卡罗法，模拟，求目标函数最小。2. MATLAB程序设计用Matlab求出关于的偏导数：Fun1.mclear allclcsyms x1 x2 x3 x4 x5 x6 x7 y = 174.42*x1/x5*(x3/(x2 - x1). 0.85*(1 - 2.62*(1 - 0.36*(x4/x2). -0.56)1.5*(x4/x2)1.16)/(x6*x7)0.5;a1=diff(y,x1);a2=diff(y,x2);a3=diff(y,x3);a4=diff(y,x4);a5=diff(y,x5);a6=diff(y,x6);a7=diff(y,x7);x01=0.1;x02=0.3;x03=0.1;x04=0.1;x05=1.5;x06=16;x07=0.75;a11=subs(a1,x1 x2 x3 x4 x5 x6 x7,0.1 0.3 0.1 0.1 1.5 16 0.75);a22=subs(a2,x1 x2 x3 x4 x5 x6 x7,0.1 0.3 0.1 0.1 1.5 16 0.75);a33=subs(a3,x1 x2 x3 x4 x5 x6 x7,0.1 0.3 0.1 0.1 1.5 16 0.75);a44=subs(a4,x1 x2 x3 x4 x5 x6 x7,0.1 0.3 0.1 0.1 1.5 16 0.75);a55=subs(a5,x1 x2 x3 x4 x5 x6 x7,0.1 0.3 0.1 0.1 1.5 16 0.75);a66=subs(a6,x1 x2 x3 x4 x5 x6 x7,0.1 0.3 0.1 0.1 1.5 16 0.75);a77=subs(a7,x1 x2 x3 x4 x5 x6 x7,0.1 0.3 0.1 0.1 1.5 16 0.75);yy=subs(y,x1 x2 x3 x4 x5 x6 x7,0.1 0.3 0.1 0.1 1.5 16 0.75);Xr1=0.05;Xr2=0.1;Xr3=0.1;Xr4=0.1;Xr5=0.1;Xr6=0.1;Xr7=0.05;C1=25;C2=20;C3=20;C4=50;C5=50;C6=10;C7=25;(105/9)*(a11*Xr1*x01)2+(a22*Xr2*x02)2+(a33*Xr3*x03)2+(a44*Xr4*x04)2+.(a55*Xr5*x05)2+(a66*Xr6*x06)2+(a77*Xr7*x07)2)+105*(yy-1.5)2+200function f=fun(x)Yx(1)=(8721*(-x(3)/(x(1) - x(2)(17/20)*(-(131*(x(4)/x(2)(29/25)*(1 - 9/(25*(x(4)/x(2)(14/25)(3/2)/50 - 1)/(x(6)*x(7)(1/2)/(50*x(5) + (148257*x(1)*x(3)*(-(131*(x(4)/x(2)(29/25)*(1 - 9/(25*(x(4)/x(2)(14/25)(3/2)/50 - 1)/(x(6)*x(7)(1/2)/(1000*x(5)*(x(1) - x(2)2*(-x(3)/(x(1) - x(2)(3/20);Yx(2)=(8721*x(1)*(24759*x(4)*(1 - 9/(25*(x(4)/x(2)(14/25)(1/2)/(31250*x(2)2*(x(4)/x(2)(2/5) + (3799*x(4)*(x(4)/x(2)(4/25)*(1 - 9/(25*(x(4)/x(2)(14/25)(3/2)/(1250*x(2)2)*(-x(3)/(x(1) - x(2)(17/20)/(100*x(5)*x(6)*x(7)*(-(131*(x(4)/x(2)(29/25)*(1 - 9/(25*(x(4)/x(2)(14/25)(3/2)/50 - 1)/(x(6)*x(7)(1/2) - (148257*x(1)*x(3)*(-(131*(x(4)/x(2)(29/25)*(1 - 9/(25*(x(4)/x(2)(14/25)(3/2)/50 - 1)/(x(6)*x(7)(1/2)/(1000*x(5)*(x(1) - x(2)2*(-x(3)/(x(1) - x(2)(3/20);Yx(4)=-(8721*x(1)*(24759*(1-9/(25*(x(4)/x(2)(14/25)(1/2)/(31250*x(2)*(x(4)/x(2)(2/5)+(3799*(x(4)/x(2)(4/25)*(1-9/(25*(x(4)/x(2)(14/25)(3/2)/(1250*x(2)*(-x(3)/(x(1)-x(2)(17/20)/(100*x(5)*x(6)*x(7)*(-(131*(x(4)/x(2)(29/25)*(1-9/(25*(x(4)/x(2)(14/25)(3/2)/50-1)/(x(6)*x(7)(1/2);Yx(3)=-(148257*x(1)*(-(131*(x(4)/x(2)(29/25)*(1-9/(25*(x(4)/x(2)(14/25)(3/2)/50-1)/(x(6)*x(7)(1/2)/(1000*x(5)*(x(1)-x(2)*(-x(3)/(x(1)-x(2)(3/20);Yx(5)=-(8721*x(1)*(-x(3)/(x(1)-x(2)(17/20)*(-(131*(x(4)/x(2)(29/25)*(1-9/(25*(x(4)/x(2)(14/25)(3/2)/50-1)/(x(6)*x(7)(1/2)/(50*x(5)2);Yx(6)=(8721*x(1)*(-x(3)/(x(1) - x(2)(17/20)*(131*(x(4)/x(2)(29/25)*(1 - 9/(25*(x(4)/x(2)(14/25)(3/2)/50 - 1)/(100*x(5)*x(6)2*x(7)*(-(131*(x(4)/x(2)(29/25)*(1 - 9/(25*(x(4)/x(2)(14/25)(3/2)/50 - 1)/(x(6)*x(7)(1/2);Yx(7)=(8721*x(1)*(-x(3)/(x(1) - x(2)(17/20)*(131*(x(4)/x(2)(29/25)*(1 - 9/(25*(x(4)/x(2)(14/25)(3/2)/50 - 1)/(100*x(5)*x(6)*x(7)2*(-(131*(x(4)/x(2)(29/25)*(1 - 9/(25*(x(4)/x(2)(14/25)(3/2)/50 - 1)/(x(6)*x(7)(1/2);Yx=Yx(1) Yx(2) Yx(3) Yx(4) Yx(5) Yx(6) Yx(7);xi=0.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.00.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.00.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.010.010.010.010.010.010.010.010.010.010.010.010.010.010.010.010.010.00.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.050.010.010.010.010.010.010.010.010.010.010.010.010.010.010.010.010.010.00.050.050.050.050.050.050.010.010.010.010.010.00.050.050.050.050.050.050.010.010.010.010.010.00.050.050.050.050.050.050.010.010.010.010.010.00.050.050.050.050.050.050.010.010.010.010.010.00.050.050.050.050.050.050.010.010.010.010.010.00.050.050.050.050.050.050.010.010.010.010.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.050.050.010.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.010.050.01;b=200 2752152902903652503252653403404156507256657407408152303052453203203952803552953703704456807556957707708453804553954704705454305054455205205958309058459209209952303052453203203952803552953703704456807556957707708452603352753503504253103853254004004757107857258008008754104854255005005754605354755505506258609358759509501025;for i=1:108 a=xi(:,i);y=105/9*(Yx(1)*a(1)*x(1)2+(Yx(2)*a(2)*x(2)2+(Yx(3)*a(3)*x(3)2+(Yx(4)*a(4)*x(4)2+(Yx(5)*a(5)*x(5)2+(Yx(6)*a(6)*x(6)2+(Yx(7)*a(7)*x(7)2)+105*(174.42*(x(1)/x(5)*(x(3)/(x(2)-x(1)0.85*(1-2.62*(1-0.36*(x(4)/x(2)(-0.56)1.5*(x(4)/x(2)1.16)/(x(6)*x(7)0.5-1.5)2+b(i);c(i)=y;endf,j=min(c);Main.maeq=;beq=;vlb = 0.0750 0.2250 0.0750 0.0750 1.1250 12.0000 0.5625;vub = 0.1250 0.3750 0.1250 0.1250 1.8750 20.0000 0.9350 ;x,fval=fmincon(fun,x0,A,b,aeq,beq,vlb,vub)七.储油罐的变位识别与罐容表标定 1. 问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。 许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向