📄 5.htm
字号:
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=gb2312">
<title>CTerm非常精华下载</title>
</head>
<body bgcolor="#FFFFFF">
<table border="0" width="100%" cellspacing="0" cellpadding="0" height="577">
<tr><td width="32%" rowspan="3" height="123"><img src="DDl_back.jpg" width="300" height="129" alt="DDl_back.jpg"></td><td width="30%" background="DDl_back2.jpg" height="35"><p align="center"><a href="http://202.112.20.132"><font face="黑体"><big><big>白云黄鹤★</big></big></font></a></td></tr>
<tr>
<td width="68%" background="DDl_back2.jpg" height="44"><big><big><font face="黑体"><p align="center"> 遗传算法 </font></big></big></td></tr>
<tr>
<td width="68%" height="44" bgcolor="#000000"><font face="黑体"><big><big><p align="center"></big></big><a href="http://cterm.163.net"><img src="banner.gif" width="400" height="60" alt="banner.gif"border="0"></a></font></td>
</tr>
<tr><td width="100%" colspan="2" height="100" align="center" valign="top"><br><p align="center">[<a href="index.htm">回到开始</a>][<a href="index.htm">上一层</a>][<a href="6.htm">下一篇</a>]
<hr><p align="left"><small>发信人: coy (我爱美好哦), 信区: Algorithm <br>
标 题: 遗传算法（5） <br>
发信站: 武汉白云黄鹤站 (Wed May 31 21:42:11 2000), 站内信件 <br>
<br>
//det_inv.c <br>
<br>
#include "f:\coy\mywork\lunwen\gaopc\thesis.h" <br>
void inverse(a,a_inverse,n) <br>
MATRIX a,a_inverse; <br>
int n; <br>
{ <br>
int i,j,k,k1,l,l1,n1; <br>
float d1, <br>
v,w; <br>
MATRIX b; <br>
b = matrix(1, n, 1, n); <br>
for(i=1; i<=n; i++) <br>
for(j=1; j<=n; j++) <br>
b[i][j] = 0.0; <br>
d1 = det(a, n); <br>
if(d1==0.0) printf("No inverse exists"); <br>
else <br>
{ <br>
{ <br>
v=-1.0; <br>
for(i=1;i<=n;i++) <br>
{ <br>
v=-v; <br>
w=-1.0; <br>
for(j=1;j<=n;j++) <br>
{ <br>
w=-w; <br>
n1=n-1; <br>
for(k=1;k<=n1;k++) <br>
{ <br>
k1=k; <br>
if(k>=i) k1=k+1; <br>
for(l=1;l<=n1;l++) <br>
{ <br>
l1=l; <br>
if(l>=j) l1=l+1; <br>
b[k][l]=a[k1][l1]; <br>
} <br>
} <br>
a_inverse[j][i]=det(b,n1)/d1*v*w; <br>
} <br>
} <br>
} <br>
} <br>
free_matrix(b,1,n,1,n); <br>
} <br>
float det(input_matrix, nl) <br>
MATRIX input_matrix; <br>
int nl; <br>
{ <br>
int ia,ib,i1,j1,n2,nz,n3,l,l1,m,m1; <br>
float d,x; <br>
MATRIX temp_input_matrix,b1; <br>
temp_input_matrix = matrix(1,nl,1,nl); <br>
b1 = matrix(1,nl,1,nl); <br>
for(i1=1;i1<=nl;i1++) <br>
for(j1=1;j1<=nl;j1++) <br>
temp_input_matrix[i1][j1]=input_matrix[i1][j1]; <br>
d=1.0; <br>
for(ia=1;ia<=nl;ia++) <br>
{ <br>
x=-1.0; <br>
n2=nl-ia+1; <br>
n3=n2-1; <br>
n3=n2-1; <br>
nz=0; <br>
for(ib=1;ib<=n2;ib++) <br>
if(temp_input_matrix[ib][1]!=0.0) nz=ib; <br>
if(nz==0.0) d=0.0; <br>
else <br>
{ <br>
for(i1=1;i1<=nz;i1++) <br>
x=-x; <br>
d*=temp_input_matrix[nz][1]*x; <br>
for(l=1;l<=n3;l++) <br>
{ <br>
l1=l; <br>
if(l>=nz) l1=l+1; <br>
for(m=1;m<=n3;m++) <br>
{ <br>
m1=m+1; <br>
b1[l][m] = temp_input_matrix[l1][m1] - temp_input_matrix[l1][1] <br>
/temp_input_matrix[nz][1]*temp_input_matrix[nz][m1]; <br>
} <br>
} <br>
} <br>
for(i1=1;i1<=n3;i1++) <br>
for(j1=1;j1<=n3;j1++) <br>
temp_input_matrix[i1][j1]=b1[i1][j1]; <br>
} <br>
free_matrix(b1,1,nl,1,nl); <br>
free_matrix(temp_input_matrix,1,nl,1,nl); <br>
<br>
return(d); <br>
} <br>
<br>
-- <br>
＼○／ <br>
┃ <br>
／\ <br>
脖子扭扭，屁股扭扭，大家一起来跳舞 <br>
<br>
※ 来源:．武汉白云黄鹤站 bbs.whnet.edu.cn．[FROM: 211.69.196.11] <br>
[阅读精华区资料] 结束Q, ← │ 上一项资料 U,↑│ 下一项资料 <Enter>,<Space>,↓ <br>
<br>
</small><hr>
<p align="center">[<a href="index.htm">回到开始</a>][<a href="index.htm">上一层</a>][<a href="6.htm">下一篇</a>]
<p align="center"><a href="http://cterm.163.net">欢迎访问Cterm主页</a></p>
</table>
</body>
</html>⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -