MATLAB程序设计大赛培训讲座文稿.docVIP

MATLAB程序设计大赛培训讲座文稿 ——李金武 前言： 在讲座开始前，给大家先介绍一下本学期讲座的安排。这学期的MATLAB培训是为MATLAB程序设计大赛做准备的，总共有四次，经过四次后，相信你们在软件使用的能力上都会有所提高。希望每次通知你们参加时，你们都能积极的参加。在MATLAB程序设计大赛中获奖的话，可以直接参加暑期培训的哦。培训大体是这样安排的：1. MATLAB的基本操作以及在数学方面的应用（主要是高数和线性代数）；2.数据可视化（其实就是作图）；3.结构化程序设计；4.多项式运算；5.函数的编写；6.真题回顾，主要是两年的程序设计的初赛题目。7.如果你们有兴趣的话，还可以给你们介绍一下一种高级算法，神经网络。 一、科学计算器功能 exp()：以e为底的幂函数； sin：求正弦值(弧度制） sind：求正弦值(角度制） asin：反正弦(弧度制) asind：反正弦(角度制) ^：幂次方 pi：。 ＋，加；－，减；＊，乘；／，除； log以e为底的自然对数,其他的对数形式用换底公式。 如：要计算的值，可以转换为求的值。 在MATLAB中输入：log(9)/log(3)，结果ans＝2 。 下面来一个综合的例子： :求的值， f=(sqrt(sin(6*pi/7)-0.001)+(log(4)^6))/(exp(2/3)-cos(1000*pi)); f= 8.1836 插入语：easy；so easy。 二、在线性代数中的应用 2.1、矩阵的输入方式 数组和矩阵是MATLAB中最基本的，所有的问题最终都将转化为矩阵的形式来求解，下面先来讲一下矩阵的输入及其运算。一般的矩阵基本上有以下三种输入的方式： 方式一： a=[1 2 3;4 5 6; 7 8 9]; 方式二： a=[1,2,3;4,5,6;7,8,9]; 方式三： a=[1 2 3 4 5 6 7 8 9]; 等比数列的产生 方法一：A=[begin:step:end] begine起始点，step步长,默认为1，end结尾点。 A=[0:2:100]; %即可产生100内的偶数 方法二：A=linspace(begine,end,number) begine起始点，end结尾点，number为点的个数，默认为100。 A=[1,100]；%将1到100平分为100段 2.2、矩阵，数组的连接 水平连接： c=[a,b]; 竖直连接：c=[a;b]; 必须注意的是：连接的前提是要满足相应维数相等的。 2.3、计算数组大小 [m,n]=size(A)：m为数组A的行数，n为数组A的列数； numel(A)：返回数组A的元素个数； a=length(A); 当A是一维数组时，a为数组A的元素个数 当A是普通二维数组时，a为m,n中较大的一个数。 2.4、矩阵的逆运算，行列式，转置，秩 逆运算：inv(A)或A^(-1); 行列式：det(A) 矩阵的转置：B=A; 矩阵的秩：rank(A); 改变矩阵的形状：resphere(A,n,m);将数组A转变成n行m列的矩阵，其中。 2.5、矩阵的特征值、特征向量和特征多项式 eig(A)返回矩阵A的n个特征值； [v,d]=eig(A)返回以矩阵A的特征值为列的矩阵V和以矩阵A的特征值为对角元素的矩阵D。 poly(A)生成矩阵A的特征多项式，矩阵A的特征多项式的根，实际就是矩阵A的特征值。 2.6、特殊的矩阵的产生 ones全部数都为1的矩阵； zeros全部数为0的矩阵； eye单位矩阵； 以ones为例： ones(n)，n阶全为1的方阵，ones(n,m),n行，m列的矩阵。 diag(A):对角方阵，A为一维数组； magic(n)：魔术数组 大家可以验证一下，魔术数组每行，每列，以及对角线的的算数和是相等的。 2.7、矩阵的旋转 rot90 矩阵旋转90度 fliplr 矩阵的左右翻转 flipud 矩阵的上下翻转 2.8、矩阵中元素的引用 A(i,j)表示矩阵A的第i行第j列元素。A(i,：)矩阵的第i行的数值；A(：,j)矩阵A的第j列的数值。 2.9、矩阵的基本运算 B=A+1：表示矩阵A中的每一个元素都加上1。 B=n*A：表示矩阵A中的每一个元素都乘以n。 B=A*C：矩阵A和矩阵C相乘。 2.10、在求解多元一次方程组上的应用 解方程组： A=[1 2;4 1]; %系数矩阵 B=[5;6]; %结果矩阵 C=A\B %求解C = [1 2] ，即 2.11、区分普通算术运算和点(.)运算 点运算就是数组或矩阵中的每一个元素都参加运算。 2.12、MATLAB的注释：% 注释的存在意义是为了让别人或者自己以后容易懂