matlab用梯形法计算河床断面积,Matlab课计算题.doc

Matlab课计算题

PAGE \* MERGEFORMAT PAGE \* MERGEFORMAT 1Matlab计算题：1，农夫老李有一个半径为10m的圆形牛栏，里面长满了草，老李要将家里的一头牛拴在牛栏边的一根栏桩上，要求只让牛吃到圆形牛栏中的一半的草，请问栓牛鼻的绳子应为多长？答案为11.5873米程序如下。clear,clc,close allsyms R;cos_r=(200-R^2)/200;cos_R=R/20;afa_r=2*acos(cos_r);afa_R=2*acos(cos_R);f=50*afa_r+R^2*afa_R/2-50*sin(afa_r)-R^2*sin(afa_R)/2-pi*100/2;R=solve(f);disp(['栓HYPERLINK "/search?word=牛鼻子&fr=qb_search_exp&ie=utf8"牛鼻子的绳长应为：',HYPERLINK "/search?word=num2str&fr=qb_search_exp&ie=utf8"num2str(eval(R)),'米。'])syms x yf1=x^2+y^2-100;f2=(x-10)^2+y^2-R^2;HYPERLINK "/search?word=ezplot&fr=qb_search_exp&ie=utf8"ezplot(f2 ,[-12,12])axis([-12,12,-12,12])axis squarehold onplot([0 10],[0 0],'*b')HYPERLINK "/search?word=ezplot&fr=qb_search_exp&ie=utf8"ezplot(f1,[-10,10])hold off2，图5.4所示,:为了在海岛I与某城市C之间铺设一条地下光缆,每千米光缆铺设成本在水下部分是C1,在地下部分是C2,为使得铺设该光缆的总成本最低,问题1 求总造价最小的P点坐标f=inline('3000*x/(225+x^2)^(1/2)-1500*(30-x)/((30-x)^2+100)^(1/2)')%对总造价函数的导函数通过二分法求零点，来求最值a=0;b=30;dlt=1.0e-3; %根据题意，误差小于10k=1;while abs(b-a)>dltc=(a+b)/2;if f(c)==0break;elseif f(c)*f(b)<0a=c;elseb=c;endfprintf('k=%d,x=%.3f\n',k,c);k=k+1;end X=c3，有一艘宽为5m的长方形驳船欲驶过某河道的直角弯，经测量知河道的宽度为10m和12m，试问，要驶过该直角湾，驳船的长度不能超过多少米？(误差<0.001m)驳船的长度问题：　g=inline(12/cos(x)+10/sin(x)-5/[sin(x)*cos(x)]) f=inline(12/cos(x)^2*sin(x)-10/sin(x)^2*cos(x)+5/sin(x)^2-5/cos(x)^2) a=0+eps;　　b=pi/2-eps;　　dlt=1.0e-3;　　k=1;　　while abs(g(b)-g(a))>dlt　　c=(a+b)/2;　　if f(c)==0　　break;　　elseif f(c)*f(b)<0　　a=c;　　else　　b=c;　　end　　fprintf('k=%d,x=%.3f\n',k,c);　　k=k+1;　　end　　g(c)P811，某农夫有一个半径10米的圆形牛栏，长满了草，他要将一头牛拴在牛栏的边界的栏桩上，但只让牛吃到一半的草，问栓牛鼻的绳子应为多长？解：clear;syms a t x;y=pi/2*a*a-pi*a*x-int('sqrt(a*a+2*a*(t-x)-t*t)-sqrt(a*a-t*t)','t',x,a);%y=subs(y,a,10);f=inline(char(y));x=fzero(char(y),1)L=sqrt((110-x)^2+10^2-x^2)2.如图所示，为了在海岛I与某城市C之间铺设一条地下光缆，每千米光缆铺设成本在水下部分是C1，在地下的部分是C2，为使的铺设该光缆的总成本最低，光缆C1的转折点P(海岸线上)应取在何处？ 如果实际测得海岛I与城市C之间水平距离L=30km，海岛距海岸垂直距离h1=15km，城市距海岸线垂直距离h=10km，C1=3000万元/km，C2=1500万元/km，求p点坐标(误差《10-3km)解：syms h1 h2 l c1 c2 x;y=c1*sqrt(h1^2+x.^2)+c2*sqrt(h2^2+(l-x)^2);dy=diff(y,x)f=subs(dy,{h1,h2,l