旋转抛物面、市场份额分配.docVIP

实验一 实验目的与要求 集合体是由旋转抛物面、抛物柱面与平面y=2所围成且在第I卦限的部分，绘出该立体的图形及其在xOy平面上的投影区域，并计算该立体的体积。 实验方案 程序源文件如下： function y=Ex001 clf,clear syms x y z; int(int(int(1,z,0,sqrt(y-x.^2)),x,0,sqrt(y)/2),y,0,2) %求体积 clear dd=0.01; n=1+2/dd; [x y]=meshgrid(0:dd:2,0:dd:2); for j=1:n for i=1:n if x(i,j)sqrt(y(i,j)./4) y(i,j)=2 z(i,j)=sqrt(y(i,j)-x(i,j).^2); else z(i,j)=0; end end end mesh(x,y,z); %绘制立体图像 view([5 1 8]); %确定观察角 hold on x=0:0.001:(1/sqrt(2)); %绘制投影曲线 y1=4.*x.^2; y2=2; plot(x,y1,b,x,y2,b); 实验结果和数据处理 运行程序得如下立体图形： 运行得到给立体体积为： 实验二 实验目的与要求 有两家公司R和S经营同类的产品，它们相互竞争。每年R公司保有1/4的顾客，而3/4转移向S公司；每年S公司保有2/3的顾客，而1/3转移向R公司。当产品开始制造时R公司占有3/5的市场份额，而S公司占有2/5的市场份额。问两年后，两家公司所占的市场份额变化怎样，五年后会怎样？十年以后如何？是否有一组初始市场份额分配数据使以后每年的市场分配稳定不变？ 实验方案 依题意可得本题数学模型为 其中，记，，x为R公司的市场份额，y为S公司的市场份额。可得程序源文件为： function y=Ex002 A=[1/4 1/3;3/4 2/3]; [v,d]=eig(A) %得到矩阵A的特征值和特征向量 syms n; %定义符号变量n X0=[3/5;2/5] %初始值 X=v*d.^n*inv(v)*X0 %第n年的结果 for i=1:10 %从第一年到第十年分别代入n值 i subs(X,i) end 实验结果和数据处理 运行上述程序，得出以下结果。由图2知，两年后，R公司占有0.3097的市场份额，S公司占有0.6903的市场份额。 图1 图2 由图4可看出，五年后R公司占有0.3077的市场份额，S公司占有0.6923 的市场份额。 图3 图4 由图6可看出，十年后R公司占有0.3077的市场份额，S公司占有0.6923 的市场份额。 图5 图6