2020/3/3
11
2、利用非线性方程组的符号法求解如下方程组：
x1 0.7sin x1 0.2cos x2 0
x2 0.7 cos x1 0.2sin x2 0
在（0.5，0.5）附近的解。
（a）
a) x1=0.5265 x2=0.5079
b) x1=0.5079 x2=0.5265
ans =
x0=[0.5 0.5] >> fsolve('fc',x0)
0.5265 0.5079
2020/3/3
10
九、方程符号求解
1 求微分方程 y(0)=0的解
x
d2 y dx2

3
dy dx

x2
，y(1)=0，
（c）
a)x^4/5 - x^3/4 b)x^3/3 - x^2/3 c) x^4/3 - x^3/3 d) x^4/5- x^3/4
>> [x,y]=dsolve('Dx=y,Dy=x','t')
2020/3/3
9
非线性方程的符号注解
语法： fsolve(‘fun’, x0) （1）fun为所要求解的函数名，常以M文件形式给出； （2）X0为初始向量或矩阵
Fc.m function y=fc(x) y(1)=x(1)-0.7*sin(x(1))0.2*cos(x(2)); y(2)=x(2)0.7*cos(x(1))+0.2*sin(x(2)); y=[y(1) y(2)]
2020/3/3
8
【例
】 求微分方程
x
d2 y dx2

3
dy dx

x2
，y(1)=0，y(0)=0的解。
>> y=dsolve('x*D2y3*Dy=x^2','x')
%求微分方程的通解
y=
C3*x^4 - x^3/3 + C2
>> y=dsolve('x*D2y3*Dy=x^2','y(1)=0,y(0)=0','x') %求微分方程的特解

-x+10y-2x=7

-2x+10z=6 的解
2020/3/3
5
2020/3/3
6
九、方程符号求解
2020/3/3
7
符号常微分方程
MATLAB提供了dsolve命令，可以用于对符号常微分方程进行求解。
语法：
dsolve('eq', 'con', 'v') dsolve('eq1,eq2…', 'con1,con2…', 'v1,v2…')
（1）当y是因变量时，微分方程'eq'的表述规定为：
y的一阶导数
dy dx
或dy
dt
表示为Dy；
%求解微分方程 %求解微分方程组
y的n阶导数 dn y 或dn y 表示为Dny。
dxn dtn
（2）微分初始条件'con'应写成'y(a)=b，Dy(c)=d'的格式；当初始条件少于
微分方程数时，在所得解中将出现任意常数符C1、C2……，解中任意常数符的数目 等于所缺少的初始条件数。
c) x1=0.6265 x2=0.7079
d) x1=0.7079 x2=0.6265
2020/3/3
12
1 对于圆周率值，若取值为3.1415，则其
有数字为
（b）
a、3位
b、4位
c、5位 d、2位
2020/3/3
13
2020/3/3
14
2下面的MATLAB语句中正确的有：（cd）
a) a=sin(x)
2020/3/3
20
3 试用一个二次多项式拟合下列数据 x 0.5 1.0 1.5 2.0 y 1.75 2.45 3.81 4.80
2.5 8.00
3.0 8.60
（d）
a) 0.4600 1.8601 0.9730 b) 0.8400 2.9401 0.85608 c) 0.4900 1.8601 0.9730 d) 0.4900 1.2501 0.8560
(a) return; (b) break (c) continue (d) keyboad
2020/3/3
17
二、多项式及其运算
2020/3/3
18
二、多项式及其运算
1 求多项式 x4 12x3 25x 116 在x=3,8处的值
a)-52 -1730
（c）
b) -54 -1732
c) [4 19 3;21 75 54; 22 72 72]
d) [6 25 63;0 36 18; 1 15 0]
2020/3/3
23
（b） 2、求数组b=[1 2 3;4 5 6;7 8 9]的指数运算
b)Record-1=3+4i
c) a=sin(pi)
d) C=1+6j
2020/3/3
15
3 符合MATLAB变量命名规则的有a) var!205 b)VAR15（b）
c)my-var2015 d)2015_VAR
2020/3/3
16
4在循环结构中跳出循环，但继续下次循环
的命令为
（c）
c) -52 -1732
d) -54 -1730
2020/3/3
19
2 计算多项式x3+11x2+55x+125 在x=[1
2;3 4]的矩阵计算值
（a）
a)[192 287;416 585]
b)[190 285;406 582]
c)[194 292;422 580]
d)[192 288;420 580]
习题课
Any 函数
判断元素是否为非零元素any(v),如果v是 非零元素返回true(即1)否则返回flase(即 0)
如果A是矩阵，则把矩阵的列当做向量来处 理，函数返回每个列向量的逻辑值
2020/3/3
2
2020/3/3
3
Magic函数
2020/3/3
4
线性方程组求解
试求数值方程组 10x-y=9
2020/3/3
21
三、矩阵及数组代数运算
2020/3/3
22
三、矩阵及数组代数运算
1、已知A=[0 9 6;1 3 0;3 5 7];B=[1 4 3;1 5 0;2 5
9], 则A.*B
（b）
a) [21 75 54; 4 19 3; 22 72 72]
b) [0 36 18; 1 15 0; 6 25 63]
y=
x^4/3 - x^3/3
【例 】 求微分方程组
dx y, dy x dt dt
的解。
>> [x,y]=dsolve('Dx=y,Dy=x') x= C8*cos(t) + C7*sin(t) y= C7*cos(t) - C8*sin(t)
程序分析：默认的自由变量是t，C1、C2为任意常数，程序也可指定自由变量， 结果相同：