MatLb实验报告.docxVIP

MatLab实验报告学院：__________________班级：__________________姓名：__________________学号：__________________一、实验内容实验一（油库问题）1.问题描述2.问题分析由于油库的横截面为圆形，故此问题可转化为三重积分问题。由于纵截面两边的抛物线方程可求得，故采用“先重后单”的积分方法。为此，先求得抛物线方程，即可得某高度对应的横截面面积，从而通过对高度的积分得出所求体积。为了便于得出表达式，采用符号函数积分的方法。通过eval得出所需V的表达式。对于标尺刻度的标定，由于人们习惯于整数刻度，故对每一整数刻度，遍历高度值寻找与之最接近的对应体积，在这个体积对应的高度处标定刻度。3.MatLab程序clear;clc;clf;syms h r S fV a b c rhrm H htrh=2; %上下底上的半径rm=3; %中心位置上的半径H=12; %高度%-----------------------------------计算库容V的表达式a=(rh-rm)/((H/2)^2); %根据油库截面参数计算侧面抛物线参数x=a*(y-b)^2+cb=H/2;c=rm;r=sym(a*(ht-b)^2+c); %高度h处的横截面圆半径S=sym(pi*(a*(ht-b)^2+c)^2); %高度h处的横截面圆面积fV=int(S,ht,0,h); %在高度h处取一高度为dh的圆形薄片，则薄片%体积dV=S*dh，总体积为dV从0到h的积分， %即int(S,ht,0,h)disp(油库内油面的深度为h时)disp(库内油量的容积V的表达式为：V(h)=)disp(eval(fV)) %在命令窗口显示V(h)的计算表达式%-----------------------------------计算并绘制库容标尺eps=0.1; %标定容量的精确度（单位：m）div=20; %标尺分度值（单位：m^3）axis equalaxis([-4,8,0,12]);hold ongridht=0:0.1:12;xl=-eval(r); %根据h产生h处的半径序列xr=eval(r);title(油库横截面示意图)xlabel(半径/m)ylabel(深度/m)plot(xl,ht,r-,xr,ht,r-,linewidth,3); %绘制油库横截面示意图plot([2,4],[0,0],b-,linewidth,2); %绘制V=0处的标尺text(4,0,[num2str(0)]);h=12;plot([2,4],[12,12],b-,linewidth,2); %绘制最大容量处的标尺text(4,12,[num2str(floor(eval(fV)))]);for i=div:div:floor(eval(fV)) %从20到最大容量绘制标尺 for j=0:eps:12.1 %搜索对应某容量的高度值 h=j;Vl=eval(fV); h=j+eps;Vr=eval(fV); if (Vli Vri) %如果i落在[V(j),V(j+eps)]内，则在高度j处标定容量i break end endht=j; plot([eval(r),4],[j,j],b-,linewidth,2); %标定容量 text(4,j,[num2str(i)]);endtext(5,11,[油量/m^3]) %显示说明文字4.运行结果V(h)的表达式油量刻度的标定实验二（按揭贷款计算）1.问题描述2.问题分析按题意，所求月还款额应为等额本息还款方式。为此，先推导计算月还款额度Rm的表达式：设需还款额为A，第n个月还款前剩余贷款为an，第n个月应偿还利息为In，月利率为im，则：此公式即为进行MatLab编程时所需的原理式。3.MatLab程序clear;clc;s=120; %建筑面积prc=4200; %单价每m^2arate=5.5/100; %年利率years=20; %还款年数total=s*prc; %房屋总价格