庄宝007-CSDN博客

原创有容量限制的车辆路径规划问题(Capacitated Vehicle Routing Problem)

在看matlab的时候发现了这篇文章仔细阅读一下。(英语渣渣，自学用)旨在以尽可能高效的方式为一定数量的客户提供服务。问题在于可能解决方案的组合爆炸，随着客户数量呈超指数级增长。经典解决方案能够提供全局最优解的良好近似。D-Wave的量子退火器是一台旨在解决优化问题的机器。与经典计算机相比，这台机器利用量子效应加快计算速度。在量子退火器上解决CVRP问题的问题在于优化问题的特定表述。

2024-05-21 21:12:21 1099

原创【MATLAB】fminsearchbnd()函数安装包下载

链接：https://pan.baidu.com/s/1apc6_f92rDnV_NwHpZZpDQ?若链接失效，可以自行注册下载，操作也很简单，得到安装包之后，添加路径即可运行。

2024-05-16 15:03:48 583

原创【matlab基础知识代码】（十九）全局最优解&可行解区域

基准测试函数——Rastrigin函数。

2024-05-16 14:49:45 853

原创用阻尼牛顿法求解无约束问题: minf(x),armijo

输入:x0是初始点，fun,gfun,Hess 分别是求目标函数值，梯度，Hesse矩阵的函数。%解方程组Gk*dk=-gk,计算搜索方向。%输出: x,val分别是近似最优点和最优值，k是选代次数。while (m<20) %用Armij0搜索求步长。%功能:用阻尼牛顿法求解无约束问题: minf(x)

2024-05-13 16:20:35 596

原创修正牛顿法求解无约束问题

输入:x0是初始点，fun, gfun,Hess分别是求目标函数值，梯度，Hesse矩阵的函数。end %检验终止准则。dk=-Ak\gk;%解方程组Gk*dk=-gk,计算搜索方向。%计算Hesse矩阵。while (m<20) %用Armij0搜索求步长。%输出:x,val分别是近似最优点和最优值，k是迭代次数。%功能:用修正牛顿法求解无约束问题:min f(x)

2024-05-13 16:17:54 418

原创【matlab基础知识代码】（十八）无约束最优化问题

message: '优化已终止:↵ 当前的 x 满足使用 1.000000e-04 的 OPTIONS.TolX 的终止条件，↵F(X) 满足使用 1.000000e-04 的 OPTIONS.TolFun 的收敛条件↵'绘制三维等高线，获得并叠印中间结果。最优化求解函数的另一种调用方法。建立最优化问题的“结构体”模型。

2024-05-11 17:32:52 505

原创【matlab基础知识代码】（十七）一般非线性方程的数值解方法

精度高得多得多，因为把这个解代回到原始的方程里去，误差是一个非常接近零的数。我们可以修改控制变量，最后能得到双精度意义下的最好的结果。

2024-05-10 17:23:50 297

原创【matlab基础知识代码】（十六）代数方程的图解法&多项式型方程的准解析解方法

1. `subs(F,{x,y,z},{x0,y0,z0})`: 这一部分使用 `subs` 函数，将方程组 `F` 中的符号变量 `x`、`y` 和 `z` 替换为它们的解 `x0`、`y0` 和 `z0`，生成了一个数值的向量。3. `[x0,y0,z0]=vpasolve(F,[x,y,z]), size(x0)`: 这一行利用 `vpasolve` 函数求解方程组 `F`，并将解保存在变量 `x0`、`y0` 和 `z0` 中。方程组的目标是找到满足这三个方程的 `x`、`y` 和 `z` 的值。

2024-05-10 15:58:28 1005

原创【matlab基础知识代码】（十五）线性方程求解（三种情况）

它将Z中的零空间向量与符号变量`[a1;所以，x是矩阵A的零空间向量和方程Ax=B的特解的线性组合。`x0=sym(pinv(A)*B)`：这行代码计算矩阵A的伪逆，将其乘以矩阵B，将结果转换为符号表示，并将其赋值给x0。本质上，它解决了方程Ax=B，找到了x的解。`Z=null(sym(A))`：这行代码计算矩阵A的零空间，将其转换为符号表示，并将其赋值给Z。在这里，`a1`、`a2`、`a3` 和 `a4` 被声明为符号变量。`pinv()` 是计算矩阵的伪逆。`pinv(A)` 会计算矩阵A的伪逆。

2024-05-10 13:51:30 519

原创【matlab基础知识代码】（十四）矩阵的相似变换与三角分解&矩阵的Jordan 变换与奇异值分解

（之前学习牛顿法有涉及到Cholesky的相关证明）Cholesky分解是一种将一个对称正定矩阵分解为一个下三角矩阵和其转置的乘积的方法。简单来讲，就是将一个矩阵分解为一个下三角矩阵和其转置矩阵相乘的结果。Cholesky分解的主要步骤如下：1. 对于一个对称正定矩阵A，可以表示为A = LL^T，其中L是一个下三角矩阵。2. 如果A的第i行第j列元素为a(i,j)，那么L的第i行第j列元素为l(i,j)。

2024-05-09 09:30:00 1422

原创【matlab基础知识代码】（十三）特征值与特征向量

V,D]=eig(H) 先构造相伴矩阵。v为特征值对应的特征向量。

2024-05-09 08:45:00 372

原创【matlab基础知识代码】（十二）逆矩阵与广义逆矩阵

矩阵维数较大，警告: 矩阵接近奇异值，或者缩放错误。结果可能不准确。5 11 10 8;9 7 6 12;4 14 15 1];警告: 矩阵接近奇异值，或者缩放错误。结果可能不准确。B =1.0e+15 *ans =伪逆矩阵。

2024-05-08 17:45:00 612

原创【matlab基础知识代码】（十一）特殊矩阵输入&矩阵性质

删除B的第i行，删除B的第j列。

2024-05-08 16:15:00 416

原创【matlab基础知识代码】（十）级数求和与序列求积

得到收敛区域，级数的和。

2024-05-07 09:30:00 455

原创【matlab基础知识代码】（九）傅里叶级数逼近&泰勒级数

原函数，自变量列表，中心点列表，order，k。尝试不同的展开阶次，6，8，10...16。在这个区间拟合效果比较好。

2024-05-06 09:30:00 384 1

原创【matlab基础知识代码】（八）不定积分的求解

【代码】【matlab基础知识代码】（八）不定积分的求解。

2024-05-06 09:00:00 447 1

原创【matlab基础知识代码】（七）单变量函数求导&偏导数的计算

蓝色的函数原函数，红色的线导数曲线>> syms t;

2024-05-05 08:45:00 772

原创【matlab基础知识代码】（六）极限计算&区间极限与多变量极限

检查对函数 'piecewise' 的调用中是否存在不正确的参数数据类型或缺少参数。分段函数：条件（范围），结果。为啥老师可以运行出来。

2024-05-05 08:30:00 561 2

原创【matlab基础知识】（四）特殊二维图像绘制

窗口分为四个区域，阶梯，火柴杆，直方图，x取对数.y不变。

2024-05-04 10:30:00 243

原创【matlab基础知识】（五）三维曲线绘制

这段代码使用meshgrid函数创建了一个二维的网格，然后根据一定的数学公式计算了一个二维的函数p。最后使用surf函数绘制了这个函数的三维图像。使用surf函数绘制了函数p的三维图像，并通过xlim([-1.5 1.5])设置了x轴的显示范围为-1.5到1.5。用点运算的形式来构造Z矩阵。

2024-05-04 10:00:00 440

原创【matlab基础知识】（三）二维曲线绘制plot

由于0.01cosx波动太小，所以plotyy绘制多纵轴曲线。左蓝为sinx,右黄为0.01cosx。颜色，线形，曲线上的标志。条件和函数值做一个点乘。

2024-05-03 16:35:24 342

原创 CPLEX＋YALMIP安装（matlab调用）

通过百度网盘分享的文件：yalmip+c…链接:https://pan.baidu.com/s/1wwYDu9bHFPmcCg5vD2pSXg?pwd=9d1r提取码:9d1r复制这段内容打开「百度网盘APP 即可获取」具体安装步骤可以参考这篇博客，非常详细在matlab命令窗口中输入“yalmiptest"，看一下运行结果，这样应该就成功了。

2024-05-03 15:41:06 822

原创【牛顿法不需要手动计算梯度黑森矩阵matlab代码】

之前说的不需要手动计算函数梯度和黑森矩阵的牛顿法，但是有的题目求出来的值可能会存在一些错误，谨慎使用。

2024-04-30 17:19:50 296 1

原创【matlab基础知识】（二）

这对函数通常用于评估算法的性能或确定代码中的瓶颈。通过在代码的关键部分之间放置。，可以精确地测量这些部分的执行时间，从而进行优化或者了解其性能。是一对用于测量代码执行时间的函数。(耗时短，尽量采用向量型运算）历时 0.014985 秒。历时 0.002803 秒。在 MATLAB 中，

2024-04-30 17:00:02 796 1

原创【python导入Matplotlib库，附代码测试，适用于新手小白】

其中，`plt.plot`函数用于绘制曲线图，参数`color`、`linestyle`、`linewidth`等可用于设定曲线的颜色、线型、宽度等细节。1. `import matplotlib.pyplot as plt`: 这行代码导入了Matplotlib库，并使用`plt`作为别名，使得在代码中可以简洁地调用Matplotlib的函数。这将产生一个关于x的立方函数的值。5. `plt.plot(x, y)`: 这行代码使用Matplotlib的`plot`函数绘制了x和y之间的曲线图。

2024-04-29 13:04:10 1311

原创【蒙特卡洛方法估计π的值】python，附代码

如果在圆内的点的数量足够大，那么这些点占正方形面积的比例将近似于圆和正方形面积的比例，即 π/4。因此，通过统计落在圆内的点的数量，并将其与总点数的比值乘以4，我们就可以得到一个近似的π值。- 如果点在单位圆内，增加 `total` 计数器，并在图上将该点绘制为红色圆圈 (`'ro'`)。这个函数的参数是`num_samples`，它代表了我们生成的随机点的数量。1. 首先导入了必要的库：`random`、`math`、`sys` 和 `matplotlib.pyplot`。这就是这段代码的要点。

2024-04-29 12:42:58 1006

原创共轭梯度法matlab程序（需要直接复制）

优点就是不需要求Hess矩阵，然后利用精确线搜索方法求步长（感觉优点还是蛮明显的）（内大最优化代码作业）老师说PRP公式最好。

2024-04-27 09:44:31 232

原创 tensorflow（1）———python

让我看看是哪个蠢货,拜托朋友安装了tensorflow(废了好几个小时)，结果现在一点也不会的（没错，正是在下）（这是干什么的，其实我也不太清楚，就安装了）（啊我真的服了我自己）（干脆我撞墙得了）（平静发疯中）好像还是安装了最新版本，朋友还给设置了一个虚拟环境，总之相当麻烦（这就是我比较喜欢matlab的原因，而且我记得之前好像还安装到了2.14版本，总之很混乱）然后再某乎复制了代码，跑了一下，感觉好牛的样子。2. 将生成的样本数据进行可视化，以散点图的形式展示在图上，x轴为房屋大小，y轴为房价。

2024-04-27 08:56:00 1340 2

2301_80118808的博客