数值计算方法实验二.pptVIP

*第1页，共21页，星期日，2025年，2月5日上课时间10-16周,周三下午14：00–16：30(课间不休)*第2页，共21页，星期日，2025年，2月5日实验报告要求实验报告在课后一周内上交；算法使用编程语言不限；实验报告可以提交书面试验报告，也可以提交电子版试验报告；电子版试验报告与程序代码发送至服务器如下文件夹：张昆\学生作业\实验x文件名称：实验x张三*第3页，共21页，星期日，2025年，2月5日实验二题目本次实验要求至少完成两个实验题目：一、自选试验题目：数值试验三：（1）、（2）、（3）、（4）注：实验题目在指定范围内可以自拟*第4页，共21页，星期日，2025年，2月5日Matlab基础（实验一）基本命令基本数学运算符号解函数文件格式函数调用文件保存位置*第5页，共21页，星期日，2025年，2月5日基本命令1.help帮助命令helpformat%查看format命令的帮助；2.clear清除变量命令clearx%清除变量x(无参数时清理工作空间）3.清理命令窗口clc%清屏4.format输出格式formatlong%15位有效数字格式输出formatshort%4-7位有效数字格式输出formatcompact%紧凑格式(输出不加空行）*第6页，共21页，星期日，2025年，2月5日基本数学运算乘法：*除法：/乘方：^根号：sqrt(x)正弦：sin(x)余弦：cos(x)自然对数：log(x)自然指数：exp(x)*第7页，共21页，星期日，2025年，2月5日1.vpa控制运算精度vpa(表达式,运算精度)vpa(pi,20)%显示?至100位有效数字2.solve方程的符号解solve(表达式,x)符号解*第8页，共21页，星期日，2025年，2月5日符号解例解方程解在MATLAB工作窗口输入命令:y=solve(x^3-sin(x)-12*x+1=0,x)vpa(y,5）y=roots([2,0,-1,-1])*第9页，共21页，星期日，2025年，2月5日函数文件格式函数文件由function语句引导，其格式为：function[输出形参表]=函数名(输入形参表)%注释说明部分函数体：例：function[k,xk,yk,p]=jhnewtonqx(x0,ddmax)%牛顿切线法求非线性方程的根…注：其中函数名的命名规则与变量名相同。输入形参为函数的输入参数，输出形参为函数的输出参数。当输出形参多于1个时，则应该用方括号括起来。*第10页，共21页，星期日，2025年，2月5日函数调用函数调用的一般格式是：[输出实参1,实参2,...]=函数名(输入实参1,实参2,...）例：[k,xk,yk,p]=jhnewtonqx(1.5,20)*第11页，共21页，星期日，2025年，2月5日文件保存位置设置当前目录位置cdC:\Matlab7\work%设置Matlab当前目录为：C:\Matlab7\work*第12页，共21页，星期日，2025年，2月5日Matlab基础（实验二）1.矩阵输入2.特殊矩阵3.误差分析4.内置函数*第13页，共21页，星期日，2025年，2月5日矩阵输入a=[] %建立空矩阵a,可在workspace中编辑a=[1,2,3;456] %[分行“,”或空格分列“;”]a=1:5 %建立1×5矩阵(向量)a=(1,2,3,4,5)’a=0:pi/5:2*pi%(pi=?)建立1×11矩阵（向量）a：初始元0，步长?/5,终止元2?【注意“：”的用法】*第14页，共21页，星期日，2025年，2月5日常用矩阵1、零矩阵z=zeros(3,4) %产生3×4零矩阵z；zeros(5) %产生5阶零矩阵；2、单位矩阵E=eye(6) %产生6阶单位矩阵E；2、幺矩阵ones(3,2); %产生3行2列元素为1的矩阵；3、随机矩阵