- 博客(123)
- 收藏
- 关注
RSS订阅原创 单源最短路径算法
单源最短路径算法是一种用于在有向图或无向图中找到从指定源节点到其他所有节点的最短路径的算法。常用的单源最短路径算法有Dijkstra算法、Floyd-Warshall算法和Bellman-Ford算法。Dijkstra算法是最常用的单源最短路径算法之一，它的基本思想是：从源节点开始，每次选择距离源节点最近的一个未访问过的邻居节点，更新其距离值，直到所有节点都被访问过。
2024-01-11 08:26:15
427
转载 MATLAB画图之自定义图片大小
默认情况下，figure大致在screen的正中间，大致为一个方形（具体参数没有深究），如果我们想要自定义figure在screen的位置和大小，有两种方法可以实现。该指令中的单位为相对值，其中0.3，0.2表示axes的左下角在figure中的距离为figure的宽度×0.3和figure的高度×0.2；axes: figure 中的图像，axes的大小取决于figure，figure大小变化，axes大小也会变化；)指令，该指令设置的是图的坐标线axes的位置及大小，并不是坐标线标注的位置。
2024-01-09 10:24:46
1043
转载 matlab单纯形法
这是目标函数：$Z = {C_B}{B^{ - 1}}b + ({C_N} - {C_B}{B^{ - 1}}N){X_N}$ (b就是增广矩阵多出来的那些)我们需要的检验数就是系数减去${Z_j}$（${\sigma _j} = {c_j} - {z_j}$，${Z_j} = {C_B}{N_j}$）这是运筹学史上最辉煌的阶段。最小比值法则($\frac{{{b_i}}}{{{a_{ik}}}}$，&${a_{ik}}$>0)，使${x_2}$，${s_1}$，${s_2}$对应的基矩阵变成。
2024-01-08 10:58:11
15
转载 矩阵QR分解的MATLAB与C++实现
矩阵的QR分解目的是将一个列满秩矩阵\(A\)分解成\(A=QR\)的形式，我们这里暂时讨论\(A\)为方阵的情况。首先假设初始方阵为\(A\)，\(\vec{x_i}\)、\(\vec{y_i}\)、\(\vec{z_i}\)都为列向量。可以看出其实矩阵的QR分解的步骤并不多，就是不断地循环进行\(A\)的正交化、标准化、求\(Q\)、求\(R\)这几步。现在主要的问题是如何由矩阵\(A\)计算得到矩阵\(Q\)和\(R\)呢？注：\(A_{:k}\)表示\(A\)的第\(k\)列向量。
2024-01-04 16:47:51
57
原创 matlab 绘图坐标轴设置
x 坐标轴设置成 - pi/2 到 pi/2，间隔为 pi/4。上面的图坐标轴默认格式，不太符合数学习惯。y 坐标轴设置成 - 1 到 1，间隔为 0.5。简单的曲线：y=sin (x);
2024-01-04 10:17:25
614
转载 Matlab-simulink仿真2ASK调制
采用的是Zero-Order Hold模块(零阶保持器)，作用是将输入信号每过一个采样时间更新一次，并保持到下一次采样，多用于将连续信号离散化(这里不太好理解，如果不加这个模块，AWGN模块就会报错)使用的是Relay模块，由于基带信号最大值为1，加上噪声后会在1左右浮动，这里设置判决值为0.3，即采样值幅度大于0.3时，认为有信号，输出1，采样值幅度小于0.3时，认为无信号，输出0，解调结束。使用的是AWGN(高斯白噪声)模块，模式设置为SNR(信噪比)，将信噪比设置为15，信噪比越小，噪声干扰越大。
2024-01-03 12:11:08
325
原创 matlab将数据读取和写入txt文档
另外，在这些字符串后添加一个“t”，如‘rt’或‘wt+’，则将该文件以文本方式打开；如果添加的是“b”，则以二进制格式打开，这也是fopen函数默认的打开方式。（注：当用‘a’时，如果文本中已经存在数据，不会清空数据，而是在数据之后写入，而‘w’会清空原本的数据，重新写入）说明：fid用于存储文件句柄值，如果fid>0，这说明文件打开成功。三，如果要手动保存到指定位置，指定文件名。一，如果要保存单行、单列数据。matlab中打开文件。二，如果要保存一个矩阵。matlab读取文档。
2024-01-02 10:41:15
790
转载 Simulink混合仿真系统实例练习——汽车行驶控制系统
位置变换器是汽车行驶控制系统的输入，其作用是将速度操纵机构的位置转换为相应的速度，速度操纵机构的位置和设定速度间的关系为：Vg=50x+45 x=(0~1)•第二部分，离散PID控制器。其数学模型为：积分环节：x(n)=x(n-1)+u(n)微分环节：d(n)=u(n)-u(n-1)控制器输出：y(n)=P*u(n)+I*x(n)+D*d(n)其中P,I,D分别是PID控制器的比例、积分和微分控制参数。例：汽车行驶控制系统是应用很广的控制系统之一，控制的目的是对汽车速度进行合理的控制。
2024-01-01 17:33:51
169
转载 基于Simulink的SSB信号调制与解调仿真
在做滤波器设计时，最先使用的巴特沃斯滤波器，由于过渡带比较宽，在进行边带滤波时，总会残留一些边带，所以最后选择了边带窄的椭圆滤波器，实现了比较好的边带滤波效果。按照书本上的原理，加性噪声为高斯白噪声，然后经过BPF后进行相干解调，但由于版本问题，没有高斯白噪声模块，所以我这里就直接用了带限白噪声，相当于高斯白噪声通过BPF后的带限白噪声，然后进行相干解调，解调出来的信号，与基带信号相比，会有不少的畸变，当然因为实际情况更为复杂，所以我的仿真模块一些参数的设计并不能满足实际，所以这是我需要提升的。
2024-01-01 16:13:34
78
转载 线性卷积 matlab 实现
当选取N = max（N1，N2）做循环卷积时，前（M-1）个样本在误差中(即有别于线性卷积)，这里M = min（N1，N2）。如果将X1 与X2 通过补零而成为N =（N1+N2-1）点序列，那么循环卷积就与线性卷积一样了。实现中的N过大时会占用较多存储空间 当N选取的值比要求值小时就会引入误差。%x3 = conv(x1,x2) %MATLAB 自带的函数。circonvt(x1,x2,9) %前面一节刚写的函数。差别在于循环移位和线性移位。
2024-01-01 12:21:31
47
转载 一种采用匹配滤波器进行静脉图像分割方法的matlab实现
如果将血管图像看做一个信号，虽然这个信号的具体情况不得而知，但是血管的特点是已知的，基于这点先验知识可以针对血管设计匹配滤波器，那么按照匹配滤波器的原理，当血管部分输入时会出现较大值，而当背景区域输入时将出现较小值，从而将静脉与血管分开。出于一些医疗上的目的，经常需要对一些医疗成像的血管图像进行分割，提取出血管部分，来进行进一步的研究，所以血管分割作为图像分割中的一类，有很多关专门解决这个问题的算法被提出。至于全局阈值t的选择，最粗犷最常用的方法就是试几次找最好的，当然也有其他找这个阈值的方法。
2024-01-01 11:17:00
37
转载 Matlab图像处理工具箱用户指南；裁剪图像及空间变换部分翻译
变换后的图像看起来和原来的图像完全一样。例如，您可以指定您自己的可分离的滤波或插值核心，用MATLAB的interp2和interp3函数构造一个自定义的重采样，甚至可以实施先进的抗混叠技术。imtransform支持几个可选输入参数，您可以通过它们来控制变换的几个方面，例如输出图像的大小和使用的填充值。像素点(1,1)在输入图像中是黑色的，输出图像中的像素点(41,41)也是黑色的。为了看到相对于原图像的变换效果，您可以使用几个可选的输入参数以指定输出图像的大小和输出图像中要包含多少输出空间。
2023-12-30 13:52:51
79
转载 遇见蒙特卡洛
可以看到，随着随机点数的增加，算得的pi精度越来越高。（注意这里，我将脚本文件命名为pi.m，系统会有冲突警告，所以最好改下文件名。昨天，呃，不，是今天凌晨，我问一位大牛有没有什么既简单又强势的算法，他说了蒙特卡洛。今天查了些资料，见识到它的强大与及简洁。这里参考一篇论文，讲解关于测度的MATLAB实现。这里以计算y=x^2在[0,1]上的定积分为例。最后，看到知乎上有句话总结了蒙特卡洛方法——现在用MATLAB实现蒙特卡洛方法的几个应用。可见结果运行还是很稳定的。多么值得思考的一句话啊。
2023-12-30 10:38:15
21
原创 Matlab曲线的颜色、线型等参数设置方法
在Matlab绘图时，如果在同一视窗内绘制的曲线特别多的时候，通常需要用不同的颜色和线型将其区分开，否则在黑白色打印输出图像时，用户无法辨别不同属性的曲线。下面是的曲线颜色的配置模板，大家可以直接修改使用，直接引用PlotStyle这个数据结构就可以了。
2023-12-29 11:54:51
863
原创 hough变换是如何检测出直线和圆的？
如(2)中说明，过点（x1,y1）的所有圆可以表示为（a1(i),b1(i),r1(i)），过点（x2,y2）的所有圆可以表示为（a2(i),b2(i),r2(i)），过点（x3,y3）的所有圆可以表示为（a3(i),b3(i),r3(i)），如果这三个点在同一个圆上，那么存在一个值（a0,b0,r0），使得 a0 = a1(k)=a2(k)=a3(k) 且b0 = b1(k)=b2(k)=b3(k) 且r0 = r1(k)=r2(k)=r3(k)，即这三个点同时在圆（a0,b0,r0)上。
2023-12-28 12:32:52
791
转载 使用matlab进行傅里叶分析和滤波
一般而言，N点fft的结果y，在$n=N/2+1$处对应的频率为最高采样率的一半，y的后一半与前一半对称。下例是将振幅为1的5Hz正弦波和振幅为0.5的10Hz正弦波相加之后，滤除8Hz以上的信号。下例 是将振幅为1的5Hz正弦波和振幅为0.5的10Hz正弦波相加之后进行傅里叶分析。下例 是将振幅为1的5Hz正弦波和振幅为0.5的10Hz正弦波相加之后进行傅里叶分析。以巴特沃斯低通滤波器为例，说明调用方法。以巴特沃斯低通滤波器为例，说明调用方法。对应的逆变换有两种，分别为。
2023-12-27 15:09:40
235
转载 计算模糊控制离线查询表
设X、Y、Z∈{-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6}，T(x)={NB,NM,NS,NZ,PZ,PS,PM,PB}，T(y)=T(z)={NB,NM,NS,ZE,PS,PM,PB}比如，假设x输入为-4，则输入的模糊集合为：$A'=\frac{0}{-6}+\frac{0}{-5}+\frac{1}{-4}+\frac{0}{-3}+...+\frac{0}{5}+\frac{0}{6}$当输入的维数较高，即有很多个模糊子句用and相连时，模糊推理的计算便比较复杂。
2023-12-27 10:09:02
33
原创 matlab-调用摄像头人脸识别
VJ算法在Matlab里面实现的时候，已经训练好了正脸、侧脸、上半身、眼睛、嘴、鼻子，这些都是可以直接检测，不需训练，直接调用CascadeObjectDetector函数即可。Installer在MATLAB里面有下划线，然后你点开它，MATLAB会提供大概13个软件包，这时候选择。至于调用摄像头进行人脸识别，肯定是 调用摄像头的过程中对每一帧图像分别进行识别，然后再在图像中框出来。VJ算法的目的是检测人脸，但是其思想同样可以用于检测其他物体，只需进行训练即可。所以检测的时候就要压缩你图像的像素了。
2023-12-25 08:55:09
388
原创 MATLAB图像处理相关函数
语法：P=impixel(I) \ P=impixel(X,map) \ P=impixel(RGB) \ P=impixel(I,c,r) \ P=impixel(X,map,c,r) \ P=impixel(RGB,c,r) \ [c,r,P]=impixel(...) \ P=impixel(x,y,I,xi,yi) \ P=impixel(x,y,RGB,xi,yi) \ P=impixel(x,y,X,map,xi,yi) \。
2023-12-25 08:33:56
919
原创 matlab saveas 和imwrite的区别
1、背景色：saveas保存的图像 背景色自动设置为白色，imwrite保存图像为所见即所得。2、图像大小： saveas无视你设置的图像大小，按默认保存，imwrite保存所见即所得。% 在当前文件夹下新建image文件夹，如果已存在会warning，不影响运行。set(gcf,'color',[1,1,1]) % 背景色设置为白色。saveas(handle，['目录','文件名'])% 将图像设置为距屏幕左下角 [80，100]像素。% 图像大小设置为400*600像素。需要与getframe连用。
2023-12-22 20:21:10
373
转载 matlab 全局变量的使用举例
全局内存是指内存的一种特殊类型，它能够被所有的工作区访问。如果一个变量在函数中被声明全局变量，那么它将占用的是全局内存，而不是本地工作区。昨天在写项目时，想要把获取到的临时变量放入一个全局变量，为以后的使用做准备，结果总是出错，今天做了一个小程序，放在这里备用。如果相同的变量在另一个函数中被声明为全局变量，那么这个变量所占有内存区域就是第一个函数中的相同变量。在一个程序，你能利用全局内存，在函数之间对大规模数据进行交换。最好把全局变量声明在函数的开头，这样可以区别于本地变量。
2023-12-22 08:55:16
215
原创 Matlab自带的曲线拟合程序
因为这种拟合方法用的比较少，所以特此总结出来，希望后来对自己的道路段的拟合会有用处！这个函数的功能是能自动搜索参数的取值，从而使得方程的误差最小。
2023-12-21 17:37:37
361
原创 matlab函数_连通区域
如果properties等于字符串'basic'，则属性：'Area'，'Centroid'和'BoundingBox'将被计算。B是一个P×1的cell数组，P为对象个数，每个cell 是Q×2的矩阵，对应于对象轮廓像素的坐标。L中不同的正整数元素对应不同的区域，例如：L中等于整数1的元素对应区域1；返回值STATS是一个 长度为max(L(:))的结构数组，结构数组的相应域定义了每一个区域相应属性下的度量。'Orientation' 与区域具有相同标准二阶中心矩的椭圆的长轴与x轴的交角（度）
2023-12-21 13:52:28
711
原创 MATLAB简易验证码识别程序介绍
对于部分特殊的验证码，需要对分割后的图片进行标准化处理，也就是说尽量把每个相同的字符都变成一样的格式，减少随机的程度。最简单的比如旋转还原，复杂点的比如扭曲还原等等。比如本文中分割后的数字1和8宽度不一致，把他们的宽度填充一致，就是标准化的一种。这一步可以用很多种方法，最简单的就是模板对比，对每个出现过的字符进行处理后把点阵变成字符串，标明是什么字符后，通过字符串对比来判断相似度。所谓二值化就是把不需要的信息通通去除，比如背景，干扰线，干扰像素等等，只剩下需要识别的文字，让图片变成2进制点阵。
2023-12-19 19:46:39
1020
原创 FFT的物理意义
假设采样频率为Fs，采样点数为N，做FFT之后，某一点n（n从1开始）表示的频率为：Fn=(n-1)*Fs/N；该点的模值除以N/2就是对应该频率下的信号的幅度（对于直流信号是除以N）；该点的相位即是对应该频率下的信号的相位。相位的计算可用函数atan2(b,a)计算。atan2(b,a)是求坐标为(a,b)点的角度值，范围从-pi到pi。要精确到xHz，则需要采样长度为1/x秒的信号，并做FFT。要提高频率分辨率，就需要增加采样点数，这在一些实际的应用中是不现实的，需要在较短的时间内完成分析。
2023-12-19 16:50:19
872
原创 Matlab文件操作
判断文件位置指针是否达到文件结束位置，若文件位置指针已在文件末尾，则返回 1，否则返回 0：eofstat=ftell(fid)7、 其中 fid 为所要关闭的文件的句柄，status 为关闭文件的返回代码，若关闭成功则为 0，否则为 -1。4、 有两个标准代码文件，不需打开就可以直接使用，分别为：fid=1 标准输出文件，fid=2 标准错误文件。6、在Windows下，打开文件的默认是二进制格式，如果要以文本方式打开，则必须在打方式中加上字符 't'。写文件，若文件已存在则原内容将被覆盖；
2023-12-19 10:55:03
787
原创 matlab特征值分解和奇异值分解
函数 eig格式 d = eig(A) %求矩阵A的特征值d，以向量形式存放d。d = eig(A,B) %A、B为方阵，求广义特征值d，以向量形式存放d。[V,D] = eig(A) %计算A的特征值对角阵D和特征向量V，使AV=VD成立。[V,D] = eig(A,'nobalance') %当矩阵A中有与截断误差数量级相差不远的值时，该指令可能更精确。'nobalance'起误差调节作用。
2023-12-19 10:48:53
611
转载 Matlab指针
由上面所述内容可看出，handle类型对象在赋值给其他变量时，不会对原始对象进行复制。此外，在函数中对handle类型对象做修改后，也不需在将其作为返回参数输出。也就是说，在使用handle类型对象时，我们可以认为所操作的就是一个指针，其指向的内容就是构造函数初始化的实例。第一印象貌似是Matlab中不存在指针，所有变量与函数的赋值都是按值传递的，不会对参数进行修改。这句代码实现了两个功能，将类实例化了一个对象（object），并将该对象的句柄赋给变量。对实例化的对象进行的操作直接作用在。
2023-12-19 09:38:22
135
原创 matlab中norm与svd函数用法
格式：n=norm(A,p)功能：norm函数可计算几种不同类型的矩阵范数,根据p的不同可得到不同的范数。
2023-12-18 14:48:08
317
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人