|
| |
 精品推荐 |
 |
|
| |
|
|
|
|
数值计算程序大放送-数学变换与滤波
|
日期：2008年6月18日 作者： 查看:[大字体
中字体 小字体]
|
e[jj]=r[jj];
for (kk=0; kk<=n-1; kk++)
{
e[jj]=e[jj]+h[i*n+kk]*a[kk*l+j];
}
}
}
js=brinv(e,m);
if (js==0)
{
free(e);
free(a);
free(b);
return(js);
}
for (i=0; i<=n-1; i++)
{
for (j=0; j<=m-1; j++)
{
jj=i*m+j;
g[jj]=0.0;
for (kk=0; kk<=m-1; kk++)
{
g[jj]=g[jj]+a[i*l+kk]*e[j*m+kk];
}
}
}
for (i=0; i<=n-1; i++)
{
jj=(ii-1)*n+i;
x[jj]=0.0;
for (j=0; j<=n-1; j++)
{
x[jj]=x[jj]+f[i*n+j]*x[(ii-2)*n+j];
}
}
for (i=0; i<=m-1; i++)
{
jj=i*l;
b[jj]=y[(ii-1)*m+i];
for (j=0; j<=n-1; j++)
{
b[jj]=b[jj]-h[i*n+j]*x[(ii-1)*n+j];
}
}
for (i=0; i<=n-1; i++)
{
jj=(ii-1)*n+i;
for (j=0; j<=m-1; j++)
{
x[jj]=x[jj]+g[i*m+j]*b[j*l];
}
}
if (ii<k)
{
for (i=0; i<=n-1; i++)
{
for (j=0; j<=n-1; j++)
{
jj=i*l+j;
a[jj]=0.0;
for (kk=0; kk<=m-1; kk++)
{
a[jj]=a[jj]-g[i*m+kk]*h[kk*n+j];
}
上一篇:在SPLUS中直接连接SQLSERVER 2005
下一篇:《雷神之锤III》里求平方根的函数
|
| 相关文章： |
|
|
|
| 相关软件： |
|
| |
|
