4例6.10 数值方法解 e-3sin(4t+2)+4e-0,5cos2t=0.5 使用slve()函数: ∈>> syms t x solve(exp(-3*t)*sin(4*t+2)+ 4*exp(-0.5*)*COS(2*)-0.5) ↓使用图解法求初始值: >>y=@(七)exp(-3*七).*sin(4*七+2)+ 4*exp(-0.5*t).*COS(2*t)-0.5; [,f]= fsolve(y,3.5203) 2/20/2021星期六, 2008-9-6,13:09:29 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:09:29 例 6.10 数值方法解 使用solve()函数： 使用图解法求初始值：
重新设置相关精度的控制变量 k > ff=optimset; ff. TolX=le-16 ff. TolFun=1e-30 It, f]=fsolve(y, 3. 5203, ff) ↓所期望的精度可能无法达到 ↓然而,在算精度制下的最好结果可以得到 2/20/2021星期六, 2008-9-6,13:09:29 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:09:29 重新设置相关精度的控制变量 所期望的精度可能无法达到 然而，在算精度制下的最好结果可以得到
61.4非线性矩阵方程求解 Riccati方程(第4章) AX+XA-XBX+C=0 更多的非线性矩阵方程,例如, e广义 Riccati方程 AX+XD-XBX+C=0 e类 Riccati方程 AX+XD-XBX+C=0 还有很多很多的矩阵方程 2/20/2021星期六, Slide 1(of 11) 高等应用数学问题的 MATLAB求解 2008-9-6,13:09:29 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:09:29 6.1.4 非线性矩阵方程求解 Riccati 方程（第4章） 更多的非线性矩阵方程，例如， 广义Riccati方程 类Riccati方程 还有很多很多的矩阵方程
Riccati方程求解 ↓函数 fsolve()只能求解出∫(m)=0,而不 是F(X)=0,其中,∫,c是向量不是矩阵 将矩阵方程转换成向量方程 向量转换成矩阵 数学表述→X ◆ MATLAB表述X= reshape(m,size(A)) 矩阵转换成向量 ◆数学表述 F 令 MATLAB表述f=F(:) 2/20/2021星期六, Slide 1(of 11) 高等应用数学问题的 MATLAB求解 2008-9-6,13:09:29 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:09:29 函数fsolve()只能求解出 ，而不 是 ，其中， 是向量不是矩阵 将矩阵方程转换成向量方程 向量转换成矩阵 数学表述 MATLAB表述 矩阵转换成向量 数学表述 MATLAB表述 Riccati方程求解
以向量的形式描述 Riccati方程的残差 function y=new_are(x, A, B,C) Ⅹ= reshape(x,Size(A) y1=A2*X+X*A-X*B*x+C;y=y1(:); 求解 Riccati方程的新函数 function x=solve are(A,b,c, xo) if nargin==3, xO=rand (size (a)); end x=fsolve(@new are, xo(:),[,A,B, c) X=reshape(x, size(A)) 2/20/2021星期六, 2008-9-6,13:09:29 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:09:29 以向量的形式描述Riccati方程的残差 求解Riccati方程的新函数