计算方法实验
文章平均质量分 83
千里澄江
冷静地思考,热情地生活。
知乎主页:https://www.zhihu.com/people/95-9-17-28-56
github主页:https://github.com/gmh1627
展开
-
计算方法实验8:快速傅里叶变换与快速傅里叶逆变换实现对给定函数的 Fourier 分析及重建
通过快速傅里叶变换与快速傅里叶逆变换实现对给定函数的 Fourier 分析,t∈01,函数f以及划分数nf1t0.7sin2π×2tsin2π×5tn2427f2t0.7sin2π×2tsin2π×5t0.3×randomt27fnggfng∣gi∣n∣gi∣f1f1f1knf2f2f2k。原创 2024-05-20 14:49:49 · 1130 阅读 · 0 评论 -
计算方法实验7:实现三次样条插值算法
point.txt文件中包含了21个压铁的位置信息。原创 2024-05-12 12:24:52 · 997 阅读 · 0 评论 -
计算方法实验9:Romberg积分求解速度、位移
现在要用数值方法求∫abf(x) dx\int_{a}^{b} f(x) \, dx∫abf(x)dx,设h=b−anh=\frac{b-a}{n}h=nb−a,已知:复化梯形积分Tn(f)=h[12f(a)+∑i=1n−1f(a+ih)+12f(b)]T_{n}\left(f\right)=h\left[\frac{1}{2}f\left(a\right)+\sum_{i=1}^{n-1}f\left(a+ih\right)+\frac{1}{2}f\left(b\right)\right]Tn(f原创 2024-05-09 21:12:54 · 833 阅读 · 0 评论 -
计算方法实验2(补充):列主元消元法解线性方程组
C++源代码#include<bits/stdc++.h>using namespace std;// 列主元消去法求解线性方程组vector<long double> Column_Elimination(vector<vector<long double>> A, vector<long double> b);int main(){ vector<vector<long double>> A(3,原创 2024-05-08 22:45:27 · 506 阅读 · 0 评论 -
计算方法实验4:Jacobi方法求矩阵特征值(C++)
【代码】计算方法实验4:Jacobi方法求矩阵特征值。原创 2024-04-16 15:32:27 · 590 阅读 · 0 评论 -
计算方法实验6:对鸢尾花数据集进行主成分分析(PCA)并可视化
iris数据集包含150条数据,从iris.txt读取,每条数据有4个属性值和一个标签(标签取值为0,1,2)。要求对这150个4维数据进行PCA,可视化展示这些数据在前两个主方向上的分布,其中不同标签的数据需用不同的颜色或形状加以区分。原创 2024-04-14 15:12:15 · 1437 阅读 · 0 评论 -
计算方法实验5:C++实现矩阵的奇异值分解
生成一个4 × 3的随机矩阵A,应用Jacobi方法求解矩阵AAT的特征值,计算矩阵A的SVD分解。要求A的每个元素均为[0;1]区间内的随机数。原创 2024-04-14 12:53:59 · 1388 阅读 · 0 评论 -
计算方法实验3:反幂法求按模最小特征值及特征向量
ALU其中L为单位下三角阵,U为上三角阵.则Axb可化为LUxLyb只需解方程组LybUxy。原创 2024-03-31 16:04:18 · 941 阅读 · 0 评论 -
计算方法实验2:列主元消元法和Gauss-Seidel迭代法解线性方程组
即需解线性方程组Ayb,其中A99×99−2ϵhϵ0⋮0ϵh−2ϵhϵ⋮00ϵh−2ϵh⋱⋯⋯⋯⋯⋱ϵ000⋮−2ϵhy00y1001yy1y2y3⋮y99bah2ah2⋮ah2ah2−ϵ−h。原创 2024-03-24 12:39:24 · 912 阅读 · 0 评论 -
计算方法实验1:圆形镜面成像问题
设Tcosθsinθ,则有PTQTPx−cosθ2sin2θQx−cosθ2Qy−sinθ2PTQTPx2−2Pxcosθ1Qx2Qy21−2Qxcosθ−2Qysinθ由费马原理,光线沿PTQT最短的路径传播,因此只需对上式求导求极小值点。关于θ求导得Px2−2Pxcosθ1Pxs。原创 2024-03-17 16:01:15 · 1041 阅读 · 0 评论