<?xml version="1.0" encoding="utf-8"?><rss version="2.0">
<channel>
<title><![CDATA[微雨天]]></title>
<link>http://blog.pfan.cn/weiyutian</link>
<description>编程爱好者博客</description>
<language>zh-cn</language>
			<item>
		<title><![CDATA[用c++写成的最小二乘法的源代码]]></title>
		<link>http://blog.pfan.cn/weiyutian/41873.html</link>
		<description><![CDATA[&nbsp;

#include &lt;stdio.h&gt;#include &lt;conio.h&gt;#include &lt;math.h&gt;#include &lt;process.h&gt;#define N 5//N个点#define T 3 //T次拟合#define W 1//权函数#define PRECISION 0.00001float pow_n(float a,int n){int i;if(n==0)return(1);float res=a;for(i=1;i&lt;n;i++){res*=a;}return(res);}void mutiple(float a[][N],float b[][T+1],float c[][T+1]){float res=0;int i,j,k;for(i=0;i&lt;T+1;i++)for(j=0;j&lt;T+1;j++){res=0;for(k=0;k&lt;N;k++){res+=a[i][k]*b[k][j];c[i][j]=res;}}}void matrix_trans(float a[][T+1],float b[][N]){int i,j;for(i=0;i&lt;N;i++){for(j=0;j&lt;T+1;j++){b[j][i]=a[i][j];}}}void init(float x_y[][2],int n){int i;printf("请输入%d个已知点：\n",N);for(i=0;i&lt;n;i++){printf("(x%d y%d):",i,i);scanf("%f %f",&amp;x_y[i][0],&amp;x_y[i][1]);}}void get_A(float matrix_A[][T+1],float x_y[][2],int n){int i,j;for(i=0;i&lt;N;i++){for(j=0;j&lt;T+1;j++){matrix_A[i][j]=W*pow_n(x_y[i][0],j);}}}void print_array(float array[][T+1],int n){int i,j;for(i=0;i&lt;n;i++){for(j=0;j&lt;T+1;j++){printf("%-g",array[i][]]></description>
		<author><![CDATA[与天比高]]></author>
		<pubDate>2009-03-30 20:56:00</pubDate>
		</item>
				<item>
		<title><![CDATA[用c++写成的最小二乘法的源代码]]></title>
		<link>http://blog.pfan.cn/weiyutian/41872.html</link>
		<description><![CDATA[&nbsp;

#include &lt;stdio.h&gt;#include &lt;conio.h&gt;#include &lt;math.h&gt;#include &lt;process.h&gt;#define N 5//N个点#define T 3 //T次拟合#define W 1//权函数#define PRECISION 0.00001float pow_n(float a,int n){int i;if(n==0)return(1);float res=a;for(i=1;i&lt;n;i++){res*=a;}return(res);}void mutiple(float a[][N],float b[][T+1],float c[][T+1]){float res=0;int i,j,k;for(i=0;i&lt;T+1;i++)for(j=0;j&lt;T+1;j++){res=0;for(k=0;k&lt;N;k++){res+=a[i][k]*b[k][j];c[i][j]=res;}}}void matrix_trans(float a[][T+1],float b[][N]){int i,j;for(i=0;i&lt;N;i++){for(j=0;j&lt;T+1;j++){b[j][i]=a[i][j];}}}void init(float x_y[][2],int n){int i;printf("请输入%d个已知点：\n",N);for(i=0;i&lt;n;i++){printf("(x%d y%d):",i,i);scanf("%f %f",&amp;x_y[i][0],&amp;x_y[i][1]);}}void get_A(float matrix_A[][T+1],float x_y[][2],int n){int i,j;for(i=0;i&lt;N;i++){for(j=0;j&lt;T+1;j++){matrix_A[i][j]=W*pow_n(x_y[i][0],j);}}}void print_array(float array[][T+1],int n){int i,j;for(i=0;i&lt;n;i++){for(j=0;j&lt;T+1;j++){printf("%-g",array[i][]]></description>
		<author><![CDATA[与天比高]]></author>
		<pubDate>2009-03-30 20:55:00</pubDate>
		</item>
		</channel>
</rss>