chenwei的博客

实现：%%Question 1 part(b)clear all;clc;t0 = 0;x0 = 1/2;dt = 0.1;tf = 1;t_range = t0:dt:tf;x_EU = zeros(1,length(t_range));x_EU(1)= x0;x_RK = zeros(1,length(t_range));x_RK(1)= x0;for k = 1:length(t_range) - 1 x_EU(k+1) = euler_scheme(x_EU...

代码前部分是求解偏微分方程，后面一部分是绘图。在绘图的过程中，主要难点有两个：第一， 如何使用python或MATLAB绘制等高线以及热力图，并把它们叠加在一张图上。第二，如何使用python或MATLAB绘制二维的向量图。MATLAB实现代码：% Solution of 2D Stokes and continuity equations with finite differences% on a regular grid using stream function - vorticit

例1clear,clc%% 数据准备x1 = [7 1 11 11 7 11 3 1 2 21 1 11 10]';x2 = [26,29,56,31,52,55,71,31,54,47,40,66,68]';x3 = [6 15,8,8,6,9,17,22,18,4,23,9,8]';x4 = [60,52,20,47,33,22,6,44,22,26,34,12,12]';y = [78.5 74.3 104.3 87.6 95.9 109.2 102.7 72.5 93.1 115.9

例1例2例3

加入高斯白噪声并作平滑处理clear,clc%% (1) 读取数据A = imread('LenaRGB.jpg');%% (2)在数据中加入高斯白噪声作平滑处理% 加入高斯白噪声B = double(A)/255;[m,n,p] = size(B); % 大小white_noise = 0 + 0.1*randn(m,n); % 白噪声C = B + white_noise; % 加白噪声D = C*255; %将像素范围扩大至0--255D2 = uint8(D);

现在很多工作都转移到了线上，有些情况下需要将签名插入word，但是我们有时候是没法打印后再签名的，因此我就写了这么个程序。一、读取图片%% 读取图片pic = imread('C:\Users\cw817615\Desktop\签名.jpg');threshold = graythresh(pic);%确定二值化阈值pic_bw = im2bw(pic,0.55);%对图像二值化，黑色的为0，白色的为1imshow(pic_bw) % 展示处理前的二值化图片pic_backup = p

文章目录数值解符号解数值解问题Matlab实现main.m 主函数clear,clcx0 = 0; % xminx1 = 1; % xmaxh = 0.1; % 步长y0 = 1/2; % y0为初始值% 欧拉法[xy{1}(:,1), xy{1}(:,2)] = forwardEuler(x0, y0, x1, h);% 改进欧拉法[xy{2}...

前言：绘制连续型的分布函数很容易，直接根据分布函数计算函数值即可。但是对于离散型随机变量而言，没有已知的分布函数，只能使用经验分布函数或者说是累积分布函数进行近似。下文以离散分布中比较经典的二项分布、泊松分布以及几何分布为例绘制它们的经验分布函数。二项分布二项分布的理论部分可以参照二项分布百度百科Len = 100000; %% 二项分布Bin(n,p)N = 50; % 二项分...

本篇文章的代码实现了以下三大方面的功能：一、计算级比和光滑比并做级比检验；二、序列的灰色预测；三、精度检验，主要做了以下内容：①相对残差Q检验（MAPE）；②关联度检验；③方差比C检验；④小误差概率P检验1 灰色预测模型本部分内容参照了王佳佳的硕士学位论文《京津冀地区高速铁路对区域经济影响研究》，我觉得这篇论文讲灰色预测讲得蛮详细的，如果想学灰色预测的可以看看。2...

前言： 找到的第一篇文献，里面的密度函数出现了明显问题，在优快云里面找也经常是只找到一个或者两个分布的估计，但是有时候得四个分布一起用，找来找去难免会有些麻烦，找到一篇文献：邓泽怀的硕士学位论文《基于实测数据的海杂波统计建模》，里面详细地讲解了上面提到的四种分布参数的矩估计、极大似然估计，并作了拟合优度检验。实践证明这里面讲到的关于参数估计的内容非常准确。文章目录1 背景知识2 典型分布模型2...

前言： 折腾了老半天，直接在优快云里面搜找不到代码，百度刚开始也就找到一个bnb20工具箱，可那已经是2009年的玩意了。不过还是很幸运，在我绝望到想用遗传算法求解的时候，终于找到了一个叫Yalmip的工具箱，讲道理，功能强大且好用。1 Yalmip工具箱的下载与安装Yalmip的下载(建议在我给的这个链接里下载，官网下载的速度实在是emmmm)Yalmip的安装2 Yalmip的使用...

前言： 第一次听说Box-Cox变换，发现原来这玩意这么厉害，值得写一写记录一下。以下内容参考了栗然等人写的论文《基于Box-Cox变换的风电场短期风速预测模型》1 问题背景风力发电是解决当前突出的能源和环境问题的有效手段，因而得到了普遍重视和大规模开发利用，是目前世界上增长速度最快的能源开发形式。尤其近十几年来其发电成本大幅下降，已具备与传统常规电源竞争的潜力。目前，国内外对于风力发电的...

**前言：**大学期间只学习过《常微分方程》，没想到有些学校竟然还学《时滞微分方程》，于是找到一本由内藤敏机(日本)等著，马万彪等译的《时滞微分方程——泛函数微分方程引论》(有需要的可以私聊，优快云貌似上传不了书籍，说侵权emmm)，看着头秃，不过受到不少启发，尤其是对Logistic方程的改进，真真是长见识了。没找到有人用欧拉法解一阶时滞微分方程的，于是一不做二不休便用MATLAB实现了一下下...

文章目录1.背景2. 单种运输工具带时间窗2.1 带时间窗车辆路径问题的描述2.2 遗传算法的求解2.2.1 编码方案的设计2.2.2 种群初始化2.2.3 约束处理与适应度函数2.2.4 选择算子2.2.5 交叉算子2.2.6 变异算子2.2.7 其他部分2.3 粒子群算法的求解3. 多种运输工具不带时间窗3.1 立体交通介绍3.2 遗传算法求解3.2.1 染色体编码设计3.2.2 适应度函数3...

前言： 最近在依据概率密度生成随机数方面遇到不少麻烦，找遍了优快云，大多是生成常见分布的随机数，很少有讲生成任意给定概率密度的随机数，即便有部分人写了，但都没给出比较好的模板。有心人天不负，在一番摸索之后，终于是写出来了，顺便还把分位数的两种解法给写出来了，以下是需要解决的问题与相关代码。1. 背景知识1.1 一元p-范分布  在研究线性回归模型(y=βx+ε)(y=\beta x+\va...