武汉大学遥感院复试上机练习节选说课讲解

1、武汉大学遥感院复试上机练习（节选）精品文档现有一数组，其元素值见文件 test.txt，要求通过编程的方法求出这组数中值最 大的元素、值最小的元素和中值元素。结果通过编程的方法写入文本文件中输 出。#include #include using namespace std;int main()int num;ifstream filein( 机试复习 2004testtest.txt );if (!filein) cout num;int *data = new int num;for ( int i = 0; i datai; filein.close();int swap;for ( in

2、t i = 0; i num-1;i+)for ( int j = i + 1; j dataj) swap = datai;datai = dataj;dataj = swap; else continue ;/ 冒泡法排序ofstream fileout( 机试复习 2004testtest2.txt );if (!fileout) cout cannot open the file! ;return 0; fileout max: datanum-1 endl;fileout min: data0 endl;fileout medium: datanum/2 endl;fileout.c

3、lose();delete data,data= NULL;return 1;收集于网络，如有侵权请联系管理员删除精品文档编写调试程序从磁盘的数据文件“ SourceData.txt”中读入10000个数据，先理解 为 12580的矩阵，使用卷积矩阵 1 进行卷积，再理解为 100100的矩阵，使 用卷积矩阵 2 进行卷积。000 1 0011 5 112a11a120 1 a10nb11b120b1k1a21a22a2nb21b220b02kam1am2卷积矩阵amn1bl1bl 2blk1002101621210lk100 cijai l ii,j k jj biiii 1 jj 1 i

4、2 ii,j 2 jj卷积矩阵 2例：由卷积模板 B (33)对 A 矩阵(56)进行卷积得到结果矩阵 Ca11a12a13a14a15a16a21a22a23a24a25a26a31a32a33a34a35a36a41a42a43a44a45a46a51a52a53a54a55a56Ab11b12b13b21b22b23b31b32b33BC 中的元素 c33 为：c33 a22 b11 a23 b12 a24 b13 a32 b21 a33 b22 a34 b23 a42 b31 a43 b32 a44 b33#include #include using namespace std;in

5、t main()int *data = new int 10000;ifstream filein( 机试复习 2005testSourceData.txt ); if (!filein) cout cannot open the file! ;收集于网络，如有侵权请联系管理员删除精品文档return 0; for ( int i = 0; i datai; filein.close();int Mat133 = 0, -1, 0, -1, 5, -1, 0, -1, 0 ;int Mat255 = 0, 0, -1, 0, 0, 0, -1, -2, -1, 0, -1, -2, 16, -

6、2, -1,0, -1, -2, -1, 0, 0, 0, -1, 0, 0 ; int Matrix112580, MatrixNew112580 = 0 ;int Matrix2100100, MatrixNew2100100 = 0 ;int k = 0;for ( int i = 0; i 125; i+) for ( int j = 0; j 80; j+) Matrix1ij = datak;k+; k = 0;for ( int i = 0; i 100; i+) for ( int j = 0; j 100; j+) Matrix2ij = datak;k+; for ( in

7、t i = 1; i 124; i+)for ( int j = 1; j 79; j+) for ( int m = -1; m = 1; m+)for ( int n = -1; n = 1; n+)MatrixNew1ij += Matrix1i + mj + n * Mat1m + 1n + 1;for ( int i = 2; i 99; i+)for ( int j = 2; j 99; j+) for ( int m = -2; m = 2; m+)for ( int n = -2; n = 2; n+)MatrixNew2ij += Matrix2i + mj + n * Ma

8、t2m + 2n + 2;ofstream fileout( 机试复习 2005testtest2.txt ); if (!fileout) cout cannot open the file!;return 0; for ( int i = 1; i 124; i+)收集于网络，如有侵权请联系管理员删除精品文档 for ( int j = 1; j 79; j+)fileout MatrixNew1ij ; fileout endl;fileout.close();delete data,data= NULL;return 1;(1) 利用随机函数 rand()生成(0，1)区间的 100

9、个均匀分布随机数；提示： 首先利用“srand(unsigned)time(NULL) ”语句设置产生随机数的开始 点(需要包含 time.h 头文件)；产生随机数的函数为 rand()，为了将该随机 数转换为( 0，1)范围内的数，利用“rand()*1.0/RAND_MAX ”可实现此目 的。(2) 利用正态分布随机数的产生函数生成正态分布随机数。正态分布随机 数生成函数如下： y 2lnxcos2 x *0.3+1(3) 计算出 100个正态分布随机数 y 的均值和方差，计算公式如下：n n n2 2 2xxi /n ， 2 ( xi /n)-( xi /n)i 1 i 1 i 1其中，

10、自变量 x 即为 100 个均匀分布的随机数。#include #include #include #include #define PI 3.1415926using namespace std;int main()double *data = new double 100;srand( unsigned )time( NULL);收集于网络，如有侵权请联系管理员删除精品文档for ( int i = 0; i 100; i+) datai = rand()*1.0 /RAND_MA; Xdouble *data2 = new double 100;for ( int i = 0; i 10

11、0; i+) data2i = 1 + 0.3*cos(2 *PI *datai)*sqrt(-2 * log(datai); double average = 0, variance = 0;for ( int i = 0; i 100; i+) average += datai;variance += datai * datai;average /= 100;variance = variance/100 - average*average;/ 输出略delete data, data =NULL;delete data2, data2 = NULL;return 1;题目：已知某多光谱遥

12、感影像一个样区的数据。分别计算该样区数据的均值向量和协方差矩阵。计算公式如下 :1）均值向量 :mixik0M m1m2T mn2）协方差矩阵11121n式中:21222nxik 表ijNk1(xik mi )(xjk mj )0n1n2nn示第 i个波段样区的第k 个像元值， N 为样区像元总个数。multispectral.h中：#include 收集于网络，如有侵权请联系管理员删除精品文档#include using namespace std;/=/ 计算多光谱影像每个波段的均值/average: 输出每个波段的均值/num:每个波段的像元总数/band: 每个波段的波段数 /data

13、: 多光谱影像的像素值vector *data);/=void Matrix( vector *matrix, double *average, vector *data, int band, int num);Multispectral.cpp 中：#include multispectral.hvoid Average( double * average , int num, int band, vector * data) for ( int i = 0; i band; i+)average i = 0; for ( int i = 0; i num; i+) for ( int j

14、= 0; j band; j+)average j += ( double )data ji / num; void Matrix( vector * matrix , double * average , vector * data ,int band, int num) for ( int i = 0; i band; i+)for ( int j = 0; j band; j+) matrix i.push_back(0); for ( int i = 0; i band; i+)for ( int j = 0; j band; j+) for ( int k = 0; k num; k

15、+) matrix ij += (double ) data ik -averagei)*( double )data jk -average j) /num; 主函数中：#include #include #include #include multispectral.husing namespace std;int main()收集于网络，如有侵权请联系管理员删除精品文档int band, num;ifstream filein( 机试复习 2006data.txt ); / 判断略 filein band;filein num;vector *data= new vector band;

16、int count;for ( int i = 0; i band; i+)for ( int j = 0; j count;datai.push_back(count);filein.close();double *average = new double band;vector *matrix = new vector num;Average(average, num, band, data); Matrix(matrix, average, data, band, num);delete data, data =NULL;delete matrix, matrix = NULL; del

17、ete average, average = NULL;/ 输出略return 1;外包矩形的概念如下图所示：（点落在内部和边缘上）用编制的程序判断哪些点落在矩形外包内部，输出落在外包内的点号#include #include #include using namespace std; struct POINT double x;double y;收集于网络，如有侵权请联系管理员删除精品文档POINT() memset( this , 0, sizeof ( POINT); ;int main() int num;ifstream filein(机试复习 2007test.txt);if (!

18、filein) cout num;vector data;POINT point;for ( int i = 0; i point.x;filein point.y; data.push_back(point);filein.close();double xmax = data0.x, xmin = data0.x, ymax = data0.y, ymin = data0.y;for ( int i = 1; i xmax) xmax = datai.x;if (datai.xymax)ymax = datai.y;if (datai.yymin)ymin = datai.y;cout 内部

19、点有： ;for ( int i = 0; i num; i+) if (datai.x != xmax&datai.x != xmin&datai.y != ymax&datai.y != ymin)cout i + 1 ; cout endl; return 1;程序改：(将判断点放在文件中，以点名(字符)， x ，y 方式储存) struct POINTJudgedouble x;double y;char name10;POINTJudge()收集于网络，如有侵权请联系管理员删除精品文档ifstream filein2( 机试复习 2007 外包矩形 欲判断 点 .dat );if (

20、!filein2)cout cannot open the file!;return 0;vector record;while (!filein2.eof()POINTJudge jugde;filein2 jugde.xjugde.y;if (jugde.x xmin&jugde.y ymin) record.push_back(jugde);filein2.close();cout 内部点有： ;for ( int i = 0; i record.size(); i+) cout recordi.x recordi.y endl ;cout

21、endl;题目：编程统计某专业的一门课程考试成绩的 最高分、最低分、平均分 以及考试成绩分布的 均方差说明：均方差计算公式如下：2X2X2NX：考试成绩；N：成绩总个数收集于网络，如有侵权请联系管理员删除精品文档#include #include #include using namespace std;int main() int num = 0;double max = 0, min = 100, sum1 = 0, sum2 = 0;ifstream filein(新建文件夹 2008 年 score.txt );if (!filein)cout num;double *data = n

22、ew double num;for ( int i = 0; i datai;if (dataimax)max = datai;if (datai min) min = datai;sum1 += datai;sum2 += datai * datai;filein.close();double average = 0.0, jun = 0.0;average = ( double )sum1 / num;jun = sqrt( double )sum2 / num - average*average);ofstream fp2( 机试复习 2008result.txt );if (!fp2)

23、cout cannot open the file!;return 0; fp2 最高分为 max endl 最低分为 min endl;fp2 平均分为 average endl 均方差为 jun endl;fp2.close();delete data,data= NULL;return 1;试题：文件 test.txt 中存贮了一串字符串，请统计该文件中各字母出现的次数 (字母的大小写不加区分，如 a与A 算同一字母，标点符号不统计)， 并按出 现次数的多少由高到低进行排列输出 ，最后将字母及出现次数的排列写入 result.txt 文件中。#include 收集于网络，如有侵权请联系管

24、理员删除精品文档#include #include using namespace std;struct alp int num; char d; alp() memset( this , 0, sizeof ( alp ); ;int main() ifstream fp( 机试复习 2009source.txt );char data3000; /char *data =new char3000;int al26 = 0 ;int w = 0;while (!fp.eof() dataw = fp.get();w+; fp.close();for ( int j = 0; j = 0 &

25、dataj - 65 = 0 & dataj - 97 = 25) aldataj - 97+;vector alph;alp me, ma;for ( int o = 0; o 26; o+) me.d = A + o;me.num = alo;alph.push_back(me); for ( int i = 0; i 25; i+) for ( int j = i + 1; j alphj.num) ma = alphi; alphi = alphj; alphj = ma; for ( int o = 0; o 26; o+) cout alpho.d alpho.num endl;/

26、 delete data,data=NULL; 输出略 return 1;题目：猫让 10 只老鼠站成一圈，它决定每数到第 n只老鼠就吃掉，请输出最后 剩下的老鼠的位置#include #include 收集于网络，如有侵权请联系管理员删除精品文档void main() int i, j, id;int num, l;printf( 请输入老鼠的个数： n );scanf( %d, &num);printf( 请输入间隔 (范围：大于 0 且小于老鼠个数 )：n ); scanf( %d, &l);int *p;p = ( int *)malloc( sizeof (int )*num);for (i = 0; in