matlab生成对角矩阵 如何用较少的语句生成这个复杂的矩阵

&MATLAB小技巧之十六:利用MATLAB绘制矩阵色块图&这福利我要幸福死了! - MATLAB的日志,人人网,MATLAB的公共主页
好久不发状态了,大家还在吗?还在学习MATLAB么?
&MATLAB小技巧之十六:利用MATLAB绘制矩阵色块图&这福利我要幸福死了!
更多技巧资料请关注微信公共账号:数据挖掘的那些事
R语言中有一个根据实值矩阵绘制色块图的程序(用于绘制相关系数矩阵图),可以用丰富的颜色和形状形象的展示矩阵元素值的大小。遗憾的是MATLAB中没有这样的函数,因此我就用MATLAB编写了一个matrixplot函数,下面给出示例和源码,与大家一起分享!
**********************************************
【例1】绘制网格线,网格中显示矩阵元素。
显示黑色文字:
&& x = [1,-0.2,0.3,0.8,-0.5&& -0.2,1,0.6,-0.7,0.2&&& 0.3,0.6,1,0.5,-0.3&&& 0.8,-0.7,0.5,1,0.7&& -0.5,0.2,-0.3,0.7,1];&& XVarNames = {'xiezhh','heping','keda','tust','tianjin'};&& matrixplot(x,'FillStyle','nofill','XVarNames',XVarNames,'YVarNames',XVarNames);
根据矩阵元素值自动设置文字颜色:
&& matrixplot(x,'FillStyle','nofill','XVarNames',XVarNames,'YVarNames',XVarNames,'TextColor','Auto','ColorBar','on');
【例2】绘制实值矩阵各元素对应的色块,通过不同的&参数/参数值&控制色块的形状、大小、颜色等属性。
方形色块,充满方格,灰白色字体
&& matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'TextColor',[0.6,0.6,0.6],'ColorBar','on');
方形色块,根据矩阵元素值自动确定色块的大小和颜色,不显示字体
&& matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'DisplayOpt','off','FigSize','Auto','ColorBar','on');
椭圆形色块,根据矩阵元素值自动确定色块的大小和颜色,不显示字体
&& matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'DisplayOpt','off','FigSize','Auto','ColorBar','on','FigShape','e');
圆形色块,根据矩阵元素值自动确定色块的大小和颜色,不显示字体,上三角形式显示
&& matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'DisplayOpt','off','FigSize','Auto','ColorBar','on','FigShape','c','FigStyle','Triu');
六边形色块,根据矩阵元素值自动确定色块的大小和颜色,显示字体,下三角形式显示
&& matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'DisplayOpt','on','FigSize','Auto','ColorBar','on','FigShape','h','FigStyle','Tril');
表盘形色块,根据矩阵元素值自动确定色块的颜色,不显示字体
&& matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'DisplayOpt','off','FigSize','Full','ColorBar','on','FigShape','d');
【例3】读取真彩图片,绘制伪红外图。
&& I = double(imread('kids.tif'));&& matrixplot(I,'DisplayOpt','off','colorbar','on','grid','off');
最后贴上matrixplot函数的源码,源码中有详细的调用格式及参数说明,感兴趣的版友可以通过设置不同参数,做出更为精彩的图形,希望大家喜欢!function matrixplot(data,varargin)%&& 根据实值矩阵绘制色块图,用丰富的颜色和形状形象的展示矩阵元素值的大小。%%&& matrixplot(data) 绘制矩阵色块图,data为实值矩阵,每一个元素对应一个色块,色%&&&&&&&&&&&&&&&&&&& 块颜色由元素值大小决定。%%&& matrixplot(data, 'PARAM1',val1, 'PARAM2',val2, ...) %&&&&&&&&& 用成对出现的参数名/参数值控制色块的各项属性。可用的参数名/参数值如下:%&&&&&&&&& 'FigShape' --- 设定色块的形状,其参数值为:%&&&&&&&&&&&&&&& 'Square'& --- 方形(默认)%&&&&&&&&&&&&&&& 'Circle'& --- 圆形%&&&&&&&&&&&&&&& 'Ellipse' --- 椭圆形%&&&&&&&&&&&&&&& 'Hexagon' --- 六边形%&&&&&&&&&&&&&&& 'Dial'&&& --- 表盘形%%&&&&&&&&& 'FigSize' --- 设定色块的大小,其参数值为:%&&&&&&&&&&&&&&& 'Full'&&& --- 最大色块(默认)%&&&&&&&&&&&&&&& 'Auto'&&& --- 根据矩阵元素值自动确定色块大小%%&&&&&&&&& 'FigStyle' --- 设定矩阵图样式,其参数值为:%&&&&&&&&&&&&&&& 'Auto'&&& --- 矩形矩阵图(默认)%&&&&&&&&&&&&&&& 'Tril'&&& --- 下三角矩阵图%&&&&&&&&&&&&&&& 'Triu'&&& --- 上三角矩阵图%%&&&&&&&&& 'FillStyle' --- 设定色块填充样式,其参数值为:%&&&&&&&&&&&&&&& 'Fill'&&& --- 填充色块内部(默认)%&&&&&&&&&&&&&&& 'NoFill'& --- 不填充色块内部%%&&&&&&&&& 'DisplayOpt' --- 设定是否在色块中显示矩阵元素值,其参数值为:%&&&&&&&&&&&&&&& 'On'&&&&& --- 显示矩阵元素值(默认)%&&&&&&&&&&&&&&& 'Off'&&&& --- 不显示矩阵元素值%%&&&&&&&&& 'TextColor' --- 设定文字的颜色,其参数值为:%&&&&&&&&&&&&&&& 表示单色的字符('r','g','b','y','m','c','w','k'),默认为黑色%&&&&&&&&&&&&&&& 1行3列的红、绿、蓝三元色灰度值向量([r,g,b])%&&&&&&&&&&&&&&& 'Auto'&&& --- 根据矩阵元素值自动确定文字颜色%%&&&&&&&&& 'XVarNames' --- 设定X轴方向需要显示的变量名(默认为X1,X2,...),其参数值为:%&&&&&&&&&&&&&&& 字符串矩阵或字符串元胞数组,若为字符串矩阵,其行数应与data的列数相同%&&&&&&&&&&&&&&& 若为字符串元胞数组,其长度应与data的列数相同。%%&&&&&&&&& 'YVarNames' --- 设定Y轴方向需要显示的变量名(默认为Y1,Y2,...),其参数值为:%&&&&&&&&&&&&&&& 字符串矩阵或字符串元胞数组,若为字符串矩阵,其行数应与data的行数相同%&&&&&&&&&&&&&&& 若为字符串元胞数组,其长度应与data的行数相同。%%&&&&&&&&& 'ColorBar' --- 设定是否显示颜色条,其参数值为:%&&&&&&&&&&&&&&& 'On'&&&&& --- 显示颜色条%&&&&&&&&&&&&&&& 'Off'&&&& --- 不显示颜色条(默认)%%&&&&&&&&& 'Grid' --- 设定是否显示网格线,其参数值为:%&&&&&&&&&&&&&&& 'On'&&&&& --- 显示网格线(默认)%&&&&&&&&&&&&&&& 'Off'&&&& --- 不显示网格线%%&& Example:%&& x = [1,-0.2,0.3,0.8,-0.5%&&&&&&& -0.2,1,0.6,-0.7,0.2%&&&&&&&& 0.3,0.6,1,0.5,-0.3%&&&&&&&& 0.8,-0.7,0.5,1,0.7%&&&&&&& -0.5,0.2,-0.3,0.7,1];%&& matrixplot(x);%&& matrixplot(x,'DisplayOpt','off');%&& matrixplot(x,'FillStyle','nofill','TextColor','Auto');%&& matrixplot(x,'TextColor',[0.7,0.7,0.7],'FigShap','s','FigSize','Auto','ColorBar','on');%&& matrixplot(x,'TextColor','k','FigShap','d','FigSize','Full','ColorBar','on','FigStyle','Triu');%&& XVarNames = {'xiezhh','heping','keda','tust','tianjin'};%&& matrixplot(x,'FigShap','e','FigSize','Auto','ColorBar','on','XVarNames',XVarNames,'YVarNames',XVarNames);%%&& CopyRight:xiezhh(谢中华),编写
% 对第一个输入参数类型进行判断if ~ismatrix(data) || ~isreal(data)&&& error('输入参数类型不匹配:第一个输入参数应为实值矩阵');end
% 解析成对出现的参数名/参数值[FigShape,FigSize,FigStyle,FillStyle,DisplayOpt,TextColor,XVarNames,...&&& YVarNames,ColorBar,GridOpt] = parseInputs(varargin{:});
% 产生网格数据[m,n] = size(data);[x,y] = meshgrid(0:n,0:m);data = data(:);maxdata = nanmax(data);mindata = nanmin(data);rangedata = maxdata -if isnan(rangedata)&&& warning('MATLAB:warning1','请检查您输入的矩阵是否合适!');&&&endz = zeros(size(x))+0.2;sx = x(1:end-1,1:end-1)+0.5;sy = y(1:end-1,1:end-1)+0.5;
if strncmpi(FigStyle,'Tril',4)&&& z(triu(ones(size(z)),2)&0) = NaN;&&& sx(triu(ones(size(sx)),1)&0) = NaN;elseif strncmpi(FigStyle,'Triu',4)&&& z(tril(ones(size(z)),-2)&0) = NaN;&&& sx(tril(ones(size(sx)),-1)&0) = NaN;endsx = sx(:);sy = sy(:);id = isnan(sx) | isnan(data);sx(id) = [];sy(id) = [];data(id) = [];
if isempty(XVarNames)&&& XVarNames = strcat('X',cellstr(num2str((1:n)')));else&&& if (iscell(XVarNames) && (numel(XVarNames) ~= n)) || (~iscell(XVarNames) && (size(XVarNames,1) ~= n))&&&&&&& error('X轴方向变量名应为字符串矩阵或字符串元胞数组,其长度与输入矩阵的列数相同');&&& endendif isempty(YVarNames)&&& YVarNames = strcat('Y',cellstr(num2str((1:m)')));else&&& if (iscell(YVarNames) && (numel(YVarNames) ~= m)) || (~iscell(YVarNames) && (size(YVarNames,1) ~= m))&&&&&&& error('Y轴方向变量名应为字符串矩阵或字符串元胞数组,其长度与输入矩阵的行数相同');&&& endend
% 绘图figure('color','w',...&&& 'units','normalized',...&&& 'pos',[0....68099]);axes('units','normalized','pos',[0.1,0.022,0.89,0.85]);if strncmpi(GridOpt,'On',2)&&& mesh(x,y,z,...&&&&&&& 'EdgeColor',[0.7,0.7,0.7],...&&&&&&& 'FaceAlpha',0,...&&&&&&& 'LineWidth',1);&& % 参考网格线endaxis([-0.1,n+0.1,-0.1,m+0.1,-0.5,0.5]);view(2);% 设置X轴和Y轴刻度位置及标签set(gca,'Xtick',(1:n)-0.5,...&&& 'XtickLabel',XVarNames,...&&& 'Ytick',(1:m)-0.5,...&&& 'YtickLabel',YVarNames,...&&& 'XAxisLocation','top',...&&& 'YDir','reverse',...&&& 'Xcolor',[0.7,0.7,0.7],...&&& 'Ycolor',[0.7,0.7,0.7],...&&& 'TickLength',[0,0]);axis off
% 绘制填充色块if strncmpi(FillStyle,'Fill',3)&&& MyPatch(sx',sy',data',FigShape,FigSize);end
% 显示数值文本信息if strncmpi(DisplayOpt,'On',2)&&& str = num2str(data,'%4.2f');&&& scale = 0.1*max(n/m,1)/(max(m,n)^0.55);&&& if strncmpi(TextColor,'Auto',3)&&n...
阅读(10377)|
人人移动客户端下载苹果/安卓/wp
积分 539, 距离下一级还需 261 积分
权限: 自定义头衔, 签名中使用图片
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡, 显身卡, 匿名卡, 抢沙发下一级可获得
权限: 隐身
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
现有一个156*4的矩阵,现在想通过for循环语句得到该矩阵中每12行中去掉第一行的数据从而组成新矩阵。比如原矩阵是156行4列,那么想得到的新矩阵为1:156中的2:12,14:24,26:36......,列数保持不变,跪求高手帮忙!
载入中......
也可以先用for,再嵌套if语句取余数判别,逢mod(m,12)==1就跳过
m=size(A,1);
n=size(A,2);
if mod(i,12)~=1
B(k,1:n)=A(i,1:n);
本帖最后由 遥远的生命 于
16:37 编辑
m=size(A,1) ;& && && && && && &&&%为待处理矩阵A的行数
n=size(A,2);& && && && && && && &%为待处理矩阵A的列数
if mod(m,12)==0& && && && &&&%判别行数是否是12的倍数
c=m/12 ;& && && && && && && && & %如果是12的倍数,进行k次循环
c=floor(m/12);& && && && && && &%如果不是12的倍数,向下取整数
cc=m-12*c;& && && && && && && &%最后剩cc行的
for i=0:c-1& && && && && && && &&&%前面12*c行的筛选
for j=2:12
B(k,1:n)=A(i*12+j,1:n);
if cc&1& && && && && && && && && && && && & %后面剩下的cc行筛选
for f=2:cc
B(k,1:n)=A(c*12+f,1:n);k=k+1;
好的意见建议
总评分:&热心指数 + 1&
也可以先用for,再嵌套if语句取余数判别,逢mod(m,12)==1就跳过
m=size(A,1);
n=size(A,2);
if mod(i,12)~=1
B(k,1:n)=A(i,1:n);
热心帮助其他会员
总评分:&学术水平 + 1&
热心指数 + 1&
信用等级 + 1&
a = randi(100,);
= size(a);
for st = 2:12:m
& & newa = ;
newa就是你需要的矩阵
无限扩大经管职场人脉圈!每天抽选10位免费名额,现在就扫& 论坛VIP& 贵宾会员& 可免费加入
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向,请联系(010-);
邮箱:service@pinggu.org
投诉或不良信息处理:(010-)
京ICP证090565号
论坛法律顾问:王进律师本帖子已过去太久远了,不再提供回复功能。MATLAB矩阵及其数值运算 - Matlab论坛 -
中国电子技术论坛 -
最好最受欢迎电子论坛!
后使用快捷导航没有帐号?
MATLAB矩阵及其数值运算
资深工程师
10:16:23  
MATLAB矩阵及其数值运算
变量和数据操作
变量与赋值
1.变量命名
在MATLAB6.5中,变量名是以字母开头,后接字母、数字或下划线的字符序列,最多63个字符。在MATLAB中,变量名区分字母的大小写。
2.赋值语句
(1) 变量=表达式
(2) 表达式
其中表达式是用运算符将有关运算量连接起来的式子,其结果是一个矩阵。
例2-1 计算表达式的值,并显示计算结果。
在MATLAB命令窗口输入命令:
y=3-sqrt(17);
z=(cos(abs(x+y))-sin(78*pi/180))/(x+abs(y))
其中pi和i都是MATLAB预先定义的变量,分别代表代表圆周率π和虚数单位。
输出结果是:
z = -0.3488 + 0.3286i
预定义变量
在MATLAB工作空间中,还驻留几个由系统本身定义的变量。例如,用pi表示圆周率π的近似值,用i,j表示虚数单位。
预定义变量有特定的含义,在使用时,应尽量避免对这些变量重新赋值。
内存变量的管理
1.内存变量的删除与修改
MATLAB工作空间窗口专门用于内存变量的管理。在工作空间窗口中可以显示所有内存变量的属性。当选中某些变量后,再单击Delete按钮,就能删除这些变量。当选中某些变量后,再单击Open按钮,将进入变量编辑器。通过变量编辑器可以直接观察变量中的具体元素,也可修改变量中的具体元素。
clear命令用于删除MATLAB工作空间中的变量。who和whos这两个命令用于显示在MATLAB工作空间中已经驻留的变量名清单。who命令只显示出驻留变量的名称,whos在给出变量名的同时,还给出它们的大小、所占字节数及数据类型等信息。
2.内存变量文件
利用MAT文件可以把当前MATLAB工作空间中的一些有用变量长久地保留下来,扩展名是.mat。MAT文件的生成和装入由save和load命令来完成。常用格式为:
save 文件名 [变量名表]&&[-append][-ascii]
load 文件名&&[变量名表] [-ascii]
其中,文件名可以带路径,但不需带扩展名.mat,命令隐含一定对.mat文件进行操作。变量名表中的变量个数不限,只要内存或文件中存在即可,变量名之间以空格分隔。当变量名表省略时,保存或装入全部变量。-ascii选项使文件以ASCII格式处理,省略该选项时文件将以二进制格式处理。save命令中的-append选项控制将变量追加到MAT文件中。
MATLAB常用数学函数
MATLAB提供了许多数学函数,函数的自变量规定为矩阵变量,运算法则是将函数逐项作用于矩阵的元素上,因而运算的结果是一个与自变量同维数的矩阵。
函数使用说明:
(1) 三角函数以弧度为单位计算。
(2) abs函数可以求实数的绝对值、复数的模、字符串的ASCII码值。
(3)用于取整的函数有fix、floor、ceil、round,要注意它们的区别。
(4)rem与mod函数的区别。rem(x,y)和mod(x,y)要求x,y必须为相同大小的实矩阵或为标量。
&&数据的输出格式
& &MATLAB用十进制数表示一个常数,具体可采用日常记数法和科学记数法两种表示方法。
& &在一般情况下,MATLAB内部每一个数据元素都是用双精度数来表示和存储的。数据输出时用户可以用format命令设置或改变数据输出格式。format命令的格式为:
& &format&&格式符
其中格式符决定数据的输出格式
MATLAB矩阵
&&矩阵的建立
1.直接输入法
& &最简单的建立矩阵的方法是从键盘直接输入矩阵的元素。具体方法如下:将矩阵的元素用方括号括起来,按矩阵行的顺序输入各元素,同一行的各元素之间用空格或逗号分隔,不同行的元素之间用分号分隔。
2.利用M文件建立矩阵
& &对于比较大且比较复杂的矩阵,可以为它专门建立一个M文件。下面通过一个简单例子来说明如何利用M文件创建矩阵。
& &例2-2&&利用M文件建立MYMAT矩阵。
(1) 启动有关编辑程序或MATLAB文本编辑器,并输入待建矩阵:
(2) 把输入的内容以纯文本方式存盘(设文件名为mymatrix.m)。
(3)在MATLAB命令窗口中输入mymatrix,即运行该M文件,就会自动建立一个名为MYMAT的矩阵,可供以后使用。
3.利用冒号表达式建立一个向量
& &冒号表达式可以产生一个行向量,一般格式是:
& &e1:e2:e3
其中e1为初始值,e2为步长,e3为终止值。
在MATLAB中,还可以用linspace函数产生行向量。其调用格式为:
linspace(a,b,n)
其中a和b是生成向量的第一个和最后一个元素,n是元素总数。
显然,linspace(a,b,n)与a:(b-a)/(n-1):b等价。
4.建立大矩阵
大矩阵可由方括号中的小矩阵或向量建立起来。
&&矩阵的拆分
1.矩阵元素
& &通过下标引用矩阵的元素,例如
A(3,2)=200
采用矩阵元素的序号来引用矩阵元素。矩阵元素的序号就是相应元素在内存中的排列顺序。在MATLAB中,矩阵元素按列存储,先第一列,再第二列,依次类推。例如
A=[1,2,3;4,5,6];
显然,序号(Index)与下标(Subscript)是一一对应的,以m×n矩阵A为例,矩阵元素A(i,j)的序号为(j-1)*m+i。其相互转换关系也可利用sub2ind和ind2sub函数求得。
2.矩阵拆分
& & (1)利用冒号表达式获得子矩阵
& && & ①A(:,j)表示取A矩阵的第j列全部元素;A(i,:)表示A矩阵第i行的全部元素;A(i,j)表示取A矩阵第i行、第j列的元素。
& && &②A(i:i+m,:)表示取A矩阵第i~i+m行的全部元素;A(:,k:k+m)表示取A矩阵第k~k+m列的全部元素,A(i:i+m,k:k+m)表示取A矩阵第i~i+m行内,并在第k~k+m列中的所有元素。
此外,还可利用一般向量和end运算符来表示矩阵下标,从而获得子矩阵。end表示某一维的末尾元素下标。
& & (2)利用空矩阵删除矩阵的元素
& &在MATLAB中,定义[]为空矩阵。给变量X赋空矩阵的语句为X=[]。注意,X=[]与clearX不同,clear是将X从工作空间中删除,而空矩阵则存在于工作空间中,只是维数为0。
&&特殊矩阵
1.通用的特殊矩阵
常用的产生通用特殊矩阵的函数有:
zeros:产生全0矩阵(零矩阵)。
ones:产生全1矩阵(幺矩阵)。
eye:产生单位矩阵。
rand:产生0~1间均匀分布的随机矩阵。
randn:产生均值为0,方差为1的标准正态分布随机矩阵。
例2-3 分别建立3×3、3×2和与矩阵A同样大小的零矩阵。
(1) 建立一个3×3零矩阵。
(2) 建立一个3×2零矩阵。
zeros(3,2)
(3)设A为2×3矩阵,则可以用zeros(size(A))建立一个与矩阵A同样大小零矩阵。
A=[1 2 3;4 56];& &%产生一个2×3阶矩阵A
zeros(size(A))& &&&%产生一个与矩阵A同样大小的零矩阵
例2-4 建立随机矩阵:
(1) 在区间[20,50]内均匀分布的5阶随机矩阵。
(2) 均值为0.6、方差为0.1的5阶正态分布随机矩阵。
命令如下:
x=20+(50-20)*rand(5)
y=0.6+sqrt(0.1)*randn(5)
此外,常用的函数还有reshape(A,m,n),它在矩阵总元素保持不变的前提下,将矩阵A重新排成m×n的二维矩阵。
2.用于专门学科的特殊矩阵
(1) 魔方矩阵
魔方矩阵有一个有趣的性质,其每行、每列及两条对角线上的元素和都相等。对于n阶魔方阵,其元素由1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方矩阵的函数magic(n),其功能是生成一个n阶魔方阵。
例2-5 将101~125等25个数填入一个5行5列的表格中,使其每行每列及对角线的和均为565。
M=100+magic(5)
(2) 范得蒙矩阵
范得蒙(Vandermonde)矩阵最后一列全为1,倒数第二列为一个指定的向量,其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在MATLAB中,函数vander(V)生成以向量V为基础向量的范得蒙矩阵。例如,A=vander([1;2;3;5])即可得到上述范得蒙矩阵。
(3) 希尔伯特矩阵
在MATLAB中,生成希尔伯特矩阵的函数是hilb(n)。
使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果。MATLAB中,有一个专门求希尔伯特矩阵的逆的函数invhilb(n),其功能是求n阶的希尔伯特矩阵的逆矩阵。
例2-6 求4阶希尔伯特矩阵及其逆矩阵。
命令如下:
formatrat& & %以有理形式输出
H=invhilb(4)
(4) 托普利兹矩阵
托普利兹(Toeplitz)矩阵除第一行第一列外,其他每个元素都与左上角的元素相同。生成托普利兹矩阵的函数是toeplitz(x,y),它生成一个以x为第一列,y为第一行的托普利兹矩阵。这里x,y均为向量,两者不必等长。toeplitz(x)用向量x生成一个对称的托普利兹矩阵。例如
T=toeplitz(1:6)
(5) 伴随矩阵
MATLAB生成伴随矩阵的函数是compan(p),其中p是一个多项式的系数向量,高次幂系数排在前,低次幂排在后。例如,为了求多项式的x3-7x+6的伴随矩阵,可使用命令:
p=[1,0,-7,6];
(6) 帕斯卡矩阵
我们知道,二次项(x+y)n展开后的系数随n的增大组成一个三角形表,称为杨辉三角形。由杨辉三角形表组成的矩阵称为帕斯卡(Pascal)矩阵。函数pascal(n)生成一个n阶帕斯卡矩阵。
例2-7 求(x+y)5的展开式。
在MATLAB命令窗口,输入命令:
矩阵次对角线上的元素1,5,10,10,5,1即为展开式的系数。
MATLAB运算
1.基本算术运算
& &MATLAB的基本算术运算有:+(加)、-(减)、*(乘)、/(右除)、\(左除)、^(乘方)。
注意,运算是在矩阵意义下进行的,单个数据的算术运算只是一种特例。
& & (1)矩阵加减运算
& &假定有两个矩阵A和B,则可以由A+B和A-B实现矩阵的加减运算。运算规则是:若A和B矩阵的维数相同,则可以执行矩阵的加减运算,A和B矩阵的相应元素相加减。如果A与B的维数不相同,则MATLAB将给出错误信息,提示用户两个矩阵的维数不匹配。
& & (2)矩阵乘法
假定有两个矩阵A和B,若A为m×n矩阵,B为n×p矩阵,则C=A*B为m×p矩阵。
& & (3)矩阵除法
在MATLAB中,有两种矩阵除法运算:\和/,分别表示左除和右除。如果A矩阵是非奇异方阵,则A\B和B/A运算可以实现。A\B等效于A的逆左乘B矩阵,也就是inv(A)*B,而B/A等效于A矩阵的逆右乘B矩阵,也就是B*inv(A)。
对于含有标量的运算,两种除法运算的结果相同,如3/4和4\3有相同的值,都等于0.75。又如,设a=[10.5,25],则a/5=5\a=[2.]。对于矩阵来说,左除和右除表示两种不同的除数矩阵和被除数矩阵的关系。对于矩阵运算,一般A\B≠B/A。
& & (4)矩阵的乘方
& &一个矩阵的乘方运算可以表示成A^x,要求A为方阵,x为标量。
& &在MATLAB中,有一种特殊的运算,因为其运算符是在有关算术运算符前面加点,所以叫点运算。点运算符有.*、./、.\和.^。两矩阵进行点运算是指它们的对应元素进行相关运算,要求两矩阵的维参数相同。
&&关系运算
& &MATLAB提供了6种关系运算符:&(小于)、&=(小于或等于)、&(大于)、&=(大于或等于)、==(等于)、~=(不等于)。它们的含义不难理解,但要注意其书写方法与数学中的不等式符号不尽相同。
& &关系运算符的运算法则为:
& & (1)当两个比较量是标量时,直接比较两数的大小。若筽arty闪ⅲ?叵当泶锸浇峁?
10:16 上传
点击文件名下载附件
下载积分: 积分 -1 分
284.5 KB, 下载次数: 38, 下载积分: 积分 -1 分
14:06:02  
下载了打不开呀。。求高手助我
12:14:50  
19:44:29  
很好 支持一下
21:02:13  
很好,支持一下
FPGA是当前很热门的技术。
发展到现在,FPGA器件由早期的纯逻辑粘合发展到如今的可编程片上系统(SOC),FPGA应用的领域愈发广泛,开发难度和复杂度也越来越大。
那么究竟如何才能高效学习好FPGA技术呢?
每天选一个需要解决的问题,大家一起来帮忙。
授人玫瑰,手有余香
USB Type-C终结了长期以来USB需要通过试错法来多次重插的历史,一个接口搞定了电能、数据、音视频数据三种传输需求,形成接口和电缆以及快速充电协议的大统一。
预计2017年全球所有新发布手机都将采用USB Type-C接口和USB PD供电协议,将将引发手机行业革命。
Powered by

我要回帖

更多关于 matlab生成对称矩阵 的文章

 

随机推荐