matlab 变分不等式,求解变分不等式的matlab程序

满意答案

kide754861

2013.07.11

采纳率：46% 等级：11

已帮助：5599人

function x=Porjection()

clc

sigama=0.5;

gama=0.8;

x0=[10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 0 0 0 0 0 0 0 0];

x=x0;

fun(x)

z=max_fun(x);

r=(x-z);

t=1;

while norm(r)>1e-3

norm(r)

m=0;

f=fun(x);

xx=x-gama^m*r;

ff=fun(xx);

while sum((f-ff).*r)>sigama*norm(r)^2

m=m+1;

xx=x-gama^m*r;

ff=fun(xx);

if m>200

break;

end

yita=gama^m;

y=(1-yita)*x+yita*z;

d=-(yita*r-yita*fun(x)+fun(y));

alfa=yita*sum(r.*(r-fun(x)+fun(y)))/(norm(d)^2);

x=max(0,x+alfa*d);

z=max_fun(x);

r=(x-z);

t=t+1;

end

function z=max_fun(x)

f=fun(x);

xx=x-f;

z=max(xx,0);

10分享举报

matlab 变分不等式,求解变分不等式的matlab程序相关推荐

Matlab遗传算法求解TSO,遗传算法matlab程序实例.doc
遗传算法matlab程序实例.doc --------------------------------------------------------------------------------- ...
matlab ode45求解齿轮动力学,Matlab拟合动力学参数遇到问题（ode45）
function [rho2,F] = rho2_F(k,r,ss,Ne,Np) sy = sum(r.^2); rho2 = 1 - ss./sy F = (sy - ss).*(Ne-Np)./( ...
matlab如何求解代数方程,用Matlab求解复杂代数方程
编写一个function,保存为KMN.m function [K1,M1,N1,K2,M2,N2,K3,M3,N3]=KMN(x,y,z,R,r,L,l) K1=(power(x,2)+power( ...
用matlab解决向前、向后和改进欧拉法求解微分方程（编写程序）
分别使用向前欧拉法.向后欧拉法与改进dy欧拉法求解微分方程(编写程序) 在matlab中:Π用pi表示编写程序 function [y1,y2,y3,y4]=exam4_6(a,b,k) h=(b- ...
matlab中存档算法代码,MATLAB 智能算法超级学习手册中程序代码
[实例简介] MATLAB 智能算法超级学习手册中程序代码 [实例截图] [核心代码] dc90ef43-7920-434e-bdb8-0636c31c0b44 └── MATLAB 智能算法超级学习 ...
一节双曲型方程基于MATLAB的求解,双曲方程基于matlab的数值解法
<双曲方程基于matlab的数值解法>由会员分享,可在线阅读,更多相关<双曲方程基于matlab的数值解法(9页珍藏版)>请在人人文库网上搜索. 1.双曲型方程基于MATLAB ...
fama matlab源码_用matlab程序做Fama-MacBeth回归的代码
求[100,200]之间第二个能被15整除的整数.(用Matlab程序做) x=100:200;>>n=find(rem(x,15)==0);>>X=x(n(2))很简单的两句 ...
列主元Gauss消元法求解线性方程组的MATLAB实现
目录顺序高斯消去法求解线性方程组的MATLAB实现列主元Gauss消元法求解线性方程组的MATLAB实现一.简介 MATLAB实现列主元Gauss消去法求解线性方程组,并与不列主元的Gauss消 ...
【优化调度】基于matlab遗传算法求解公交车调度排班优化问题【含Matlab源码 2212期】
⛄ 一. 遗传算法简介 1 引言公交排班问题是城市公交调度的核心内容,是公交调度人员.司乘人员进行工作以及公交车辆正常运行的基本依据.行车时刻表是按照线路的当前客流量情况,确定发车频率,提供线路车辆 ...

matlab 变分不等式,求解变分不等式的matlab程序

matlab 变分不等式,求解变分不等式的matlab程序相关推荐

最新文章

热门文章