数学建模 计算机程序 MATLAB

主要讲下思路

s=at+b 直接有 函数使用，计算结果。樱郑侍s是y的函数，t是x的函数。

1)

a) y = [a b c]*[x^2 x 1]^T 矩阵形式

令 t0=x^2 t1=x

s=y-1 有 s =at0+bt1 这是丛祥一个关于t0 和t1的最小二乘

b) lny = ln(a)+n ln(x)

t=ln(x) s=ln(y) b=ln(a) 有 s=b+nt 线性形式

注意 如果a小于0,需要做个处理 对所有y取相反数 然后按上面形式去做

2)

a) ln(y) = ln(b)+nln(x) y、b、x均脊吵为正数

s=ln(y) t=ln(x) c=ln(b)

s=c+nt

b) ln(y) = ln(b)+cx

s=ln(y) a=ln(b) t=x

s=a+ct

c) y = a ln(x)+b

s=y t=ln(x)

s=at+b

d) ln(y)=ln(a)+2x

s=ln(y) b=ln(a) t=x

s=b+2t

e) ln(y)=ln(a)+3x

s=ln(y) b=ln(a) t=x

s=b+3t

计算s和t的相关系数后，就可以推出原来的y x形式的系数

题主给出的算例是斐波那契数列计算。其中题主框出的代码的意思是这样的：

X=[] %对变量初始化，清零

for k=1:7 %循册仿或环语句计州伍算每列斐波那契数列

X=[X,myfib(k)]%调用斐波那契数列计算函数，并储存在X向量数组里

end %结束循环

X %显示斐波那契数列

如x=1，2，3……7分别带入myfib(k)函数，是可以得到一个斐波大凯那契数，而不是一个数列

运行效果如图所示。

----------------

function dx=appollo(t,x)

mu=1/82.45

mustar=1-mu

r1=sqrt((x(1)+mu)^2+x(3)^2)

r2=sqrt((x(1)-mustar)^2+x(3)^2)

dx=[x(2)

2*x(4)+x(1)-mustar*(x(1)+mu)/r1^3-mu*(x(1)-mustar)/r2^3

x(4)

-2*x(2)+x(3)-mustar*x(3)/r1^3-mu*x(3)/r2^3]

--------------------------

x0=[1.200-1.04935751]%x0(i)对应与xi的初值

options=odeset('reltol',1e-8)%该命令的另一种写法早耐是options=odesetoptions.reltol=1e-8

[t,y]=ode45(@appollo,[0,20],x0,options)%t是时间点，y的第i列对应xi的值，t和凯汪y的行数相同

plot(y(:,1),y(:,3))%绘制x1和盯睁仔x3，也就是x和y的图形

title('Appollo卫星运动轨迹')

xlabel('X')

ylabel('Y')

---