Big Bug Ban

数值分析–4–高斯消元法(2)

/**
+———————————————————-
* 计算矩阵—普通算法
+———————————————————-
* @access public
+———————————————————-
* @return void
+———————————————————-
*/
void cout_n() {

double x_muti;
//下列减法
double down_total=0;
//上列加法
double up_total=0;

int r;

int s;

//下列算法
for (int j = 0; j < this->cols; j++) {

x_muti=1;

r=j;
//遍历行
for (int i = 0; i < this->rows; i++) {

if(++r<=this->cols){
s=r;
}
else{
s=r–this->cols;
}
//cout<<"i:"<<i<<"s:"<<s<<endl;
// 对角相乘
x_muti=x_muti*this->DIAGRAM[i][s–1];

}

down_total += x_muti;

}
int d;

//上列算法
for (int j = 0; j < this->cols; j++) {

x_muti=1;

r=j–1;
//遍历行
d=this->rows–1;

for (int i = 0; i < this->rows; i++) {

if(++r<this->cols){
s=r;
}
else{
s=r–this->cols;
}
//cout<<"s:"<<s<<"d:"<<d<<endl;
// 对角相乘
x_muti=x_muti*this->DIAGRAM[d–][s];

}

up_total += x_muti;

}

cout<<up_total–down_total<<endl;
}
/**
+———————————————————-
* 计算矩阵—高斯消元法
+———————————————————-
* @access public
+———————————————————-
* @return void
+———————————————————-
*/
void cout_z(){

for(int j = 0; j < this->cols; j++){

int r=j+1;

for(int g=r;g<this->rows;g++)
{
double xn=this->DIAGRAM[g][j]/this->DIAGRAM[j][j];

cout<<"xn:"<<xn<<endl;

for(int s = 0; s < this->cols; s++){

this->DIAGRAM[g][s]=this->DIAGRAM[j][s]-this->DIAGRAM[g][s]/xn;
cout<<"xn:"<<(double)this->DIAGRAM[g][s]<<endl;
}
this->print_g();
}

}

}
};

}

#endif

Written by princehaku

10月 15th, 2009 at 10:25 上午

Posted in c/c++

without comments

数值分析–4–高斯消元法(1)

/*
* For: Histogiagram
* File: diagram.h
* Author: princehaku
*
* Created on 2009年10月3日, 下午3:31
*/
#ifndef diagram_H
#define diagram_H

namespace diagram {
class DIAGRAM_T {
private:
//行数
int rows;
//列数
int cols; //二维矩阵
double DIAGRAM[30][30];
public:
DIAGRAM_T(int rows, int cols) {
this->rows = rows;
this->cols = cols;
//申请矩阵存储空间
//this->DIAGRAM=(int **)(new int[rows*cols]);
//~ int **DIAGRAM= new int*[rows+3]; //第一维，
//~
//~ for( int i=0; i<rows; i++)
//~ {
//~ DIAGRAM[i] = new int[cols+3]; //分配第二维，每一行的空间。
//~ }
//~ //溢出
//~ if (!DIAGRAM) {
//~ cout << "Out of Memory";
//~ }
}

/**
+———————————————————-
* 得到单元格值
+———————————————————-
* @access public
+———————————————————-
* @param integer rows  行
* @param integer rows  列
+———————————————————-
* @return int
+———————————————————-
*/
double getValue(int row, int col) {
return this->DIAGRAM[row - 1][col - 1];
}
/**
+———————————————————-
* 设定某单元格的值
+———————————————————-
* @access public
+———————————————————-
* @param integer rows  行
* @param integer rows  列
* @param integer value 值
+———————————————————-
* @return bool
+———————————————————-
*/
bool setValue(int rows, int cols, double value) {
this->DIAGRAM[rows - 1][cols - 1] = value;
//cout<<"rows"<<rows<<"cols"<<cols<<"cg to"<<value<<endl;
return true;
}
/**
+———————————————————-
* 转置矩阵
+———————————————————-
* @access public
+———————————————————-
* @return bool
+———————————————————-
*/
bool transpose() {
double tmp;
for (int i = 0; i < this->rows; i++) {
for (int j = 0; j < this->cols; j++) {
if (i < j) {
tmp = this->DIAGRAM[i][j];
this->DIAGRAM[i][j] = this->DIAGRAM[j][i];
this->DIAGRAM[j][i] = tmp;
}
}
}
return true;
}

/**
+———————————————————-
* 列交换
+———————————————————-
* @access public
+———————————————————-
* @return bool
+———————————————————-
*/
bool l(int rowA,int rowB) {
double tmp;

for (int j = 0; j < this->cols; j++) {
tmp=this->DIAGRAM[rowA–1][j];
this->DIAGRAM[rowA–1][j]=this->DIAGRAM[rowB–1][j];
this->DIAGRAM[rowB–1][j]=tmp;
}

return true;
}
/**
+———————————————————-
* 行交换
+———————————————————-
* @access public
+———————————————————-
* @return bool
+———————————————————-
*/
bool h(int colA,int colB) {
double tmp;

for (int i = 0; i < this->rows; i++) {
tmp=this->DIAGRAM[i][colA–1];
this->DIAGRAM[i][colA–1]=this->DIAGRAM[i][colB–1];
this->DIAGRAM[i][colB–1]=tmp;
}

return true;
}
/**
+———————————————————-
* 打印矩阵
+———————————————————-
* @access public
+———————————————————-
* @return void
+———————————————————-
*/
void print_g() {

for (int i = 0; i < this->rows; i++) {
for (int j = 0; j < this->cols; j++) {
cout << this->DIAGRAM[i][j] << " ";
}
cout << endl;
}
}

Written by princehaku

10月 15th, 2009 at 10:22 上午

Posted in c/c++

without comments

我想说的..

最近心情着实不好…
不知道原因..
或许是因为事情很多..
或许是因为很难受..

几日前在公车上..
一个妇女和司机超起来..
司机差点都动手打她..
起因是那个女的没刷卡..
但她坚持说她刷了..而且证明也确实是她没刷..
那女的骂那司机..
快下车前..一个男的想提前下车
因为没进站..司机没开门..那男的下车后就破口大骂司机

司机不好做.

XX天..我和3在学校外遇到两个人
"同学..不好意思..打扰一下..我们是外地的..我和我妹妹两天没吃东西了..能给点钱不?"
我想了哈..给了他们10元..然后3说我觉得他们是骗子..
我说应该不是吧..然后我觉得至少他们没有倒在地上求人..至少还留有尊严..于是也没计较什么

第二天….回寝室的路上..遇到一个中年妇女
"同学..不好意思..打扰一下..我是外地来的..我几天没吃东西了..你看能帮我下不?"
"昨天给过了"
"那今天也再给点嘛"
我晕…

第三天….我从一个女生旁边过..
"同学..不好意思..打扰一下..我是外地来的..我几天没吃东西了..你看能帮我下不?"
"你不是有手的嘛?"

freestu….内部矛盾激化了…不知道会走向何方

人..活着..不知道什么时候会死..
什么时候会离开自己最喜欢的人..
但是也请珍惜周围的人!所有人

3….我现在就只想你…
你对我很重要..
你经常说不嫁给我..
三年后..也指不定会怎么样..
就算我会追你到广西..
但是…人这一辈子真的好短…..
所以我想有更多的时间和你在一起
好好的珍惜这段时间..
我不想等到失去时才后悔..

我只想和你在一起…

唯一的快乐..

Written by princehaku

10月 13th, 2009 at 4:51 上午

Posted in technology

without comments

想继续看XNA….

Written by princehaku

10月 7th, 2009 at 12:14 下午

Posted in technology

without comments

netbeans Qt开发环境配置

咳咳….就这个..东西…弄了我很久….

这个是官方的说明

MinGW 编译器和工具

经测试，NetBeans IDE 6.7 可以与 Windows 上的 Minimalist GNU (MinGW) 和类似 Unix 的 Minimal System (MSYS) 环境配合使用。下面显示了已测试的版本和安装说明。

软件或资源	已测试版本	描述
`gcc`	3.4.5	MinGW C 编译器
`g++`	3.4.5	MinGW C++ 编译器
`gdb`	6.8	MinGW GNU 调试器
`make`	3.79.1	MSYS `make` 实用程序请注意，不支持 MinGW make

从 mingw.org 安装 GNU 编译器、make 和 gdb 调试器：

使用具有“计算机管理员”权限的帐户登录 Windows。
从 mingw.org 上的 HOWTO Install the MinGW (GCC) Compiler Suite 中下载自解压安装程序，或为方便起见，直接从此 MinGW-5.1.4.exe 链接进行下载。
在 Windows 中运行 MinGW-5.1.4.exe 安装程序。
在 MinGW 安装程序中，接受缺省设置直至转入 “Choose Components” 页。
在 “Choose Components” 页中，选择 “g++ compiler”。gcc 编译器将自动包括在安装中，因此该编译器不是一个可选组件。
不要选择 MinGW make，因为您需要使用 MSYS make，此组件将单独下载。单击 “Next”。
如果可能，请使用缺省的 C:\MinGW 作为目标文件夹，以最大程度地减少使用其他位置的编译器可能出现的问题。
单击 “Install” 安装 MinGW 工具，然后依次单击 “Next” 和 “Finish”。
按照 MinGW wiki 中的 Environment Settings（环境设置）所述，将 C:\MinGW\bin 目录添加到 Path 变量中。
请参见 MinGW MSYS wiki 中的 MSYS 安装说明。您只需安装 MSYS 1.0 文件。无需安装 wiki 页中提到的 DTK 或核心文件。
为方便起见，可以直接从此 MSYS-1.0.10.exe 链接进行下载。
在安装 MSYS 时，请确保您使用的是具有“计算机管理员”权限的 Windows 帐户。如果使用有限权限的用户帐户安装 MSYS，以后尝试生成并运行项目时，NetBeans IDE 可能会出现问题。
运行 MSYS-1.0.10.exe 安装程序并接受缺省设置。
要安装 gdb 调试器，可直接从此 gdb-6.8-mingw-3.tar.bz2 链接进行下载。您必须具有 Windows zip 实用程序（如 gzip、7-zip 或 WinZip）才能解压缩此 gdb 归档文件。
将 gdb-6.8-mingw-3.tar.bz2 解压缩到 C:\MinGW 目录中，以便 gdb 可执行文件会被安装到 C:\MinGW\bin 目录中。
请参见验证安装以验证是否为 NetBeans IDE 正确安装了这些工具。