24小时热门版块排行榜

返回列表

Honey155752

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 27.6
帖子: 14
在线: 1.9小时
虫号: 6850541
注册: 2017-06-26
专业: 计算机软件

[求助] 非线性方程组求解

1.解决目标：方程组如下

CODE:

I(1) * (2 * 0.875+ R(1)) - I(2) * R(1) - 36=0 I(2) *(2 * 0.875 + R(1) + R(2))- I(1) * R(1)=0 %灯具功率表达式 I(2)^2 * R (2) - Power=0 (I(2) - I(1))^2 * R(1) - Power=0

2.经过各种资料查询，现想采用下述遗传算法解决上述非线性方程组求解问题，如下所示：

CODE:

 %定义遗传算法参数 NIND=50; %个体数目(Numbe of individuals) MAXGEN=1000; %最大遗传代数(Maximum number of generations) NVAR=4; %变量的维数 PRECI=10; %变量的二进制位数(Precision of variables) GGAP=0.9; %代沟(Generation gap) % lbi1=0; %i1取值范围[0.138,0.25] % lbi2=0; %i2取值范围[0.138,0.25] % lbr1=80; %r1取值范围[80,259.2] % lbr2=80; %r2取值范围[80,259.2] % ubi1=10; % ubi2=10; % ubr1=259.2; % ubr2=259.2; lbi1=0; %i1取值范围[0.138,0.25] lbi2=0; %i2取值范围[0.138,0.25] lbr1=50; %r1取值范围[80,259.2] lbr2=50; %r2取值范围[80,259.2] ubi1=10; ubi2 =10; ubr1=300; ubr2=300; trace=zeros(MAXGEN, 2); %建立区域描述器(Build field descriptor) FieldD=[rep([PRECI],[1,NVAR]);lbi1 lbi2 lbr1 lbr2;ubi1 ubi2 ubr1 ubr2;rep([1;0;1;1],[1,NVAR])]; %FieldD Chrom=crtbp(NIND, NVAR*PRECI); %创建初始种群 %Chrom gen=0; %代计数器 X=bs2rv(Chrom,FieldD);%种群二进制到十进制转换 %X % % X*X % X.^2 ObjV=ff(X); %目标函数 % ObjV while gen<MAXGEN %迭代 FitnV=ranking(ObjV);%分配适应度值(Assign fitness values) %FitnV SelCh=select('sus', Chrom, FitnV, GGAP); %选择 %SelCh SelCh=recombin('xovsp', SelCh, 0.7); %重组 %SelCh SelCh=mut(SelCh); %变异 %SelCh X=bs2rv(SelCh,FieldD); %种群二进制到十进制转换 %X ObjVSel=ff(X); %目标函数 [Chrom ObjV]=reins(Chrom, SelCh, 1, 1, ObjV, ObjVSel); %重插入 %X=bs2rv(Chrom,FieldD); gen=gen+1; %代计数器增加 trace(gen, 1)=min(ObjV); %遗传算法性能跟踪 trace(gen, 2)=sum(ObjV)/length(ObjV); end plot(trace(:,1));hold on; plot(trace(:,2),'-.');grid; legend(' 种群均值的变化','解的变化') %输出最优解及其对应的20个自变量的十进制值,Y为最优解,I为种群的序号 [Y, I]=min(ObjV); Y X=bs2rv(Chrom, FieldD); X(I,:)

调用的ff()函数代码如下：

CODE:

function [ y ] = ff(x) f1=2*x(:,1)*0.875+(x(:,1)-x(:,2)).*x(:,3)-36; f2=2*x(:,2)*0.875+(x(:,2)-x(:,1)).*x(:,3)+x(:,2).*x(:,4); f3=5-(x(:,1)-x(:,2)).*(x(:,1)-x(:,2)).*x(:,3); f4=5-x(:,2).*x(:,2).*x(:,4); % f1(1,1)^2+f2(1,1)^2+f3(1,1)^2+f4(1,1)^2 y=f1.^2+f2.^2+f3.^2+f4.^2; end

3.如若要求解1中的非线性方程组，可以将问题转化为求解f1^2 +f^2+f3^2+f4^2的最小值问题。并采用2中的遗传算法求解最小值，出现的问题如下：
  1）适应度函数不收敛；
  2）适应度函数不收敛，变量求解结果一直不稳定；

4.方案反思：
  1）非线性方程组转化为多元函数的方法是否不对？即是否可转化为求解f1^2+f^2+f3^2+f4^2的最小值问题？
  2）如若1）的转换方法是正确的，那么ff()函数的代码部分是否正确？
  3）如若上述的求解方案是不正确的，是否还有其他可行办法求解非线性方程组（方程组个数很多，未知量个数很多，这个才是我的终极目标，代码只是以4个
   方程展开尝试的）；

非常感谢。已被该问题折磨很多天，求帮忙解决。
目前就9枚硬币，全部奉上。。着急！@@

回复此楼

» 猜你喜欢

中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学已经有0人回复
《电磁学》教材推荐已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有138人回复
【急招】合肥工大核聚变材料计算方向2026级工程博士生已经有4人回复
大豆异黄酮分离已经有0人回复
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标已经有10人回复
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名已经有1人回复
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生已经有3人回复
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究已经有12人回复

» 本主题相关价值贴推荐，对您同样有帮助:

急求急求帮助！！！在线等！！！如何用matlab实现以下图片中非线性方程组求解？？已经有6人回复
关于方程数多于未知量的非线性方程组求解已经有0人回复
非线性方程组求解已经有0人回复
带限制的多变量非线性方程组求解已经有5人回复
非线性方程组求解已经有1人回复
牛顿迭代法接非线性方程组求解已经有1人回复
不知道matlab程序为什么错了（错误应该存在非线性方程组求解中）已经有0人回复
matlab 计算非线性方程组求解已经有4人回复
matlab非线性方程组求解，急！！！已经有5人回复
非线性方程组求解已经有0人回复
非线性方程组求解问题已经有0人回复
有偿求助非线性方程组求解收敛性证明已经有0人回复
MATLAB非线性方程组求解~~~反应釜模拟已经有9人回复
【求助】非线性方程组求解已经有4人回复
【求助】非线性方程组求解已经有3人回复
【求助】一个非线性方程组求解问题已经有6人回复

1楼 2017-07-10 16:06:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Honey155752

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 27.6
帖子: 14
在线: 1.9小时
虫号: 6850541
注册: 2017-06-26
专业: 计算机软件

引用回帖:

1楼: Originally posted by Honey155752 at 2017-07-10 16:06:39
1.解决目标：方程组如下
I(1) * (2 * 0.875+ R(1)) - I(2) * R(1) - 36=0 I(2) *(2 * 0.875 + R(1) + R(2))- I(1) * R(1)=0 %灯具功率表达式 I(2)^2 * R (2) - Power=0 (I(2) - I(1))^2 * R(1) - Power=0

2.经过 ...

这个是matlab跑出来的适应度函数图，不收敛。