1、 复习matlab的基本编程知识，画图，插值，拟合，方程求解。阅读参考资料。
2、 已知速度曲线v(t) 上的四个数据点下表所示。
思路：利用插值（即求过已知有限个数据点的近似函数）的基本原理，用多项式作为研究插值的工具，进行代数插值。
**：clc
t=[0.15 0.16 0.17 0.18];
v=[3.5 1.5 2.5 2.8];
tt=0.15:0.00001:0.18; %细分已知点。
v1=spline(t,v,tt三次样条插值。
plot(t,v,'r*',tt,v1,'g');画出图像。
s= sum(v1.*(0.18-0.15)/3000)%利用积分定义求位移。
a = 2.8-v1(2999))/0.18-0.15)/3000)%加速度为速度的导数，利用倒数的定义求0.18处的加速度。
结果：s = 0.0687
a= -255.7872
图形：3.计算**文件 给出的两个圆 a，b 的圆心，和两个圆的两条外公切线和两条。
内公切线的切点的坐标。
**：data=imread('c:\documents and settings\administrator\桌面\')
data2=double(data(:,1));
image(data2找出圆的位置；
x,y]=find(data2==4);
hold on设置分界点。
yc=mean(y)；
data3=data2(:,1:221);[x1,y1]=find(data3==4);
x1c=mean(x1);y1c=mean(y1);plot(y1c,x1c,'rp');
data4=data2(:,221:512);[x2,y2]=find(data4==4);
x2c=mean(x2);y2c=mean(y2)+221;plot(y2c,x2c,'rp');
r1=x1c-min(x1)
r2=x2c-min(x2求半径。
m=solve('(a*x1c+y1c+b)^2/(1+a^2)=r1^2','a*x2c+y2c+b)^2/(1+a^2)=r2^2','a,b')
a=eval(
xx=0:0.01:350;
yy1=-a(1)*xx-b(1);yy2=-a(2)*xx-b(2);yy3=-a(3)*xx-b(3);yy4=-a(4)*xx-b(4);
plot(yy1,xx,'r')
plot(yy2,xx,'b')
plot(yy3,xx,'g')
plot(yy4,xx,'y求切点。
for i=1:4
f1=solve('(x-x1c)^2+(y-y1c)^2-r1^2','y+a(i)*x+b(i)',x','y');
q2x=eval('m*')
endfor i=1:4
f1=solve('(x-x2c)^2+(y-y2c)^2-r1^2','y+a(i)*x+b(i)',x','y');
q2x=eval('m*')
end结果：
a的半径：r1 =79.7500；
b的半径：r2 = 80；
切线方程：内公切线：y=-0.0031x+166.5626；
y= 2.8972x-475.3679；
外公切线：y= 0.7090x -88.9699；
y= 0.7118x+ 106.3791；
a的切点：（109．81，166.22）、（184.95，60.45）、（155.68 ,21.40）、（63.31，151.44）；
b的切点：（333.21，165.52）、（257.84,271.59）、（379.73,180.22）、（287.08,310.67）。
图形：4.求微分方程组的数值解，并画出解曲线。
程序**：建立m文件。
function dy=w(t,y)
dy=zeros(3,1);
dy(1)=-10*y(1)+10*y(2);
dy(2)=28*y(1)-y(2)-y(1)*y(3);
dy(3)=-8/3)*y(3)+y(1)*y(2);
取t0=0,tf=1,输入命令[t,y]=ode45('z',[0 1],[1 1 0]);plot(t,y(:,1),'o',t,y(:,2),'t,y(:,3),'d');
数值解：y1 = y2y3=
数学建模作业
1992年全国数学建模a题数据分析。对土豆和生菜，分别建立施肥量和产量之间的多元关系，运用excel和matlab软件依次采用散点法和拟合关系进行绘图显示。在确认模型具有完美的基础上，进行线性相关 交互作用 最佳响应水平 强影响变量等的分析。同时，将两种作物进行比较，得出一系列颇有实用价值的结论。分...
数学建模作业
数学建模作业 第二章。微分方程稳定性分析。1 由解得平衡点，特征值，是不稳定点。matlab作图程序 x 5 0.1 5 fora 5 0.1 5 y a x plot x,y hold on endxlabel x ylabel y 得到图形 2 由解得平衡点，特征值，是不稳定点。matlab作图...
数学建模作业
1.试说明下面的数学规划 mp 是一个凸规划。解 1 目标函数。一阶顺序主子式为，行列式为 二阶顺序主子式为，行列式为 三阶顺序主子式为，行列式为。hessian矩阵正定，目标函数是凸函数。2 约束条件。一阶顺序主子式为，行列式为 二阶顺序主子式为，行列式为 三阶顺序主子式为，行列式为。hessia...