Holdon_d-CSDN博客

原创 C语言实现 IFFT 运算

参考《C》C语言实现FFT算法_杨贵安的博客-CSDN博客_c语言 fft 实现了 fft 运算；如果要实现 ifft 运算，只需要在输入和输出时取共轭就好了，故 ifft 的 C 语言实现如下（已经通过matlab 验证）。void ifft(int n, int k, long double pr[], long double pi[], long double fr[], long double fi[]){ int it, m, is, i, j, nv, l0; long dou.

2021-09-17 10:53:32 2697

原创 exit code -1073740791 (0xC0000409) 解决

问题pycharm 运行后什么都没有发生，也没有报错，只是出现了 exit code -1073740791 (0xC0000409) ，在网上寻找答案发现说是显卡内存不够，大多是采用改cpu的方法，即import osos.environ["CUDA_VISIBLE_DEVICES"] = "-1"但这个方法不适用于我，并没有解决我的问题解决多方寻求无果，在这个问答 https://ask.csdn.net/questions/874727 下面找到一个方法pycharm打开，run &

2021-06-09 00:14:12 3812 1

原创无 Anaconda 双击打开 ipynb 文件

无 Anaconda 双击打开 ipynb 文件在将 Github 一个项目下载到本地后，发现有 ipynb 文件，在 Github 上可以直接打开查询后发现是要用 Anaconda 自带的 Jupyter Notebook 打开的，可是我不想再下 Anaconda 了所以寻求了新的解决方法安装 jupyter命令行运行pip3 install jupyter找到 jupyter notebook一般可以在 python 的 Scripts 文件夹中找到，即C:\Users\User

2020-08-07 21:21:25 615

原创数值计算方法第七章—常微分方程初值问题的数值解法

常微分方程初值问题的数值解法本文参考书为马东升著《数值计算方法》引言未知函数为一元函数的微分方程叫常微分方程，讨论一阶常微分方程的初值问题{y′=f(x,y)y(x0)=y0\begin{cases}y'=f(x,y)\\y(x_0)=y_0\end{cases}{y′=f(x,y)y(x0)=y0李普希兹条件∣f(x,y1)−f(x,y2)∣≤L∣y1−y2∣|f(x,y_1)-f(x,y_2)|\le L|y_1-y_2|∣f(x,y1)−f(x,y2)∣≤L∣y

2020-07-04 15:58:12 5684 1

原创数值计算方法第六章—数值积分和数值微分

数值积分和数值微分数值积分概述基本思想用简单的函数代替复杂的 f(x)f(x)f(x) 求得复杂的积分梯形公式I=∫abf(x)dx≈(b−a)f(a)+f(b)2I=\int_a^bf(x)dx\approx(b-a)\dfrac {f(a)+f(b)}2I=∫abf(x)dx≈(b−a)2f(a)+f(b)中矩形公式I=∫abf(x)dx≈(b−a)f(a+b2)I=\int_a^bf(x)dx\approx(b-a)f(\dfrac {a+b}2)I=∫abf(

2020-06-24 19:03:01 1726

原创数值计算方法第五章—曲线拟合的最小二乘法

曲线拟合的最小二乘法本文参考书为马东升著《数值计算方法》最小二乘法与插值法的区别：只需“逼近” f(x)f(x)f(x) ，而不用满足插值原则（即经过插值节点）最小二乘原理设 y=f(x)y=f(x)y=f(x) 在区间 [ a ,b ][\ a\ ,b\ ][ a ,b ] 上有定义，求 φ(x)\varphi(x)φ(x) ，使得 φ(xi)\varphi(x_i)φ(xi) 与 f(xi)f(x_i)f(xi) 的误差

2020-06-24 19:02:04 8101

原创解决RuntimeError: The Session graph is empty. Add operations to the graph before calling run().

问题与解决今天根据 https://blog.csdn.net/Cs_hnu_scw/article/details/79695347（在此感谢博主，虽然因为时间问题根据教程出现了很多问题，但还是帮了我很多忙）这篇文章安装 Tensorflow 时，为了尝试安装是否成功，编译了博主提供的 pycharm 测试代码import tensorflow as tfimport osos.environ[‘TF_CPP_MIN_LOG_LEVEL’]=‘2’hello=tf.constant(‘1222

2020-06-19 16:48:27 4930 13

原创 win10解决 Your CPU supports instructions that this TensorFlow binary was not compiled to use: AVX AVX2

问题在使用pycharm编译一个简单的tensorflow测试文件时，出现了以下的错误Your CPU supports instructions that this TensorFlow binary was not compiled to use: AVX AVX2此外还有一大段dll错误大意是：你的电脑是支持AVX的，但是编译好的TensorFlow并不支持它搜索过程于是在网上搜索了这个问题，大多数都是通过掩耳盗铃的方法，即在顶端加入以下代码import osos.environ

2020-06-19 16:25:40 527

原创数值计算方法第四章—插值法

插值法本文参考书为马东升著《数值计算方法》代数插值插值已知某些（有限）点的函数值求其余点的函数值代数多项式插值n+1n+1n+1 个互异节点处满足插值条件 P(xi)=yi ,i=0,1,⋯ ,nP(x_i)=y_i\ ,i=0,1,\cdots,nP(xi)=yi ,i=0,1,⋯,n 的 nnn 次代数多项式 P(x)=a0+a1x+a2x2+⋯+anxnP(x)=a_0+a_1x+a_2x^2+\cdots+a_nx^nP(x)=a0+a1

2020-06-11 14:27:30 2486

原创数值计算方法第二章—非线性方程的数值解法

非线性方程的数值解法初始近似值的搜索方程的根f(x∗)=0f(x^*)=0f(x∗)=0 ，x∗x^*x∗ 为 f(x)f(x)f(x) 的根单根和重根f(x∗)f(x^*)f(x∗) 有 mmm 重根 x∗x^*x∗ 的充要条件是f(x∗)=f′(x∗)=⋯=f(m−1)(x∗)=0f(x^*)=f'(x^*)=\cdots=f^{(m-1)}(x^*)=0f(x∗)=f′(x∗)=⋯=f(m−1)(x∗)=0f(m)(x∗)≠0f^{(m)}(x^*)\ne0f(m)(

2020-06-07 13:05:52 1070

原创数值计算方法第一章—数值计算引论

数值计算引论本文参考书为马东升著《数值计算方法》误差的来源观测误差（测量误差）观测或实验所得到的参数与真值的误差模型误差（描述误差）建立的数学模型与实际事物的差距。如自由落体的数学模型 s=12gt2s=\dfrac 12gt^2s=21gt2 忽略了空气阻力截断误差（方法误差）许多数学运算是通过极限过程来定义的，而计算机只能完成有限次运算。故实际应用需要将解题方案加工成有限序列，即表现为无穷过程的截断例：数学模型是无穷级数∑k=0+∞1k!f(k)(x0)

2020-06-07 11:37:41 634

原创数值计算方法第三章—线性方程组的数值解法知识点总结

线性方程组的数值解法高斯消去法顺序高斯消去法通过初等变换消去方程组系数矩阵主对角线以下的元素，而使方程组化为等价的上三角形方程组列主元高斯消去法在每一次消元之前，将绝对值大的元素交换到主对角线的位置上高斯—若尔当消去法每次消元时利用主元将所在列其余元素全部消为0矩阵三角分解法杜利特尔分解uri=ari−∑k=1r−1lrkuki , i=r,r+1,⋯ ,nu_{ri}= a_{ri}-\sum_{k=1}^{r-1}l_

2020-05-17 22:53:40 2752

Holdon_d的博客