北师MATLAB课件实验课二.pptVIP

数学模型实验课（二） 符号运算、矩阵 问题1.讨论收益模型 在售猪问题中， P(t) = (0.65-rt)(200+gt)-kt 1.对每天的饲养花费做灵敏度分析，分别考虑对最佳售猪时间和相应收益的影响。 2.如果有新的饲养方式，每天的饲养花费为60美分，会使猪按7磅/天增重。那么是否值得改变饲养方式？ 3.求出使饲养方式值得改变的最小的增重率。 一、符号运算 syms x %定义符号变量 sin(x^3)+3*sin(x^3) %函数代数运算 diff(sin(x^2)) %函数求导 diff(sin(x^2),3) %求3阶导数 int(-2*x/(1+x^2)^2) %函数积分 int(x/(1+z^2),z) %对变量z的含参量积分 int(x*log(1+x),0,1) %定积分 int(2*x, sin(t), 1) %变上限积分 x=1:3; subs(x)%带入算数 句柄函数的优点： g = @(c) (quad(@(x) (x.^2 + c.*sin(x)), 0, 1)); g(2) syms t k f=(0.65-0.01*t)*(200+5*t)-k*t diff(f,t) ans = 5/4-1/10*t-k t0=(5/4-k)*10 % k对最佳售猪时间的影响。 compose(f,t0,t,k)% k对净收益f的影响 ans = (21/40+1/10*k)*(525/2-50*k)-k*(25/2-10*k) 问题2 在量纲分析中，求解齐次线性方程组 x1 – 2x4 = 0 x2 + x4 = 0 x3 = 0 x5 任意 二、矩阵 1. 矩阵的生成 a=[1 2 3 -4; 5 6 7 8; -1 3 -5 7] %生成矩阵 b=[a; 2 4 6 8] % 添加一行 c=[b,b(4,:)’] % 将b的第四行转置为一列添在后面 d=ones(2) %2×2的全 1 矩阵 f=ones(3;1) %全1向量 I=eye(3) %3×3的单位矩阵 v=repmat(0,[3,5,7]) ％三维矩阵 C=cat(3,zeros(3,5),ones(3,5),randn(3,5),rand(3,5)) 2. 矩阵的运算 h=A’ %矩阵转置 A±B %矩阵（与数）加减 A.*B , A./B %矩阵对应元素相乘，除 X=A*B %矩阵相乘 X=A.^2 %矩阵各元素乘方 X=inv(A) %方阵A求逆 X=A\B %左除，A*X=B, X=inv(A)*B, 解方程 X=A/B %右除，X*B=A, X=A*inv(B) size(A) %矩阵A的维数 length(a) %向量 a 元素的个数 3 矩阵函数 sin(a)+2.*a b.^2 b^2 det(b) rank(b) a1=sum(a) %矩阵 a 各列元素的和 sum(a1) %向量 a1 各分量的和 a1=max(a) %矩阵 a 各列元素的最大元 max(a1) %向量 a1 各分量的的最大元 3.矩阵的特征值，特征向量 AX=λX eig(A); %返回值就是特征值矩阵和特征向量矩阵 解问题1 A=[1 0 0 -2 0;0 1 0 1 0;0 0 1 0 0 ]; C=zeros(2,5); B=[A; C]%B是方阵 [v,d]=eig(B) % d是 5阶对角方阵，对角线上元素为B的特征值，v是 5阶方阵，其列向量为相应的特征向量。 v = 1.0000 0 0 0.8165 0 0 1.0000 0 -0.4082 0 0 0 1.0000 0 0 0 0 0 0.4082 0 0 0 0 0 1.0000 d = 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0