义章(绪论).docVIP

绪论 1.1 概述 MATLAB---Matrix Laboratory是矩阵实验室的意思。 1.2 MATLAB的使用方法 clear , clear x,y,z ,clear clc (清除窗口命令。) 文件扩展名为：.m 1.3 创建变量名 用户创建的变量名不可超过31个字符。变量名必须以字母开头，后面可跟大小写字母、数字或下划线，字符间不可有空格、减号等。 1.4 命令窗口管理 format compact---剔出显示中多余的空行或空格。 format long e----将系统默认的5位数字的显示格式转换位16位数字加3位指数位的显示格式。 format short-----恢复统默认设置。 ^c ----用于终止程序或函数的执行。 ^p ----用于将最近一次键入的信息显示在MATLAB命令窗口中，按Enter键可再次执行该命令。按Enter键执行前，可修改调入的命令行。 1.5 MATLAB的基本语法 MATLAB要求表达式在使用前一定要先对其进行赋值操作。注意三点： 1、无逗号（，）或分号（；） 2、有逗号（，）无分号（；） 3、只有分号（；） 例： 在MATLAB中可写为： t=(1/(1+p*x))^k 在执行本句之前首先要对p、x、k赋值。 MATLAB的语法举例： 1-d*c^(x+2) d*c^x 或 2+d*c^x (2/d)*c^(x+2) 或 2/d*c^(x+2) 或 2*c^(x+2)/d (d*c^x+2)/g^2.7 sqrt(d*c^x+2) 或 (d*c^x+2)^0.5 y=sqrt(abs(pi-sin(x)/cosh(a)-ln(x+a))) 练习 1.1 一个长管，其内表面半径为a，温度为Ta；外表面半径为Tb；则其径向和切向应力可分别表示为： 式中r为管子的径向坐标，E为管子材料的弹性模量，ɑ为热膨胀系数。沿径向管壁的温度分布为： 用MATLAB表达形式写出该方程并用如下值进行验证：ɑ=1.2x10-5，E=3x107，v=0.3,Ta=500，Tb=300 ,a=0.25，b=0.5，r=0.375。 [答案：= -8011.5，= 5231.9, =383.0075] 1.2 压力为po的气体在可逆绝热条件下从容器中溢出,其质量与流量之比如下: 式中pe为容器出口的外部压力,k为可逆绝热的气体常数。用MATLAB表达形式写出该方程并用如下给出值进行验证：k=1.4，pe/p0=0.3。 [答案：=0.4271] 1.3 依据如下给出的公式，当n=0时可在10-7的误差范围内算出的近似值，而当n取两项（n=0,1）时其误差范围则为10-15。实际上，n的取值每增加一个，近似值的精度将提高10-8。因此，n取4项（n=0,1,2,3）时值可计算到31位，可用MATLAB符号工具箱进行验证。 注意：阶乘的???算可采用下述gamma函数：gamma(n+1)=n!；gamma(4n+1)=(4n)!,n=0,1,…。