C++博客-NoNews-随笔分类-改进http://www.cppblog.com/hannosogno/category/17679.htmlzh-cnThu, 08 Sep 2011 15:52:15 GMTThu, 08 Sep 2011 15:52:15 GMT60.Net多线程总结http://www.cppblog.com/hannosogno/archive/2011/09/08/155385.htmlNoNewsNoNewsThu, 08 Sep 2011 14:46:00 GMThttp://www.cppblog.com/hannosogno/archive/2011/09/08/155385.htmlhttp://www.cppblog.com/hannosogno/comments/155385.htmlhttp://www.cppblog.com/hannosogno/archive/2011/09/08/155385.html#Feedback0http://www.cppblog.com/hannosogno/comments/commentRss/155385.htmlhttp://www.cppblog.com/hannosogno/services/trackbacks/155385.html阅读全文

NoNews 2011-09-08 22:46 发表评论
]]>C++封装好的曲线拟合http://www.cppblog.com/hannosogno/archive/2011/09/08/155387.htmlNoNewsNoNewsThu, 08 Sep 2011 14:46:00 GMThttp://www.cppblog.com/hannosogno/archive/2011/09/08/155387.htmlhttp://www.cppblog.com/hannosogno/comments/155387.htmlhttp://www.cppblog.com/hannosogno/archive/2011/09/08/155387.html#Feedback0http://www.cppblog.com/hannosogno/comments/commentRss/155387.htmlhttp://www.cppblog.com/hannosogno/services/trackbacks/155387.html用C++封装好的曲线拟合算法，现在已经封装的非常好用而具有很大的伸缩性了！
封装好的类用起来很方便，用一个vector初始化一下，直接就可以得到拟合的曲线的系数了。
//file : curve_fitting.h

curve_fitting.h

#pragma once
#include <math.h> 
#include <iostream>
#include <vector>
#include <algorithm>
#include <cassert>
#include <exception>
using namespace std;
namespace curve_fit
{
struct cf_Point
{
double x;
double y;
cf_Point(double _x, double _y) : x(_x), y(_y)
{}
};
class cf_Matrix2D
{
public:
explicit cf_Matrix2D(int _height, int _width) : width(_width), height(_height)
{
for (int i = 0; i < _height; ++i)
{
double* d = new double[_width]();
data.push_back(d);
}
}
~cf_Matrix2D()
{
for (int i = 0; i < height; ++i)
{
delete[] data[i];
}
}
double& at(int x, int y)
{
assert(x < height && y < width && "out of range in Matrix2D");
return data[x][y];
}
double* operator[](int i)
{
return data[i];
}
const double* const operator[](int i) const
{
return data[i];
}
int get_width() const
{
return width;
}
int get_height() const
{
return height;
}
private:
vector<double*> data;
int width;
int height;
};
class curve_fitting
{
public:
explicit curve_fitting(const vector<cf_Point>& init_points, int _fit_scale = 3) 
: points(init_points)
, fit_scale(_fit_scale)
, W(1), PRECISION(0.00001)
{
number_of_points = init_points.size();
assert(number_of_points > 0 && "要拟合的原始点个数不能为0");
}
vector<double> get_coefficient_vector()
{
if (coefficient.size() == 0)
{
calc_coefficient();
}
assert(coefficient.size() == fit_scale + 1 && "返回的系数个数应该等于阶数+1但实际上不等于，这里发生了未知错误");
return coefficient;
}
protected:
private:
void   calc_coefficient();
void   mutiple(cf_Matrix2D& a, cf_Matrix2D& b, cf_Matrix2D& c) ;
void   matrix_trans(cf_Matrix2D& a, cf_Matrix2D& b) ;
void   get_A(cf_Matrix2D& matrix_A) ;
void   convert(cf_Matrix2D& argu,int   n) ;
void   compute(cf_Matrix2D& argu,int n, vector<double>& root);
void   get_y(cf_Matrix2D& trans_A,vector<double>& y,int n);
void   cons_formula(cf_Matrix2D& coef_A, vector<double>& y, cf_Matrix2D& coef_form);
int number_of_points; //用于曲线拟合的已知点的个数
int fit_scale; //拟合曲线的阶数
int W; //权函数
double PRECISION; //精度
vector<cf_Point> points; //已知的顶点
vector<double> coefficient; //保存得到的系数
};
}
//file : curve_fitting.cpp

curve_fitting.cpp

#include "stdafx.h"
#include "curve_fitting.h"
namespace curve_fit
{
void   curve_fitting::mutiple(cf_Matrix2D& a, cf_Matrix2D& b, cf_Matrix2D& c) 
{ 
double   res=0; 
int   i,j,k; 
for(i=0;i <fit_scale+1;i++) 
for(j=0;j <fit_scale+1;j++) 
{ 
res=0; 
for(k=0;k <number_of_points;k++) 
{ 
res += a.at(i, k) * b.at(k, j);
c.at(i, j)=res; 
} 
} 
} 
void   curve_fitting::matrix_trans(cf_Matrix2D& a, cf_Matrix2D& b) 
{ 
int   i,j; 
for(i=0;i <number_of_points;i++) 
{ 
for(j=0;j <fit_scale+1;j++) 
{ 
b.at(j, i) = a.at(i, j);
} 
} 
} 
void   curve_fitting::get_A(cf_Matrix2D& matrix_A) 
{ 
int   i,j; 
for(i=0;i <number_of_points;i++) 
{ 
for(j=0;j <fit_scale+1;j++) 
{ 
matrix_A.at(i, j) = W * pow(points[i].x, j); 
} 
} 
} 
void   curve_fitting::convert(cf_Matrix2D& argu,int n) 
{ 
int   i,j,k,p,t; 
double   rate,temp; 
for(i=1;i <n;i++) 
{ 
for(j=i;j <n;j++) 
{ 
if(argu.at(i-1, i-1) == 0) 
{ 
for(p=i;p <n;p++) 
{ 
if(argu.at(p, i-1)!=0) 
break; 
} 
if(p==n) 
{ 
throw exception("方程组无解! curve_fitting::convert");
} 
for(t=0;t <n+1;t++) 
{ 
temp=argu.at(i-1, t); 
argu.at(i-1, t)=argu.at(p, t); 
argu.at(p, t)=temp; 
} 
} 
rate=argu.at(j, i-1)/argu.at(i-1, i-1); 
for(k=i-1;k <n+1;k++) 
{ 
argu.at(j, k)-=argu.at(i-1, k)*rate; 
if(fabs(argu.at(j, k)) <=PRECISION) 
argu.at(j, k)=0; 
} 
} 
} 
} 
void   curve_fitting::compute(cf_Matrix2D& argu,int n, vector<double>& root) 
{ 
int   i,j; 
double   temp; 
root.clear();
for (int i = 0; i < fit_scale + 1; ++i)
{
root.push_back(0);
}
for(i=n-1;i>=0;i--) 
{ 
temp=argu.at(i, n); 
for(j=n-1;j> i;j--) 
{ 
temp-=argu.at(i, j) * root[j]; 
} 
root[i]=temp/argu.at(i, i); 
if (root[i] < PRECISION)
{
root[i] = 0;
}
} 
} 
void   curve_fitting::get_y(cf_Matrix2D& trans_A,vector<double>& y,int n) 
{ 
int   i,j; 
double   temp; 
for(i=0;i <n;i++) 
{ 
temp=0; 
for(j=0;j <number_of_points;j++) 
{ 
temp+=trans_A.at(i, j) * points[j].y; 
} 
y[i]=temp; 
} 
} 
void   curve_fitting::cons_formula(cf_Matrix2D& coef_A, vector<double>& y, cf_Matrix2D& coef_form) 
{ 
int   i,j; 
for(i=0;i <fit_scale+1;i++) 
{ 
for(j=0;j <fit_scale+2;j++) 
{ 
if(j==fit_scale+1) 
coef_form.at(i, j) = y[i]; 
else 
coef_form.at(i, j)=coef_A.at(i, j); 
} 
} 
} 
void curve_fitting::calc_coefficient()
{
vector<double> y(fit_scale + 1);
cf_Matrix2D matrix_A(number_of_points, fit_scale+1);
cf_Matrix2D trans_A(fit_scale + 1, number_of_points);
cf_Matrix2D coef_A(fit_scale + 1, fit_scale + 1);
cf_Matrix2D coef_formu(fit_scale + 1, fit_scale + 2);
get_A(matrix_A); 
matrix_trans(matrix_A , trans_A);
mutiple(trans_A , matrix_A , coef_A); 
get_y(trans_A, y, fit_scale+1); 
cons_formula(coef_A , y , coef_formu); 
convert(coef_formu , fit_scale+1); 
compute(coef_formu , fit_scale+1 , coefficient); 
}
}
//file : main.cpp

main.cpp

#include "stdafx.h"
#include <iostream>
#include <vector>
#include "curve_fitting.h"
using namespace std;
using namespace curve_fit;
void   main() 
{ 
vector< cf_Point > v;
v.push_back(cf_Point(-2, -4));
v.push_back(cf_Point(-1, 0));
v.push_back(cf_Point(0, 0));
v.push_back(cf_Point(1, 2));
v.push_back(cf_Point(2, 12));
v.push_back(cf_Point(3, 36));
curve_fitting cf(v);
vector<double> coefficient = cf.get_coefficient_vector();
cout << "f(x) = " << coefficient[0] << " * X";
for (size_t i = 1; i < coefficient.size(); ++i)
{
cout << "  +  " << coefficient[i] << " * X^" << i << " ";
}
cout << endl;
}


NoNews 2011-09-08 22:46 发表评论
]]>救急代码：简易的Socket通信http://www.cppblog.com/hannosogno/archive/2011/08/01/152227.htmlNoNewsNoNewsMon, 01 Aug 2011 10:11:00 GMThttp://www.cppblog.com/hannosogno/archive/2011/08/01/152227.htmlhttp://www.cppblog.com/hannosogno/comments/152227.htmlhttp://www.cppblog.com/hannosogno/archive/2011/08/01/152227.html#Feedback0http://www.cppblog.com/hannosogno/comments/commentRss/152227.htmlhttp://www.cppblog.com/hannosogno/services/trackbacks/152227.html以下的四段代码完成了使用简单的TCP/UDP协议来进行通信的全部所需！


  1 #include <string>
  2 #include<stdlib.h>
  3 #include<stdio.h>
  4 #include<conio.h>
  5 #include<winsock2.h>
  6 #include<windows.h>
  7 #pragma comment(lib,"ws2_32.lib")
  8 using namespace std;
  9 
 10 void sbBroadcast(void)
 11 {
 12     SOCKET listener;
 13     sockaddr_in listen_adress;
 14     WSADATA wsdata;
 15     BOOL optval;
 16     //启动SOCKET库，版本为2.0
 17     WSAStartup(0x0202,&wsdata);
 18     optval=TRUE;
 19 
 20     //然后赋值给两个地址，一个用来绑定套接字，一个用来从网络上的广播地址接收消息；
 21     listen_adress.sin_family=AF_INET;
 22     listen_adress.sin_addr.s_addr=INADDR_ANY;
 23     listen_adress.sin_port=htons(9290);
 24     listener=socket(AF_INET,SOCK_DGRAM, 0);
 25     //设置该套接字为广播类型，
 26     //setsockopt(listener,SOL_SOCKET,SO_BROADCAST,(char FAR *)&optval,sizeof(optval));
 27     bind(listener,(sockaddr *)&listen_adress,sizeof(sockaddr_in));
 28     listen(listener, 5);
 29 
 30     char buf[256];
 31     while(1)
 32     {//从广播地址接收消息，注意用来绑定的地址和接收消息的地址是不一样的
 33         sockaddr_in from;
 34         int fromlength = sizeof(SOCKADDR);
 35         int rev = recvfrom(listener,buf, sizeof(buf),0,(struct sockaddr FAR *)&from,(int FAR *)&fromlength);
 36         if (rev != SOCKET_ERROR)
 37         {
 38             static int num = 0;
 39             printf("%s(%d) : %s\n",inet_ntoa(from.sin_addr), num, buf); //在此处输出检查到的IP,对吗????
 40             ZeroMemory(buf, sizeof(buf));
 41             num++;
 42         }
 43     }
 44 }
 45 
 46 void fbBroadcast(void)
 47 {
 48     SOCKET s;
 49     sockaddr_in addrto;
 50     WSADATA wsdata;
 51     BOOL bsocket;
 52     char *smsg="127.0.0.1&9290&255.255.255.255&9290&1&HelloMessage";
 53     WSAStartup(0x0202,&wsdata);
 54 
 55     s=WSASocket(AF_INET,SOCK_DGRAM, IPPROTO_UDP,NULL,0,WSA_FLAG_OVERLAPPED);
 56     bsocket=TRUE;
 57     //然后赋值给地址，用来从网络上的广播地址接收消息；
 58     addrto.sin_family=AF_INET;
 59     addrto.sin_addr.s_addr=INADDR_BROADCAST;
 60     addrto.sin_port=htons(9290);
 61     bool opt=true;
 62     //设置该套接字为广播类型，
 63     setsockopt(s,SOL_SOCKET,SO_BROADCAST,(char FAR *)&opt,sizeof(opt));
 64     int nlen=sizeof(addrto);
 65     while(1)
 66     {
 67         Sleep(10);
 68         static int num = 0;
 69         //从广播地址发送消息
 70         int ret=sendto(s,smsg,256,0,(sockaddr*)&addrto,nlen);
 71         if(ret==SOCKET_ERROR)
 72         {
 73             printf("GetLastError: %d \n",WSAGetLastError());
 74         }
 75         else
 76         {  
 77             printf("ok, num = %d \n", num);    
 78         }
 79         num++;
 80         Sleep(1000);
 81     }
 82 }
 83 
 84 void sbTcp(void)
 85 {
 86     SOCKET listener;
 87     sockaddr_in listen_adress;
 88     WSADATA wsdata;
 89     //启动SOCKET库，版本为2.0
 90     WSAStartup(0x0202,&wsdata);
 91 
 92     //然后赋值给两个地址，一个用来绑定套接字，一个用来从网络上的广播地址接收消息；
 93     listen_adress.sin_family=AF_INET;
 94     listen_adress.sin_addr.s_addr=INADDR_ANY;
 95     listen_adress.sin_port=htons(9290);
 96     listener=socket(AF_INET,SOCK_STREAM, 0);
 97 
 98     bind(listener,(sockaddr *)&listen_adress,sizeof(sockaddr_in));
 99     listen(listener, 5);
100 
101     char buf[256];
102     while(1)
103     {//从广播地址接收消息，注意用来绑定的地址和接收消息的地址是不一样的
104         sockaddr_in from;
105         int fromlength = sizeof(SOCKADDR);
106         SOCKET temp = accept(listener, (sockaddr *)&from, &fromlength);
107         
108         recv(temp, buf, 256, 0);
109         printf("%s : %s\n",inet_ntoa(from.sin_addr), buf); //在此处输出检查到的IP,对吗????
110 
111         
112         string msg = "ok~~~~~~~~~~";
113 
114         printf("preper to send\n");
115 
116         send(temp, msg.c_str(), msg.size() + 2, 0);
117         ZeroMemory(buf,256);
118         closesocket(temp);
119     }
120 }
121 
122 
123 void FbTcp(void)
124 {
125     sockaddr_in server_addres;
126     WSADATA wsdata;
127     char *smsg="Hello Socket!";
128     WSAStartup(0x0202,&wsdata);
129 
130     //然后赋值给地址，用来从网络上的广播地址接收消息；
131     server_addres.sin_family=AF_INET;
132     server_addres.sin_addr.s_addr = inet_addr("127.0.0.1");
133     server_addres.sin_port=htons(9290);
134     
135     int nlen=sizeof(server_addres);
136     while(1)
137     {
138         SOCKET s = socket(AF_INET,SOCK_STREAM, 0);
139         int ret = connect(s, (LPSOCKADDR)&server_addres, sizeof(server_addres));
140         
141         static int num = 0;
142         char buf[256];
143         //从广播地址发送消息
144         ret=send(s,smsg,256,0);
145         ret = recv(s, buf, 256, 0);
146         printf("%s  @  %l", buf, &buf);
147         if(ret==SOCKET_ERROR)
148         {
149             printf("GetLastError: %d \n",WSAGetLastError());
150         }
151         else
152         {  
153             printf("ok, num = %d \n", num);    
154         }
155         num++;
156         closesocket(s);
157 
158         Sleep(100);
159     }
160 }
161 
162 
163 
164 #include <string>
165 #include <iostream>
166 using namespace std;
167 
168 
169 
170 int main()
171 {
172     char c = getchar();
173     if (c == 's')
174     {
175         sbBroadcast();
176     }
177     else if(c == 'f')
178     {
179         fbBroadcast();
180     }
181     else if (c == 'z')
182     {
183         sbTcp();
184     }
185     else
186     {
187         FbTcp();
188     }
189 }


NoNews 2011-08-01 18:11 发表评论
]]>