matlab数学问题实例,matlab数学实验实例

matlab数学实验实例

数学实验与 Matlab 11《数学实验与 Matlab》程序周晓阳华中科技大学数学系我将程序按书中的顺序排列，这样便于读者利用。 本书程序均通过了调式。可直接拷贝到命令窗口运行或编制 M 文件运行。如出现问题，可能是中英文标点的缘故(排版有可能使用了中文标点) ，请将中文标点换为英文标点试试。数学实验与 Matlab 22实验１．矩阵运算与 Matlab 命令1.1 知识要点与背景：日常矩阵及其运算【 A=[4 2 3;1 3 2;1 3 3;3 2 2], % 表 1-1、 表 1-2的 数 据 分 别 写成 矩 阵 形 式B=[35 20 60 45;10 15 50 40;20 12 45 20] 】【 C=A*B %矩 阵 乘 法 ， 求 各 订 单 所 对 应 的 原 材 料 和 劳 动 力 。 】 【 whos % 查 看 Matlab工 作 空 间 中 变 量 及 其 规 模 】 1.2 实验与观察：矩阵和 Matlab 语言1.2.1 向量的生成和运算【 x=linspace(0,4*pi,100); %将 [0,4π ]区 间 100等 分 ， 产 生 了 一 个100维 向 量y=sin(x); %计 算 函 数 值 ， 产 生 了 一 个 与 x同 维 的 100维 函 数 向 量 y y1=sin(x).^2; %计 算 函 数 向 量 ， 注 意 元 素 群 运 算y2=exp(-x).*sin(x); %以 x为 横 坐 标 ， y为 纵 坐 标 画 函 数 的 图 用 不 同 的 线 型 将 函 数 曲 线 绘 制 在 一 个 图上plot(x,y, - ,x,y1, - ,x,y2, .- ) 】 1. 向量的创建◆直接输入向量。 【 x1=[1 2 4],x2=[1,2,1],x3=x1 】◆冒号创建向量 。【 x1=3.4:6.7x2=3.4:2:6.7数学实验与 Matlab 33x3=2.6:-0.8:0 】◆生成线性等分向量。 【 x=linspace(0,1,5) 】2. 向量的运算【 y=sin(x) 】【 y1=sin(x).^2; y2=exp(-x).*sin(x); 】 1.2.2.矩阵创建和运算1.创建矩阵(1)数值矩阵的创建◆直接输入法创建简单矩阵。【 A=[1 2 3 4; 5 6 7 8; 9 10 11 12] 】 【 B=[-1.3,sqrt(3);(1+2)*4/5,sin(5);exp(2),6] 】 (2)符号矩阵的创建◆【 syms a11 a12 a13 a14 a21 a22 a23 a24 a31 a32 a33 a34 …b11 b12 b13 b14 b21 b22 b23 b24 b31 b32 b33 b34A1=[a11 a12 a13 a14 ;a21 a22 a23 a24; a31 a32 a33 a34],B1=[b11 b12 b13 b14 ;b21 b22 b23 b24; b31 b32 b33 b34] 】 2.矩阵的运算【 C=A1+B1,D=A1-B1 】【 syms ccA=c*A1 】数学实验与 Matlab 44【 C=A1*B1 】{ ??? Error using ==> sym/mtimes, Inner matrix dimensions must agree. }【 A2=A1(:,1:3), B1 】【 G=A2*B1 】【 g11=A2(1,:)*B1(:,1) 】【 A, A_trans=A 】【 H=[1 2 3 ; 2 1 0 ; 1 2 3 ], K=[1 2 3 ; 2 1 0 ; 2 3 1]h_det=det(H), k_det=det(K),H_inv=inv(H),K_inv=K^-1 】【 A=[3 0 1; 1 1 0;0 1 4];B=inv(A-2*eye(3))*A, B=(A-2*eye(3))\A 】3.分块矩阵：矩阵的裁剪、分割、修改与抽取(1)【 A=[1 0 1 1 2;0 1 -1 2 3;3 0 5 1 0;2 3 1 2 1], vr=[1,3];vc=[1,3];A1=A(vr,vc) %取 出 A的 1、 3行 和 1、 3列 的 交 叉 处 元 素 构 成 新 矩 阵A1 】◆将上面的矩阵A分为 四块，并把它们赋值到矩阵B中， 观察运行后的结果。【 A11=A(1:2,1:2),A12=A(1:2,3:5),A21=A(3:4,1:2),A22=A(3:4,3:5)B=[A11 A12;A21 A22] 】(2)矩阵的修改和提取◆ 【 A=[1 0 1 1 2;0 1 -1 2 3;3 0 5 1 0;2 3 1 2 1]A(1,:)=[0 0 0 0 0]; A 】数学实验与 Matlab 55◆ 观察：【 B(:,[2,4])=[ ] %删 除 矩 阵 B的 第 2、 4列 】(3)矩阵元素的抽取4.生成特殊矩阵。◆【 y1=rand(1,5), y2=rand(1,5),rand( seed ,3), x1=rand(1,5), rand( seed ,3), x2=rand(1,5) 】5. 常用矩阵函数6. 数据的简单分析◆【 rand( seed ,1);A=rand(3,6),Asort=sort(A), Amax=max(A), Asum=sum(A) 】1.2.3 Matlab 工作环境和编程2.Matlab 的基本设计1.3 应用、思考与练习1.3.1 关系矩阵1.3.2 投入产出 1.3.3 循环比赛的名次【 A=[0 1 1 0; 0 0 1 1; 0 0 0 1; 1 0 0 0],数学实验与 Matlab 66e=ones(4,1); c=A*e; s=c 】★ 画矩阵结构图的g plot指令。◆(3) 【 clf, A=[0 1 1 0;0 0 1 1;0 0 0 1;1 0 0 0]; xy=[0 1;0 0;-1 –0.5;1 –0.5];graphy_plot(A,xy,1,0.5), % gplot(A,xy) 】1.3.4 参考程序graphy_plot.m【 function y=graphy_plo