
【Matlab代码】数学实验
数学实验
裴来凡
一个爱大笑的女孩。
展开
-
【Matlab代码】分形树木
程序设计的思路:1,我们先设定主干;2,我们将主干设定的比例缩小;3,将第二步得到的缩小的主干旋转一定的角度,然后接到主干上面预定的分叉的位置。原创 2020-10-16 13:46:38 · 1678 阅读 · 0 评论 -
【Matlab代码】Sierpinsk地毯
1,将一个大正方形作为最基本的图形;2,将第一步得到的数组缩小到三分之一;3,将第二部得到的数组复制八份,平移到相应的位置。原创 2020-10-16 11:22:40 · 2134 阅读 · 5 评论 -
【Matlab代码】Koch曲线
步骤:1,将基本的图形即线段用数组表示出来。由于是线段,我们只需要列出端点坐标即可;2,将第一步得到的数组做如下处理:首先将线段缩小到原来的三分之一。然后将得到的线段复制四份:将第一段放在最左边,将第二段逆时针旋转60°,左端点与第一段的右端点相接;将第三段顺时针旋转60°,左端点与第二段的右端点相接;最后将第四段左端点与第三段右端点相接。这样我们得到了最初屏幕上的第五个图。...原创 2020-10-16 11:08:11 · 1861 阅读 · 0 评论 -
【Matlab代码】斑点猫头鹰的生存问题的演化过程
A=[0 0 0.33;0.18 0 0;0 0.71 0.94];[pc,lambda]=eig(A);%求A的特征值和对应的特征向量[Y,I]=sort(diag(abs(lambda)),'descend')%对特征值的绝对值降序排列temp=diag(lambda);lambda=temp(I)%输出按特征值的绝对值降序排列的特征值Lambda_norm=[norm(lambda(1);norm(lambda(2));norm(lambda(3))]%对三个特征值的模pc=pc(:,I)原创 2020-10-16 10:24:15 · 690 阅读 · 0 评论 -
【Matlab代码】捕食者与被捕食者问题的演化过程
求A的特征值与特征向量。A=[0.5 0.4;-0.104 1.1];[pc,lambda]=eig(A);%求A的特征值和对应的特征向量[Y,I]=sort(diag(abs(lambda)),'descend');%对特征值的绝对值降序排列temp=diag(lambda);lambda=temp(I)%输出按特征值的绝对值降序排列pc=pc(:,I)%与特征值对应的特征向量...原创 2020-09-25 21:19:51 · 1714 阅读 · 0 评论 -
【Matlab代码】特征值与特征向量基本概念
函数:[V,D]=eig(A)或[V,D]=eig(A,‘nobalance’)功能:求矩阵A的特征值和特征向量说明:(1)生成特征矩阵D和特征向量构成的矩阵V,使得AV=VD。矩阵D由A的特征值在主对角线构成的对角矩阵。V是由A的特征向量按列构成的矩阵。(2)[V,D]=eig(A)中,先对A作相似变换再求A的特征值和特征向量;而[V,D]=eig(A,‘nobalance’)中,直接求矩阵A的特征值和特征向量。...原创 2020-09-25 20:40:59 · 1312 阅读 · 0 评论 -
【Matlab代码】有持续外力下单自由度阻尼系统的数值解
function dx=odefun4(t,x)global c;%全局变量dx=zeros(2,1);%存储x1,x2的导数dx(1)=x(2);%第一个方程dx(2)=20*sin(5*t)-20*c*x(2)-100*x(1);%第二个方程注意:修改为:dx(2)=20sin(5t)-20cx(2)-100*****x(1);%第二个方程原创 2020-09-25 20:13:11 · 645 阅读 · 0 评论 -
【Matlab代码】外力单自由度阻尼系统的数值解
function dx=odefun2(t,x)dx=zeros(2,1);%存储x1,x2的导数dx(1)=x(2);%第一个方程dx(2)=-20*x(2)-100*x(1);%第二个方程[t,x]=ode45('odefun2',[0,4],[1;0])原创 2020-09-25 19:39:14 · 2574 阅读 · 0 评论 -
【Matlab代码】无外力下单自由度阻尼系统的解析解
x=dsolve('D2x+2*c*w*Dx+w^2*x=0')原创 2020-09-25 19:02:10 · 846 阅读 · 0 评论 -
【Matlab代码】微分方程的解析解和数值解
Matlab求微分方程的解析解用法:dsolve(‘equation’,‘condition;,‘v’)说明:(1)equation是方程式,condition是条件,v是自变量(确省为t)(2)若不带条件,则解中带积分常数。(3)如果没有显示解,则系统尝试给出隐氏解。(4)如果无隐氏解,则返回空符号。格式:(1)y’表示为Dy,y’‘表示为D2y,依次类推。(2)有多个方程或多个条件时,写多个对应的参数即可。求微分方程y’’-5y’+6y=e^axy=dsolve('D2y-5*Dy+6原创 2020-09-25 18:13:15 · 24566 阅读 · 1 评论 -
【Matlab代码】生成物浓度散点图
化学反应生成物浓度原创 2020-09-25 14:38:57 · 3294 阅读 · 0 评论 -
【Matlab代码】斐波拉契数列
实验:显示数列的函数类型function fib1(n)%显示数列前n项fn=[1,1];%将前两项放到数组fn中for i=3:n fn=[fn,fn(i-2)+fn(i-1)];%将第i项添加到数组fn中endplot(fn)%前n项显示成曲线原创 2020-09-25 11:10:28 · 4463 阅读 · 0 评论 -
【Matllab代码】数据拟合简介
做数据拟合的方法:在matlab软件中,做数据拟合的函数是polyfit。函数:polyfit(x,y,n)功能:用n阶多项式拟合数据(x,y),使得在数据点处误差的平方和最小。说明:参数x和y都是数组,分别存放数据的横坐标和竖坐标;参数n是多项式的阶,是人为指定的。函数polyfit的返回值,只是多项式的系数。如果需要计算多项式的值,调用函数polyva即可。函数:polyval(p,x)功能:根据多项式的系数,计算多项式的值。说明:参数p,是函数polyfit的返回值,参数x,是指定的横原创 2020-09-25 10:57:25 · 378 阅读 · 0 评论 -
【Matlab代码】程序设计
(1)for循环:for i=范围表达式可执行语句end(2)while循环:while 逻辑表达式可执行语句end(3)if-else-end分支:if 逻辑表达式可执行语句else if 逻辑表达式可执行语句else可执行语句end说明:上述结构中,可以没有else if段和else段。(4)switch-case-end分支:switch 表达式case 常量表达式1可执行语句块1case 常量表达式2可执行语句块2case 常量表达式n可执行语句块n原创 2020-09-24 21:47:45 · 179 阅读 · 0 评论 -
【Matlab代码】M文件
将要执行的命令,写在1个文件中,就是M文件。如果需要执行的命令较多,或者需要重复执行命令,就应该采用M文件方式。MATLAB自带M文件编辑器Notebook,M文件的扩展名是.m。M文件应保存在current directory(当前工作目录)下。...原创 2020-09-21 10:52:58 · 1727 阅读 · 0 评论 -
【Matlab代码】二维特殊图形
常见的统计图,都有对应的函数。原创 2020-09-21 10:39:29 · 364 阅读 · 0 评论 -
【Matlab代码】hold指令
将多次绘制的曲线,显示在同一图形窗口中。函数:hold on,hold off功能:hold on,开启图形保持功能;hold off,关闭图形保持功能。说明:MATLAB默认的配置是hold off。原创 2020-09-21 10:37:14 · 961 阅读 · 0 评论 -
【Matlab代码】用ezplot函数绘制
函数:ezplot(F,[xmin,xmax])功能:绘制符号函数F在区间[xmin,xmax]内的图像。参数说明:参数F是函数表达式,只含有一个变量,参数区间可缺省,默认为[-2pi,2pi]。注释:可以用ezplot函数,绘制参数方程与隐函数的图形,不能设置图形的格式。ezplot('cos(4*x)*sin(10*x)*exp(-abs(x))',[-3,3,-1,1])说明:参数[-3,3,-1,1]是指定x坐标和y坐标的范围,这是为了得到更好的效果。...原创 2020-09-21 10:33:58 · 5752 阅读 · 0 评论 -
【Matlab代码】plot函数
绘制二维图形,最基本的函数是plot。函数:plot(X,Y,‘s’)功能:连接参数所给点列,形成一条或多条折线。参数说明:参数(X,Y)是点列的坐标,二者是同长度的数组或矩阵。参数s是图形格式字符串,可以设定曲线的线型和颜色。注释:横坐标紧密时,效果就是曲线图。格式字符串可以缺省,由系统自动设置格式。可以有多组参数。增加一组参数,就多画一条曲线。X=-4:0.1:4;Y=(exp(X)+exp(-X))*0.5;plot(X,Y)X=-4:0.1:4;Y1=(exp(X)+原创 2020-09-21 10:20:04 · 1417 阅读 · 1 评论 -
【Matlab代码】数组及其运算
数组就是,只有一行或一列的矩阵,数组具有矩阵的所有特性,同时还有创建数组的特殊指令,还有异于矩阵运算的数组运算。特殊数组的创建:linespace(a,b,n)给出区间[a,b]内的n个等分点构成的数组。数组的运算:...原创 2020-09-21 09:40:25 · 598 阅读 · 0 评论 -
【Matlab代码】矩阵的运算
A’:A的转置;A*B:矩阵乘法;A-B:矩阵减法;A+B:矩阵加法;A/B:矩阵的右除;A\B:矩阵的左除;det(A):A的行列式;rank(A):A的秩。注意:左除与右除的区别,左除,是我们常规理解的除法。...原创 2020-09-21 09:34:28 · 1589 阅读 · 0 评论 -
【Matlab代码】对矩阵元素的操作
A(i,j):表示矩阵A的第i行第j列元素;A(:,j):表示矩阵A的第j列;A(i,:):表示矩阵A的第i行;A(i):表示矩阵A(:)的第i个元素。A=[1,2,3;4,5,6;7,8,9]A(2,3)原创 2020-09-21 09:24:15 · 4223 阅读 · 0 评论 -
【Matlab代码】矩阵的操作与运算
操作符”:“,用于增量操作。j:k:步长为1,从j到k,按等差数列构造数组;j:i:k:步长为i,按等差数列构造数组。1:51:2:78:-2:0原创 2020-09-21 09:12:53 · 1049 阅读 · 0 评论 -
【Matlab代码】矩阵的创建
1.直接输入创建矩阵:(1).“[]”作为首尾;(2).“,”或空格隔开;(3).";"或Enter键分隔。A=[1,2,3;4,15,60;7,8,9]2.用矩阵函数生成矩阵:MATLAB提供了大量的函数,来创建一些特殊的矩阵。eye(3)zeros(3,4)...原创 2020-09-21 09:06:00 · 2136 阅读 · 0 评论 -
【Matlab代码】处理函数句柄的函数。
fhda=@expfhdb=@myprocessfunctions(fhdb) isa(fhda,'function_handle')isequal(fhda,fhdb)原创 2020-02-22 20:23:51 · 271 阅读 · 0 评论 -
【代码】图像处理基础
1、字符串 A='MATLAB 7.0'whos2、命令和语句输入2*sin(0.3*pi)/(1+sqrt(5))3、矩阵的输入 A=[1,2,3 4,5,6 7,8,9]4、语句与变量(1)1996/18(2)s=1-1/2+1/3-1/4+1/5-1/6+... 1/7+1/8;s5、矩阵的下标x=4:6...原创 2019-11-29 21:32:42 · 288 阅读 · 0 评论 -
【代码】图像处理基础
1、图像捕获工具箱out1 = imaqhwinfo原创 2019-11-30 19:06:21 · 307 阅读 · 0 评论 -
【Matlab代码】形状识别及定位
A=imread('coins.png');imshow(A)[centers,radius]=imfindcircles(A,[15 30]);centersStrong5=centers(1:5,:);radiiStrong5=radius(1:5);viscircles(centersStrong5,radiiStrong5,'EdgeColor','B')rgb=imrea...原创 2020-02-10 20:29:07 · 4055 阅读 · 4 评论 -
【Matlab】图像处理——图片的缩小
https://mp.weixin.qq.com/s/evv0RJdumfshwQkXJe0Q6w原创 2020-01-28 21:14:07 · 865 阅读 · 0 评论 -
【代码】图像处理——边缘检测算子的算法
也可以是邻接矩阵。原创 2020-01-15 16:38:07 · 752 阅读 · 0 评论