Matlab函数与方程求解1.线性方程求解2.非线性方程求解3.多项式求解4.函数创建和应用5.一元函数的零根6.符号微积分7.Taylor级数展开8.符号公式的替换与数值演算9.常微分方程10.高阶微分方程

77 阅读 0 评论 51 点赞

我是靠谱客的博主心灵美秋天，最近开发中收集的这篇文章主要介绍Matlab函数与方程求解1.线性方程求解2.非线性方程求解3.多项式求解4.函数创建和应用5.一元函数的零根6.符号微积分7.Taylor级数展开8.符号公式的替换与数值演算9.常微分方程10.高阶微分方程，觉得挺不错的，现在分享给大家，希望可以做个参考。

概述

关于Matlab的一些基本操作在另一篇博客里有提及

https://blog.csdn.net/Netceor/article/details/95043332

想查找相关内容也可在页面按ctrl+F搜索文字

1.线性方程求解

2.非线性方程求解

3.多项式求解

4.函数创建和应用

5.一元函数的零根

6.符号微积分

7.Taylor级数展开

8.符号公式的替换与数值演算

9.常微分方程

10.高阶微分方程

1.线性方程求解

指令x=Ab可求线性方程组Ax=b的特解。
(假如A不是可逆阵,即方程组有可能是欠定或超定的,则x是众多解中(含0最多)的一个或最小二乘解。)

用x=lsqnonneg(A,b)可求出方程组Ax=b的一个非负最小二乘解。

$2x+y+3z=5$
$x+3y+2z=10$

A=[2 1 3;1 3 2]; b=[5;10];
x=lsqnonneg(A,b)

求通解可对增广矩阵用化最简阶梯形rref([A,b])实现高斯消元法。

2.非线性方程求解

MATLAB求‘Fun(x)=0’解的指令一般格式：

[x,fval,exitflag,output]=fsolve(Fun,x0,options,p1,p2,...)

示例：求解方程组

$sin(x)+y+z^{2}e^{x}-4=0$
$x+xy=0$
$xyz=0$

%创建函数funs
function F=funs(x)
F(1)=sin(x(1))+x(2)+x(3)^2*exp(x(1))-4;
F(2)=x(1)+x(1)*x(2);
F(3)=x(1)*x(2)*x(3);

%解方程
funs=@(x) [sin(x(1))+x(2)+x(3)^2*exp(x(1))-4; x(1)+x(1)*x(2); x(1)*x(2)*x(3)]
[x,fv]=fsolve(funs,[1,1,1])

%输出
x =
0.0000
0.1291
1.9675
fv =
1.0e-14 *
-0.1332
0.0110
0.0025

3.多项式求解

p=[1 -7 10]
%表示x²-7x+10
r=roots(p)
%求根
p=poly(r)
%由根创建多项式
A=ones(3) , poly(A)
%方阵A的特征多项式
d=eig(A),[v,d]=eig(A)
%求方阵A的特征值d与特征向量
polyval(p,[5 0;1 2])
%多项式求值
polyval(p,x)
将x代入p
d=polyder(p)
%多项式的的导数
a=[1 -2],b=[1 -5], c=polyder(a,b)
%多项式a和b乘积的导数
[q,d]=polyder(a,b)
%多项式的商a/b的导数
%用最小二乘法对定点(xi,yi)按n阶多项式拟合
p=polyfit(x,y,n)
x=[1 2 3 4 5],y=[5.5 43.1 128 290.7 498.4]
p=polyfit(x,y,3)

4.函数创建和应用

匿名函数y=@(x,a) ... 是Matlab为定义较简单的用户函数(对应一个形式返回参数的表达式)所提供的方便法门;调用时用y(x,a)取得实际返回值。另外，Matlab中对多元函数习惯用向量函数形式表示。

%MATLAB中求y=f(x)在指定区间[a,b]上的局部极小值
fun=@(x)
2*sin(x)-1
%创建与x变量有关的名为fun的函数
[x,f]=fminbnd(fun,3,6)
%返回取极小值时自变量值x与函数值f
%求多元函数在向量x0所表示的点附近的极小值
%[x,f]=fminsearch(fun,x0) 或[x,f]=fminunc(fun,x0)
funf=@(x) x(1)^2+2.5*sin(x(2))-x(1)*x(2)^2*x(3)^2;
[x,f]=fminsearch(funf,[1,-1,0])
%Matlab中对多元函数习惯用向量函数形式表示
z=@(x,a,b) a*sin(x(1))+b*cos(x(2))
[x,f]=fminsearch(@(x) z(x,2,1),[0,0])
%将一个新匿名函数做参数传入

5.一元函数的零根

查找fun函数在x0附近的零点x及对应的函数值fval，可返回

当x0=[a,b]则表示求[a,b]内的根

[x,fval,exitflag,output] = fzero(fun,x0,options)

x = fzero(@cos,[1 2]) %求cos(x)在1和2之间的零点
x=fzero(@(x)x.^3-2*x-5,2) %求f(x) = x^3-2*x- 5在2附近的根
y=@(x) fzero(@(y) (exp(y)+x^y)^(1/y)-x^2*y,1)
f=arrayfun(y,1:20)

6.符号微积分

syms x y c
% 定义了三个基本的符号对象(变量)
f=2*x^2+y^3-c*x*y
% 由基本符合对象组成符号表达式
limit(f,y,1)
%求极限 limit(f,v,a)
diff(f,x), diff(f,y,2)
%求导数、偏导数 diff(f,v,n)，默认一阶
int(f,y)
%对f关于变量y的原函数
s=int(f,x,1,2)
%求定积分
S = solve(eqn,var)
%非线性方程组符号求解
%示例
syms a u v x
S= solve(a*u^2 + v^2 == 0, u - v == 1, a^2 + 6 == 5*a)
%S为结构变量
S.u,S.v,S.a
S=solve(x^3-a,x)
S=solve(x^3-1,'Real',true)

7.Taylor级数展开

taylor(f) 求符号函数f对默认变量的5阶马克劳林展开

taylor(f,v,‘ExpansionPoint’,a,‘Order’,n) %展开到(v-a) ^(n-1)项

syms x y t u v a b
M5=taylor(exp(-x))
f=a*sin(x)*y^x+u*cos(v);
g=exp(v)+b*v^u;
T7=taylor(f,'Order',8,'ExpansionPoint',1)
T2= taylor(g,v,'ExpansionPoint',a,'Order',3)

8.符号公式的替换与数值演算

为了使表达式简洁易读,常用用符号变量替换子表达式， MATLAB提供指令subexpr、 subs实现对表达式的替换

syms x y;
f=x^2*y*sin(x-y);
ES=diff(diff(f,y),x,2)
[RS,ssub]=subexpr(ES,'ssub') %%用符号变量ssub置换相同的表达式,重写ES为RS
RES=subs(ES,x-y,'ttt')	% 用ttt置换ES中x-y,产生RES。

对符号推演得到的公式，常需代入具体参数计算其数值，这是可用串演算指令eval('字符串表达式或符号表达式')，其功能是将字符串或符号公式当作Matlab指令、表达式执行。

x=1,y=1
dzdxy11=eval(ES)

9.常微分方程

常微分方程：含有未知函数及其导数或微分的方程。若未知函数是一元函数则称常微分方程。微分方程中出现的最高阶导数为n阶，则称n阶微分方程。

微分方程的解：代入微分方程后，使微分方程成为恒等式的已知函数。n阶微分方程的解函数中含n个独立任意常数，则称通解；若通解中的任意常数都被确定，则称特解。

Matlab求常微分方程解析解指令：

dsolve('deq1==0,...,deqn==0','...定解条件...'[,'指定自变量'])

y=dsolve('D2y+2*Dy-3*y-cos(x)==0','y(0)=4/5,Dy(0)=11/10','x')
S=dsolve('Dx==2*y,D2y=Dy/t'), S.x,S.y
%函数默认自变量为t

求公式解较困难，大多无公式解；实用中常求近似的数值解，即用数表法表示解函数

常微分方程的数值解

通常用法：

[T,Y] = solver(@odefun,tspan,y0)

sol = solver(odefun,[t0 tf],y0...)

其中：求解函数是 ode45, ode23, ode113, ode15s, ode23s, ode23t, or ode23tb中的一个。

odefun：微分方程y’=f(t,y)的右端函数.

tspan：指定数值解中时间t的范围[t0，tf]，若要获取指定时刻的数值解y的值，用 tspan = [t0,t1,...,tf]来描述指定时刻。（可以是区间，也可以是点）

y0 ：描述初始条件的向量。（tspan左端点对应的函数值，即y1(0)=0,y2(0)=1,y3(0)=1）

[T,Y]：T为返回数值解中的自变量取值（向量），Y为返回数值解中相应时间T的Y值（向量）

dy=@(t,y) [y(2) * y(3);-y(1) * y(3);-0.51 * y(1) * y(2)]
[T,Y] = ode45(dy,[0 12],[0 1 1]);
%t取值范围：0-12，端点函数值分别0、1、1
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
%figure1
%[T,Y] = ode45(@odef,[0 12],[0 1 1]);
%plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')