Matlab精讲多练5

subplot(3,1,1);
plot(t,x)
%位置
hold on;
subplot(3,1,2);
plot(t,v)
%速度
hold on;
subplot(3,1,3);
plot(t,a)
%加速度
hold on;
end
结果如下：
%起始状态x v a
xe = [te^5 te^4 te^3 te^2 te 1]; ve = [5*te^4 4*te^3 3*te^2 2*te 1 0]; ae = [20*te^3 12*te^2 6*te 2 0 0];
%结束状态x v a
A=[xs;vs;as;xe;ve;ae]; t=ts:0.01:te; for i=-5:5
可以根据状态量位置速度加速度的起始和结束值列出6个方程组成方程组解该问题
matlab练 习 程 序 （ 五 次 多 项 式 轨 迹 规 划 ）
可以根据状态量(位置，速度，加速度)的起始和结束值列出6个方程，组成方程组解该问题。 1. 列出起始状态：
2. 列出终止状态：
3Hale Waihona Puke 写成矩阵形式：求解c即可。
下面是从横向-5米到5米的生成的路径。
代码如下：
clear all;close all;clc;
ts = 0;te = 10;
%起始结束时间
xs = [ts^5 ts^4 ts^3 ts^2 ts 1]; vs = [5*ts^4 4*ts^3 3*ts^2 2*ts 1 0]; as = [20*ts^3 12*ts^2 6*ts 2 0 0];
B=[0 0 0 i 0 0]'; %[xs vs as xe ve ae]

语句组m
otherwise
语句组n
end
（exswitch.m）
第五讲 MATLAB程序设计
18
（3）try语句 语句格式为： try
语句组1 catch
语句组2 end
try语句先试探性执行语句组1，如果语句组1在执行 过程中出现错误，则将错误信息赋给保留的lasterr 变量，并转去执行语句组2。
第五讲 MATLAB程序设计
14
2、选择结构
（1） 条件分支语句——if语句 在MATLAB中，if语句有3种格式。 1) 单分支if语句： if 条件 语句组
end
第五讲 MATLAB程序设计
15
2) 双分支if语句： if 条件
语句组1 else
语句组2 end
第五讲 MATLAB程序设计
16
第五讲 MATLAB程序设计
24
三、程序调试
1 错误分类
一般来说，应用程序的错误有两类：
一类是语法错误，例如函数名的拼写错、表达式 书写错等。
另一类是运行时的错误。指程序的运行结果有错 误，这类错误也称为逻辑错误。
第五讲 MATLAB程序设计
25
2、查找逻辑错误的方法：
◆ 删去语句行末的分号，使显示其运行中间结果 ◆ 利用keyboard 命令实现，return继续程序执行 ◆ 注释掉M 函数文件的函数定义行，使函数文件转
第五讲 MATLAB程序设计
19
例： 矩阵乘法运算要求两矩阵的维数相容，否则会 出 错。先求两矩阵的乘积，若出错，则自动转去 求两矩阵的点乘。（extry.m）
第五讲 MATLAB程序设计
20
3、 循环结构
（1）硬循环语句——for语句

2020/4/19
28
六、文件操作
matlab文件输入、输出（I/O）函数允许把其它格式 的数据文件直接读进matlab，它是以c语言文件I/O 函数为基础的。
1. 打开和关闭文件（fopen, fclose) 格式：fid=fopen(‘数据文件’，‘r’) ‘r’ —— 读 ‘w’ —— 写 ‘a’ —— 附加 ‘rt’ —— 读写
2020/4/19
27
break —— 中断
break 中断for, while循环语句的执行，在嵌套循 环结构中， break从最里层循环退出。
continue —— 继续
与for, while循环语句一起使用，结束本次循环，即 跳过循环体中尚未执行的语句，接着进行下一次是否 执行循环的判断。
各个省市对应的人口，第4列为人均GDP，第5列为各个省市对应的农村人口， %第6列为人均农业收入，第7列为各个省市对应的城镇人口，第8为人均城镇收
入，第9列为人均总收入. 第10列为宏观区域号,第11列为三大带号.
%****以下几步GDPPC/TPCI/RPCI/UPCI/为0,让其对应的人口也为0; data_belt_macro_mean_gem(find(data_belt_macro_mean_ge
2020/4/19
2
一、matlab的工作模式
指令驱动模式
通常matlab以指令驱动模式工作，即在 matlab 窗口下当用户输入单行指令时， matlab 立即处理这条指令，并显示结果，这 就是matlab命令行方式。
命令行操作时，matlab窗口只允许一次执 行一行上的一个或几个语句。
2020/4/19
3
例如： x1=0:10 x1 =

西南科技大学本科生课程备课教案计算机技术在安全工程中的应用——Matlab入门及应用授课教师：徐中慧班级：专业：安全技术及工程第四章课型：新授课教具：多媒体教学设备，matlab教学软件一、目标与要求掌握矩阵与数组的相关运算，及matlab中矩阵运算的相关函数，包括三角分解、正交变换、奇异值分解、特征值分解、矩阵的秩的运算等。

二、教学重点与难点本堂课教学的重点在于引导学生在编写matlab程序时能够熟练运用矩阵运算的相关函数实现相应的功能。

三、教学方法本课程主要通过讲授法、演示法、练习法等相结合的方法来引导学生掌控本堂课的学习内容。

四、教学内容一、课后习题的解说。

（1）在计算器发明（约1974年）之前，人们需要用数学用表来计算正弦、余弦和对数值。

创建正弦值数学用表的步骤如下：①创建角度矢量、范围在0～3600之间，步长为180。

②计算正弦值，用角度和计算出来的正弦值创建表格。

③分别用两个disp语句给表格加上标题和表头。

④用fprintf显示数据，要求小数点后有两位有效数字。

解：angle=0:18:360; sine=sin(angle/180*pi);disp(' SINE TABLE ')disp(' Angle Sine ')fprintf(' %4.2f %4.2e\n',[angle;sine])（2）使用搜索引擎或浏览器搜索英镑、日元、欧元和人民币对美元的汇率，并把输出结果绘制成表。

要求用disp在表格中添加标题和表头，用fprintf输出格式化数据。

①创建日元和美元的汇率表，表中共有25行，从5日元开始，步长为5日元②创建人民币和美元的汇率表，表中共有30行，从5元开始，步长为5元③创建数据表格，表中有5列，第一列是美元，第二列是欧元，第三列是英镑，第四列是人民币，第五列是日元。

计算与1到10美元等价的其它货币值。

（将结果输出到.txt文件中，此步骤属选做）解：①jpy=5:5:25*5;usd1=jpy*0.01301;disp(' JPY &USD TABLE ')disp(' JPY USD ')fprintf(' %4.2f %4.2f\n',[jpy;usd1])②cny=5:5:30*5;usd2= cny *0.1567;disp(' CNY &USD TABLE ') disp(' RMB USD ') fprintf(' %4.2f %4.2f\n',[cny;usd2])③usd=1:1:10;eur=usd* 0.7323; gbp=usd* 0.6405; cny=usd* 6.3816; jpy=usd*76.358;disp(' AS Exch')disp(' USD EUR GBP RMB JPY')fprintf(' %4.2f %4.2f %4.2f %4.2f %4.2f \n',[ usd;eur;gbp;cny;jpy])二、矩阵的相关知识掌握矩阵与数组的相关运算，及matlab 中矩阵运算的相关函数，包括三角分解、正交变换、奇异值分解、特征值分解、矩阵的秩的运算等。