好的，以下是 C++ 设计的程序。请注意，我已经按照您的要求将输入文件名设置为“input.txt”，并将所有输出打印在屏幕上。程序中的注释应该能够帮助您理解代码。
```cpp
#include <iostream>
#include <fstream>
#include <vector>
#include <algorithm>
#include <cmath>
using namespace std;
// 二维坐标结构体
struct Point {
double x;
double y;
};
// 输出所有坐标
void printAllPoints(vector<Point>& points) {
for (const auto& p : points) {
cout << "(" << p.x << "," << p.y << ") ";
}
cout << endl;
}
// 按照 x 坐标排序并输出
void sortByX(vector<Point>& points) {
sort(points.begin(), points.end(), [](const Point& p1, const Point& p2) {
return p1.x < p2.x;
});
printAllPoints(points);
}
// 按照 y 坐标排序并输出
void sortByY(vector<Point>& points) {
sort(points.begin(), points.end(), [](const Point& p1, const Point& p2) {
return p1.y < p2.y;
});
printAllPoints(points);
}
// 输出所有坐标的数量
void printNumPoints(vector<Point>& points) {
cout << "Number of points: " << points.size() << endl;
}
// 输出所有坐标的质心位置
void printCentroid(vector<Point>& points) {
double sumX = 0.0;
double sumY = 0.0;
for (const auto& p : points) {
sumX += p.x;
sumY += p.y;
}
double centroidX = sumX / points.size();
double centroidY = sumY / points.size();
cout << "Centroid: (" << centroidX << "," << centroidY << ")" << endl;
}
// 最小二乘法拟合
void leastSquaresFit(vector<Point>& points) {
double sumX = 0.0;
double sumY = 0.0;
double sumXY = 0.0;
double sumX2 = 0.0;
double sumY2 = 0.0;
for (const auto& p : points) {
sumX += p.x;
sumY += p.y;
sumXY += p.x * p.y;
sumX2 += p.x * p.x;
sumY2 += p.y * p.y;
}
double n = points.size();
double slope = (n * sumXY - sumX * sumY) / (n * sumX2 - sumX * sumX);
double intercept = (sumY - slope * sumX) / n;
double e
rror = 0.0;
for (const auto& p : points) {
double residual = p.y - (slope * p.x + intercept);
e
rror += residual * residual;
}
e
rror = sqrt(e
rror / (n - 2));
cout << "Slope: " << slope << endl;
cout << "Intercept: " << intercept << endl;
cout << "E
rror: " << e
rror << endl;
}
// 输出第 n 个坐标
void printNthPoint(vector<Point>& points, int n) {
if (n >= 1 && n <= points.size()) {
cout << "(" << points[n-1].x << "," << points[n-1].y << ")" << endl;
} else {
cout << "Invalid index." << endl;
}
}
// 沿 x 轴进行镜像变换并输出
void
mirrorX(vector<Point>& points) {
for (auto& p : points) {
p.y = -p.y;
}
printAllPoints(points);
}
// 沿 y 轴进行镜像变换并输出
void
mirrorY(vector<Point>& points) {
for (auto& p : points) {
p.x = -p.x;
}
printAllPoints(points);
}
// 绕原点逆时针旋转 a 度并输出
void rotate(vector<Point>& points, double a) {
double radians = a * M_PI / 180.0;
for (auto& p : points) {
double newX = p.x * cos(radians) - p.y * sin(radians);
double newY = p.x * sin(radians) + p.y * cos(radians);
p.x = newX;
p.y = newY;
}
printAllPoints(points);
}
// 沿 x 轴和 y 轴分别平移 x0 和 y0 并输出
void translate(vector<Point>& points, double x0, double y0) {
for (auto& p : points) {
p.x += x0;
p.y += y0;
}
printAllPoints(points);
}
int main() {
// 读取数据
ifstream inputFile("input.txt");
if (!inputFile.is_open()) {
cout << "E
rror: unable to open input file." << endl;
return 0;
}
vector<Point> points;
double x, y;
while (inputFile >> x >> y) {
Point p = { x, y };
points.push_back(p);
}
inputFile.close();
// 控制程序
char input;
while (true) {
cout << "Enter a command: ";
cin >> input;
switch (input) {
case 's':
printAllPoints(points);
break;
case 'x':
sortByX(points);
break;
case 'y':
sortByY(points);
break;
case 'n':
printNumPoints(points);
break;
case 'c':
printCentroid(points);
break;
case 'l':
leastSquaresFit(points);
break;
case 'i':
int n;
cout << "Enter an index: ";
cin >> n;
printNthPoint(points, n);
break;
case 'v':
mirrorX(points);
break;
case 'h':
mirrorY(points);
break;
case 'r':
double a;
cout << "Enter an angle of rotation: ";
cin >> a;
rotate(points, a);
break;
case 'm':
double x0, y0;
cout << "Enter x and y translations: ";
cin >> x0 >> y0;
translate(points, x0, y0);
break;
case 'e':
return 0;
default:
cout << "Invalid command." << endl;
break;
}
}
return 0;
}
```
2747
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
