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枚硬币，全部奉上。。着急！@@

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有260人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复

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

急求急求帮助！！！在线等！！！如何用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跑出来的适应度函数图，不收敛。