高等数学Matlab解法幻灯片.pptVIP

　　3. 解下列方程组。 　　(a) 　　　　　　　　　　　(b) 　　 　　4. 设 　　　　　　　　　　　　　　　　　　　　　， 求能满足Ax=0的全部的x。 　　5. 设向量［v1，v2，v3］张成的空间为H，判断向量x1和x2是否在这个空间内，其中 　　6. 设下列向量张成一个子空间，试找出哪些向量是该子空间的基向量。 　　(a) 　　(b) 　　7. 根据开普勒第一定律，当只考虑一个大天体的重力吸引时，一个小天体应该取椭圆、抛物线或双曲线轨道。在适当的极坐标中，小天体的位置(r，θ)应满足下列方程: r=β+e(r·cosθ) 其中β为常数，e是轨道的离心率。对于椭圆， 0≤e＜1； 对于抛物线，e=1; 对于双曲线，e＞1。对一个新发现的天体的观测得到了如下数据: 试确定其轨道的性质，并预测此天体在θ=0.46弧度时的位置。 　　8. 求通过(-1，0)，(0，2)，(2，3)，(2，-2)，(0，-3)五点的二次圆锥截面曲线。设其方程为 ax2+bxy+cy2+dx+ey+1=0 求其中的系数a、b、c、d、e。如果还要它通过点(1，2)，该如何处理? 试作讨论。 　　用MATLAB的pinv函数可以求出此超定方程的最小二乘解,键入 a1=pinv(A1)*b1 则 p1(t)=3.1714-3.3571t-0.8214t2+0.7500t3 要画出它的图形，可以继续键入 t=-1∶0.1∶4; plot(t，a1(1)+a1(2)*t+a1(3)*t.^2+a1(4)*t.^3，′: r′) plot(-1，5，′x′) 　　这条曲线也画在图5-42上，用实线表示。可见它虽然不通过给定的5个点中的任何一个，但比原来的曲线靠近点(-1，5)很多。各个方程的具体误差向量E及五项误差的平方和可以用矩阵方程计算如下: E=A1*a1-b1，EE=norm(E) 求得结果为 说明五个方程都有误差，但误差都已减少到很小，其平方和为最小。 　　【例5-5-5】 用有理整数来配平下列小苏打与柠檬 酸反应的化学方程。 NaHCO3+H3C6H5O8→Na3C6O7+H2O+CO2 　　解: 建模 　　将五种物质的量作为五个未知数，每种物质含有的元素Na、H、C、O写成列向量。 然后按四种元素左右平衡列出四个方程，将各变量项移至方程左端，得  其中 ? 　　现在的问题就归结为解出x向量，要求其各个分量都以最小整数表示。 　MATLAB程序 A=［1，0，-3，0，0; 1，8，0，-2，0; 1，6，-6，0，-1; 3，8，-7，-1，-2］, b=［0; 0; 0; 0］ U=rref(［A，b］) 　　程序运行结果 　　程序运行的结果为 　　在化学方程中，不允许出现小数，用手工来凑整数不大容易，因此我们将MATLAB命令窗设置为有理整数的显示格式，再运行这个程序，即键入 format rat， U=rref(A) 得到 　　将系数矩阵复原为其同解方程，可得 ? ? 由此可见，x5最小必须取30，才能保证各个分量xk为整数，于是得到 x1=48，x2=13，x3=16，x4=76，x5=30 最后可以得出配平后的化学反应方程为 48NaHCO3+13H3C6H5O8→16Na3C6O7+76H2O+30CO2 　　这个待配平的化学方程有四个方程和五个未知数，所以它是一个欠定的线性方程组，其解乘以任何常数仍然为方程的解。不过当我们另加了一个条件——系数取最小的正整数，结果就是唯一的了。 5.5.2 向量相关性 　　【例5-5-6】 对于由三个在R4空间的四维基向量v1、v2、v3张成的子空间，问w1和w2是否在此子空间内，其中 　　解: 本题的要点在于研究w是否能由v1、v2、v3以线性组合的方式组成，即是否能找到三个常数c1、c2、c3，以便得到三个基向量的线性组合: c1v1+c2v2+c3v3=w 能实现这个关系的w与向量组v1、v2、v3线性相关; 反之就线性无关。 　　(1) 方法一 利用向量组的最大无关组的概念，可以解决这个问题。把五个向量列成一个矩阵A，对它作行阶梯变换，从最简行阶梯的构成来判断其相关性。键入 A=［v1，v2，v3，w1，w2］， U0=rref(A) 得到 可见最大无关组是由v1、v2、v3、w1四个向量组成的，也就是说这四个向量线性无关，w1不可能由v1、v2、v3组合而成。而w2与v1、v2、v3构成的矩阵却只有三行，其秩为3，所以它与向量组v1、v2、v3线性相关，即可以由v1、v2、v3的