matlab subs eval,【荐】MATLAB求解符号表达式数值的方法：subs函数

大家看一下这个：function [xmin,funval]=TR(fun,x0,delta0,x)

delta=delta0;

eps=1.0e-4;

n=length(x0);

x0=x0';%x0为列向量

H=eye(n);

mu=0.25;

yita=0.75;

normg=1;

while normg>eps

fx=subs(fun,x,x0);

gradfun=jacobian(fun,x);

g=subs(gradfun,x,x0);%g为行向量

normg=norm(g);

s0=x0';%s0为行向量

[s,fs]=trustP(fun,g,H,delta,x,s0);%s应该为行向量

s=s';%s变为列向量

x1=x0+s;%x1为列向量

fx1=subs(fun,x,x1);

rho=(fx-fx1)/(-fs);

if rho>yita

x0=x0+s;%x0为列向量

else

x0=x0;

end

if rho<=mu

delta=0.5*delta;

else

if rho>=yita

delta=2*delta;

else

delta=delta;

end

g1=subs(gradfun,x,x1);%g为行向量1

y=g1-g;%y为行向量

y=y'; %ye为列向量

dyT=y';%dyT为行向量

dsT=s';%dsT为行向量

A=s*dsT;

B=H*y*dyT*H;

a=dsT*y;

b=dyT*H*y;

H=H+A/a+B/b;

end

xmin=x0;

funval=subs(fun,x,xmin)

end

function [s,fs]=trustP(fun,g,H,delta,x,s0)

epsi=1.0e-4;

r0=8;

beta=0.5;

dxT=x';

B1=x*H*dxT;

c=g*dxT;

x0=s0;

fx=subs(fun,x,x0);

sfun=fx+c+B1;%信赖域子问题函数

N=length(x);

qfun=0;

for k=1:N

qfun=qfun+x(k)^2;

end

sqfun=sqrt(qfun);

cfun=delta-sqfun;%罚函数

tol=1;

while tol>epsi

pfun=r0*cfun;

gfun=sfun+pfun;%障碍函数

[s1,sf]=NT(gfun,x,s0)

cf=subs(cfun,x,s1);

tol=r0*cf;

r0=r0*beta;

end

s=s1;

fs=subs(fun,x,s);

end

function [s1,sf]=NT(gfun,x,s0)

eps=1.0e-4;

tol=1;

s0=s0';%s0为列向量

while tol>eps

gradf=jacobian(gfun,x)

jacf=jacobian(gradf,x)

v=subs(gradf,x)

tol=norm(v);

pv=subs(jacf,x,s0

p=-inv(pv)*v';

ss1=s0+p;

s0=ss1;

end

s1=ss1;

sf=subs(fun,x,s1);

end

这三个函数是我自己编写的，但是运行时subs函数会报错！有谁能帮我解决一下！

matlab subs eval,【荐】MATLAB求解符号表达式数值的方法：subs函数相关推荐

单位斜变函数matlab,实验二用MATLAB实现线性系统的时域分析
实验二基于MATLAB的线性系统时域分析 [实验目的] 1．研究线性系统在典型输入信号作用下的暂态响应: 2．熟悉线性系统的暂态性能指标: 3．研究二阶系统重要参数阻尼比ξ对系统动态性能的影响: 4． ...
MATLAB解隐函数方程时符号表达式转化为数值的方法-用vpa函数
今天在解决一个小问题时,遇到解隐函数方程,中间涉及一个解的传递问题,才好好研究了一下 syms这个语句的一些语法规则,最终用vpa这个函数解决了符号表达式到double数值的转化. syms是符号函数 ...
matlab subs eval,matlab syms subs eval 实例
clc; close all; clear; syms T; syms t; syms vmax; len = 150; y = vmax*(1 - exp(-t/T)) y2 = vmax y3 = ...
求解非线性方程f (x)= 0的MATLAB数值法指令介绍（solve、fzero的方法与实例）
一.符号方程求解在MATLAB中,求解用符号表达式表示的代数方程可由函数solve实现,其调用格式为: solve(s)(求解符号表达式s的代数方程,求解变量为默认变量,当方程右端为0 ...
隐式函数matlab,matlab隐函数求解的几种方法
1. fzero() 2. fsolve()+ inline() 3. solve() 4. 先求解析解,后用替代法求数值解. A=[sinxcosy, e^siny], B=[siny,1+e^(x ...
matlab多项式加法运算,matlab多项式运算与代数方程求解解析.ppt
* 多项式运算与代数方程求解数学软件 Matlab Matlab基础及应用 * 多项式转化为符号表达式:poly2sym 四则运算:conv.deconv 导数与积分:ployder.polyint ...
fick第二定律matlab求解,(最新整理)实验二MATLAB数值计算常微分方程(组)的求解
<(最新整理)实验二MATLAB数值计算常微分方程(组)的求解>由会员分享,可在线阅读,更多相关<(最新整理)实验二MATLAB数值计算常微分方程(组)的求解(19页珍藏版)> ...
matlab解带参数方程,matlab求解变参数方程
求解含参数方程组以及带入数值一.求解含参方程组举个简单例子,解方程组 x+A*y=10 x-B*y=1? ??其中 x,y 为变量,A,B为字母系数. 只要在 Matlab中...... (三)用 ...
matlab绘制圆，且求解两个圆的交点坐标
前言这里简单记录下对于matlab绘制圆,以及求解两圆交点的坐标的实现(包括C++). 一绘制圆 circle.m文件: function [] = circle(x, y, r, color) ...
matlab微分方程组边值,matlab求解常微分方程边值问题的方法
matlab求解常微分方程边值问题的方法 Matlab 求解常微分方程边值问题的方法:bvp4c 函数常微分方程的边值问题,即 boundary value problems ,简称 BVP 问题, ...

matlab subs eval,【荐】MATLAB求解符号表达式数值的方法：subs函数

matlab subs eval,【荐】MATLAB求解符号表达式数值的方法：subs函数相关推荐

最新文章

热门文章