matlab squrt,MATLAB数值计算和符号运算

实验五 MATLAB 数值计算和符号运算

一、实验目的

1．掌握并理解Maltab 在数值计算的基本用法。

2．理解matlab 在工程领域解题的一般过程。

3. 掌握Matlab 符合运算的基本方法。

二、实验的设备及条件

计算机一台(带有MATLAB7.0以上的软件环境)。

三、实验内容

1．线性系统方程：分别使用左除(\)和LU 分解求解下面系统方程的解：

??

???=+=+=++377251463c b b a c b a

2. 使用quad 和trapz 求解?-5

03/dx xe x 的数值积分，并与其解析解相比较(解析解利

用符号运算进行求解，参考函数int )；(要求：使用quad 求积分时，请分别用函数文件和匿名函数的方式求解)

提示：这里trapz 是离散积分，需先分别对x 和y 数据进行采样，即需要产生相应的数据点模拟采样数据)；

3.求下面函数在[0.5,4]区间内的过零点。(用fzero 函数)

提示：首先用函数ezplot 绘制曲线f(x)在[0.5,4]上的图像，命令如下 ezplot('x^3-2*x^2*sin(x)+5*x*cos(x)+1/x',[0.5,4])；其次，根据图像得到函数上接近零点的坐标，利用fzero 函数求其零点。

4.求常微分方程的数值解

提示：需先撰写函数，再用ode45或ode23求解。

5.对边长为3m 的正方形铁板，在4个角剪去相等的正方形以制成方形无盖水槽，问何种剪法可使水槽的容积最大？(提示：设剪去的正方形边长为x ，则水槽的体积应为(3-x)2

x ，利用fminbnd 求解，考虑实际情况，注意fminbnd 的求解区间的设置) ???

????-='==++-1)0(1)0(02)1(22y y y dx dy x dx y d