MATLAB
文章平均质量分 84
知行zhe
这个作者很懒,什么都没留下…
展开
-
MATLAB---初识实例1
MATLAB中一些有趣函数logo(MATLAB的logo),earthmap(地球仪),travel(旅行商问题),penny,fifteen(小游戏),teapotdemo(茶壶)clear 清空工作区 clc清空命令行窗口预定义变量名:pi ans epseps==eps(1) 0.5+eps(x)界限虚数单位:i 无穷大...原创 2018-09-24 18:07:07 · 507 阅读 · 0 评论 -
MATLAB--非线性、线性方程组的求解
%%牛顿法和牛顿下山法% h=@(x) x.^3+x.^2-1;% x = newton_eq(h,1,1000);function result=newton_eq(h,x,n)%x为初始值%n为迭代次数f = matlabFunction(h);%转化为句柄形式f1=matlabFunction(diff(h));%weishixianX(1)=x;i=2;while 1...原创 2018-10-11 23:08:28 · 4631 阅读 · 0 评论 -
MATLAB矩阵
1.矩阵建立(1)常规表示:eg. A=[1,2,3:4,5,6:7,8,9] //此处表示三行三列,同一行用逗号或者空格分隔,不同行元素用分号分隔(2)矩阵的拼接 eg. A=[1,2,3:4,5,6:7,8,9]; B=[-1,-2,-3;-4,-5,-6;-7,-8,-9];C=[A,B;B,A]复数矩阵 C=A+i*B;2.冒号表达式e1:e2:e3 //e...原创 2018-09-21 23:58:38 · 2051 阅读 · 0 评论 -
MATLAB基础
1.数值数据类型分类 (1)整型(有无符号均八位) eg x=int8(129) x=127(范围为-128---127) x=uint8(129) x=129(范围为0---256) (2)浮点型(single函数:转为单精度;double函数:转为双精度) (3)复型(r...原创 2018-09-21 21:57:31 · 417 阅读 · 0 评论 -
MATLAB--条件插值、二维插值、拟合
%% 样条边界条件插值x=linspace(0,2*pi,15);y=sin(x);plot(x,y,'o')hold on% completep=csape(x,y,'complete',[2,2]);x1=linspace(0,2*pi,150);y1=ppval(p,x1);plot(x1,y1,'b','linewidth',1.5)% secondp=csape(x,y...原创 2018-10-15 18:00:41 · 13252 阅读 · 1 评论 -
MATLAB--插值
%多项式插值x=[1 2 3 4 5 6];y=[2 5 6 9 14 26];a=duoxiangshi(x,y)function D=duoxiangshi(X,Y)A=vander(X)Y=Y'B=[A,Y]C=rref(B)D=C(:,end)fprintf('\n 这是一个%d阶的多项式,从高阶到低阶排列!\n',length(Y)-1);plot(X,Y,...原创 2018-10-15 17:57:24 · 2149 阅读 · 0 评论 -
MATLAB--符号与数值
%符号表达式的整理%% 化简 syms a by=(cot(a/2)-tan(a/2)*(1+tan(a)*tan(a/2)));y=simplify(y);%y=simple(y) 已经删除%% 因式分解 syms m n xy=-24*m^2*x-16*n^2*x;y=factor(y);y1=m^3-n^3;y1=factor(y1);y2=factor(104);%% ...原创 2018-10-14 20:48:48 · 3292 阅读 · 0 评论 -
MATLAB--曲面图、三维图
%% mesh 网格图[X,Y]=meshgrid(-5:0.5:5);R=sqrt(X.^2+Y.^2)+eps;%加上eps防顶点缺口Z=sin(R)./R;mesh(Z)mesh(X,Y,Z)h=meshc(Z) 下图一meshz(X,Y,Z) 下图二%% 固定化纵横比 axis vis3d%% 透明hidden on%看不...原创 2018-10-03 22:24:54 · 23177 阅读 · 1 评论 -
MATLAB--绘图4
x=[0 1 1 0 0];y=[0 0 1 1 0];t=[1 1 1 1 1];plot(x,y)h=patch(x,y,'b');h.FaceAlpha=0.5;hold on% axis([-4 4 -4 4])% axis equal% M=[1 0 2;% 0 1 3;% 0 0 1];% M1=power(M,1/50);axis([-6 6 -...原创 2018-11-03 19:56:47 · 249 阅读 · 0 评论 -
MATLAB--图形旋转、压缩、镜像、剪切
%% 旋转x=[0 1 1 0 0];y=[0 0 1 1 0];plot(x,y)h=patch(x,y,'b');h.FaceAlpha=0.5;hold onaxis([-2 2 -2 2])axis equalfor i=1:40 %旋转40次 M=[cos(pi/40*i) -sin(pi/40*i); sin(pi/40*i) cos(pi/...原创 2018-10-02 16:15:58 · 1039 阅读 · 0 评论 -
MATLAB--绘图3
%% pie(x) 饼图% x=[0.1 0.2 0.4];%小于1% pie(x) 下图一%% 大于一x=[3 6 7]di=[0 0 1]%di=zeros(1,length(x));di(3)=1;label={'第一产业','第二产业','第三产业'}p=pie(x,di)%p=Patch Text Patch Text Patc...原创 2018-10-02 12:43:10 · 385 阅读 · 0 评论 -
MATLAB--绘图2
%% logspaceA=linspace(1,6)%从1到6一共100个数A1=logspace(1,6)%从10^1到10^6一共50个数,默认生成100个数A2=logspace(1,6,100)%从10^1到10^6一共100个数A3=logspace(1,pi)%生成的是10到pi的一共50个数,若写成3.1415926535……生成的不是pi%% loglogY=[1 5 1...原创 2018-10-02 09:22:23 · 445 阅读 · 0 评论 -
MATLAB--绘图1
%多图绘制x=0:0.1:2*pi;y1=sin(x);y2=cos(x);y3=sin(x).*cos(x);y4=1./x;% %subplot(n1,n2,n3)n1表示行数n2,表示列数,n3表示第几幅图% subplot(2,2,1)% plot(x,y1)% title('sin')% subplot(2,2,2)% plot(x,y2)% title('cos...原创 2018-10-01 17:06:55 · 302 阅读 · 0 评论 -
MATLAB--实例3(矩阵)
%矩阵的加减,数乘,转置 %矩阵的加法,矩阵与矩阵相加,维度必须相同% A=[1 2 3;4 5 6]% B=[3 5 2;7 4 1]% C=A+B; %矩阵与数相加,矩阵的每一个数都相加% D=A+1; %矩阵的减法同加法 %矩阵的数乘就是每一个元素都乘以这个数% A=[4 5;2 4]% % B=2*A;% %矩阵的转置,行与列交换% B=A';%对...原创 2018-09-28 11:04:28 · 4426 阅读 · 0 评论 -
MATLAB--实例2(矩阵)
%特殊矩阵-空矩阵、全1矩阵、单位矩阵%空矩阵% A=[2 3 4;3 1 5;7 8 4;]% A(2,:)=[]%删除第二行数%全一矩阵 ones()% A=ones(4)% B=ones(4,3)% C=zeros(3)% C(:,:)=1%单位矩阵 eye()% A=eye(3);%对角线上都是1% B=eye(4,5);% D=zeros(3)% for m=1...原创 2018-09-25 12:47:09 · 2477 阅读 · 0 评论 -
MATLAB--plot详解
% x=0:0.1:2*pi;% y=sin(x);% plot(x,y)默认绘图% plot(x,y,'--')%虚线% plot(x,y,':')% plot(x,y,'-.')%点画线% plot(x,y,'y')% plot(x,y,'b:')% plot(x,y,'-.bo')%圆形标记% plot(x,y,'-.bh')%六角星标记%线宽LineWidth% plo...原创 2018-10-01 11:18:31 · 7948 阅读 · 2 评论