Mathematica实验六 线性方程组.pptVIP

实验六 线性方程组 实验目的：理解齐次线性方程组解的几何意义，了解线性方程组的一些应用 实验1 求解方程组 ，并画出 的图形 实验六 线性方程组 Needs[Graphics`Colors`] arrow[{a_,b_},color_]:=Graphics[{color,Line[{{0,0}, {a,b}}],Line[{{a,b},0.9{a,b}+0.04{-b,a}}], Line[{{a,b},.9{a,b}-.04{-b,a}}]}] A={{-1,3},{4,-12}};RowReduce[A];x={1,-3}; NullSpace[A] y={3,1}; Show[arrow[x,Red],arrow[y,Green],AspectRatio-Automatic,Axes-True]; 实验六 线性方程组 练习1 求解线性方程 ，并画出 与X的图形。 观察图形并给出你的结论 实验六 线性方程组 实验2 求出通过平面上三点(0,7),(1,6)和(2,9)的二次多项式ax2+bx+c，并画出其图形 根据条件有 0*a+0*b+c=7 1*a+1*b+c=6 4*a+2*b+c=9 实验六 线性方程组 Clear[x]; A={{0,0,1},{1,1,1},{4,2,1}};y={7,6,9} p=LinearSolve[A,y] Clear[a,b,c,r,s,t]; {a,b,c}.{r,s,t} f[x_]=p.{x^2,x,1}; Plot[f[x],{x,0,2}]; Plot[f[x],{x,0,2},GridLines-Automatic,PlotRange-All] 实验六 线性方程组 练习2 求出通过平面上四点(-2,6),(1,4),(2,3),(3,2)的3次多项式ax3+bx2+cx+d，并画出其图形 实验六 线性方程组 实验3 求出通过平面上三点(0,0),(1,1),(-1,3)以及满足f(-1)=20,f(1)=9的4次多项式f(x),并画出其图形 分析：确定一个4次多项式需要知道它的5个系数。通过3个插值条件和两个导数条件可以建立5个线性方程 实验六 线性方程组 设f(x)= ax4+bx3+cx2+dx+e，则有 e=0 a+b+c+d+e=1 a-b+c-d+e=3 -4a+3b-2c+d=20 4a+3b+2c+d=9 实验六 线性方程组 LinearAlgebra`MatrixManipulation` Clear[a,b,c,d,e]; q[x_]=a*x^4+b*x^3+c*x^2+d*x+e; eqs={q[0]==0,q[1]==1,q[-1]==3,q[-1]==2,q[1]==9}; {A,y}=LinearEquationsToMatrices[eqs,{a,b,c,d,e}]; p=LinearSolve[A,y]; f[x_]=p.{x^4,x^3,x^2,x,1}; Plot[f[x],{x,-1,1},GridLines-Automatic,PlotRange-All] 最小二乘法 实验5 给定数据点(0.7,4.0),(3.3,4.7),(5.6,4.0), (7.1,1.3),(6.4,-1.1),(4.4,-3.0),(0.3,-2.5),(-1.1,1.3).试找一个圆，这些数据点尽可能落在该圆上 xs={0.7,3.3,5.6,7.1,6.4,4.4,0.3,-1.1} ys={4.0,4.7,4.0,1.3,-1.1,-3.0,-2.5,1.3} pts=Transpose[{xs,ys}]; Show[Graphics[{Red,AbsolutePointSize[5],Map[Point,pts],Violet,Circle[{3,1},3.5]}],Axes-Automatic,AspectRatio-Automatic]; Transpose[A]：求矩阵A的转置阵 Map[f,expr]映射将f分别作用到expr第一层的每一个元上得到的列表 Map[f,{a,b,c,d,e}] {f[a],f[b],f[c],f[d],f[e]} Point[{coords1,coords2,}] represents a collection of points. Graphics[Point[Table[{t,Sin[t]},{t,0,2p,2p/10}]]] violet？ 最小二乘法 设圆的中心在(c1,c2)，半径为r，则圆的方程为