MATLAB 求解常微分方程的函数是dsolve()。应用此函数可以求得常微分方程(组)的通解，以及给定边界条件(或初始条件)后的特解。S=dsolve(eqn)S=dsolve(eqn,cond)S=dsolve(eqn,cond,Name,Value)eqn1,eqn2,…为给定的常微分方程(组)。② Value为给定的常微分方程(组)的指定符号自变量， …
编程文章
在农业中，通常希望获取不同土地的面积。虽然获取这些土地的面积操作相对容易，但是却涉及高额的费用。另外，如果对于不规则形状的土地，测量土地面积的大小就变得相对困难。幸运的是，有大量以卫星图像的形式公开提供的农场土地数据。如下图所示就是得到的土地的图片。根据这样的图片，应用图像处理算法，就可以计算得到每块土地的面积。我们可以使用导数滤波器，因为它似乎与人类如何区 …
IFFT是傅里叶逆变换的一种实现，它将时域的离散信号转换为频域的离散信号。具体而言，ifft将复数向量作为输入，然后计算其离散傅里叶变换的逆变换。ifft的输出是一个复数向量，其大小与输入向量相同。以下是使用C语言实现ifft函数的示例代码：#include <stdio.h>#include <stdlib.h>#include & …
相信各位学习自动化或者电气的同学都遇到过用matlab来绘制传递函数的波特图的问题，下面我给大家介绍一种绘制波特图的简单方法。使用Matlab来绘制波特图的过程中，一般需要用到tf()函数。具体计算程序如下则在窗口直接输入程序a = [1,2];b = [2,3,1];x= tf(a,b);bode(x);即可得到传递函数的幅频特性曲线有了幅频特性曲线即可分 …
今天主要是讲解MATLAB的牛顿法求多元函数的极值程序加实例。实例1求f(x,y)= sin(x^2+y^2)*exp(-0.1*(x^2+y^2+x*y+2*x))，在-2<=x<=2，-2<=y<=2上的极值点和极值。主程序clc;clear all;close all;syms x y;%定义函数变量 x yf = sin(x^ …
实例程序clc;clear all;close all;%计算二元函数的极值点 并进行判断syms x y%定义二元变量 x yz = (6*x-x^2)*(4*y-y^2);%定义二元变量函数f1 = simplify(diff(z,x));%求z对x的一阶偏导f2 = simplify(diff(z,y));%求z对y的一阶偏导%求f1 = 0 f2 = …
数值积分的实现：（一）比较quad与quadl函数的差别基于自适应辛普森方法：[I,n]=quad(filename,a,b,tol,trace)基于自适应Gauss-Lobatto方法：[I,n]=quadl(filename,a,b,tol,trace)其中，filename是被积函数名；a和b分别是定积分的下限和上限，积分限[a,b]必须是有限的，不能 …
在C语言中实现快速傅里叶变换（FFT）可以使用Cooley-Tukey算法，这是一种递归的算法。#include <math.h>#include <complex.h>#define PI 3.14159265358979323846void fft(double complex *x, int N){if (N == 1)retu …
热烈庆祝《大家一起晒显卡》办到第10期，感谢大家的热情参与，小菜我在这里给大家深深地鞠一个躬。这一期我们为大家带来的是一款奠定NVIDIA显卡王朝基石的显卡，是NVIDIA在当年登顶性能之王的一块卡，一些老玩家一定猜到了，没错，这块卡就是大名鼎鼎的NVIDIA Riva TNT2，由网友“yan220550@sina.cn”友情提供。华硕 AGP-V3800 …
3月3日，在第54个学雷锋活动日前夕，兰州铁路局迎水桥机务段“雷锋号”机车延续至第四代，启用和谐号D3D型0035号机车，北部战区某部“雷锋班”副班长张阳（左）和和谐号“雷锋号”机车组司机长（右）共同为新式机车揭幕。“雷锋号”机车组成立于1989年5月，机车为前进型2690号蒸汽机车。1990年2月5日，该机车组被原银川铁路分局团委正式命名为“雷锋号”,成为 …