随笔分类 -  数值算法

摘要：数学和物理中常见的误区 阅读全文

摘要：在工程物理的研究中，很多现象（例如：波动方程）的定量描述需要用到二阶微分方程，这里把相关的知识整理一下，以便随时使用，内容都是高数教材[1][2]上搬过来的，谈不上原创，但是自己整理一遍确实可以加深印象。 对于微分方程，如果在方程中只出现一个未知函数的导数，那么就是常微分方程，反之就是偏微分方程。如 阅读全文

摘要：在matlab一般可以通过函数的输入和输出接口来传递变量，但是有的时候，我们因为函数已经写好，并在比较多的脚本中被调用，此时再修改输入接口的工作量比较大，如果可以通过其它方式传递变量，就可以省去很多的工作量。初步想了一下，可以通过两种形式实现： 1. 在主函数中定义全局变量 但是这里需要注意的是全局 阅读全文

摘要：在对信号进行傅立叶变换时得到的频谱，会有频率为负数的部分，对于一维FFT如果画出正负频率的频谱，并以幅度的绝对值表示，可以发现在频谱上正负频率是对称分布的，这是为什么呢？对于二维傅立叶变换也有相同的规律吗？ 1. 对于一维傅立叶变换，系数是对称的，如果从数学理论上做些分析，可以发现对于任意信号g(t 阅读全文

摘要：在信号处理中常常需要用到曲线拟合，这里介绍一下利用最小二乘拟合一般曲线的方法，并对滤掉信号中白噪声的方法作些介绍。 为了测试拟合算法的好坏，先模拟出一个信号作为检验算法的例子： 1. 用白噪声产生模拟信号： 对于理论信号y=y(x)，一般可用rand(size(x))和randn(size(x))生 阅读全文

摘要：下面的matlab代码通过窗口傅立叶变换计算了信号的功率谱和互功率谱(power spectrum for time frequency analysis) 阅读全文

摘要：这里为MHD（磁流体力学）模拟做一些准备，可能现在学习物理的人从某种程度上也得变成程序猿吧。MHD模拟面临的问题是求解一个三维的偏微分方程组，其中涉及的流体元格点非常多。所以希望先从简单的一维模型出发，做一些简单的数值解。希望可以在编写程序的过程中用到面向对象的思想，可能利用上NAVIDIA显卡的G 阅读全文