sys 数值型和符号型在MATLAB时不同的。用sym转化
矩阵
1.全零阵
zeros(n)
B=zeros(size(A)) 生成与矩阵A相同大小的全零阵
2.单位阵
y=eye(n)
Y=eye(size(A))
3.全1阵
Y=ones(n)
4.均匀分布随机矩阵呢个的生成
rand
Y=rand(n) 生成n*n的随机矩阵,元素(0,1)
s=rand('state') 产生35个元素的向量
5.正态分布随机矩阵的生成
T=rand(n)
6.随机排序
函数randperm
格式 p=randperm(n)
7.产生线性等分向量(就是1*n的矩阵)
函数 linspace
y=linspace(a,b) 在(a,b)产生100个线性等分点
y=linspace(a,b,n) 生成n个线性等分点
8.产生对数等分向量
logspace
y=logspace(a,b) 在(10^a,10^b)产生50个等分向量
9.矩阵中元素的个数
numel
n=numel(A)
10.产生已输入为对角线元素的矩阵
blkding
11.友矩阵的生成 //???
compan(u) u为多项式系统向量,A的
12.hadamard矩阵的生成
hadamard
H=hadamard(n)
平面和曲面的分割图
clc,clear,close all
clf,a=-20;
eps0=1;%双精度
[x,y]=meshgrid(-10:0.2:10); %生成平面网格
v=[-10 10 -10 10 -100 100]; %设定空间坐标系的范围
%colormap(gray)
z1=(x.^2-2*y.^2)+eps; %马鞍面
z2=a*ones(size(x));
r0=abs(z1-z2)<=eps0; %EPS是MATLAB中的函数,表示浮点相对精度
zz=r0.*z2;xx=r0.*x;yy=r0.*y;
subplot(2,2,2), %在第2图形窗口
h1=plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'+');
set(h1,'markersize',2),hold on,
axis("vis3d"),
grid on,
subplot(2,2,1);
mesh(x,y,z1);
grid,
hold on;
mesh(x,y,z2);
%{
h2=plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'.');
set(h2,'markersize',6),hold on,axis(v),
for i=1:5
a=70-i*30;
z2=a*ones(size(x));
r0=abs(z1-z2)<=1;
zz=r0.*z2;
yy=r0.*y;
xx=r0.*x;
subplot(2,2,3),
mesh(x,y,z2);
grid,hold on;
mesh(x,y,z2);
hidden off %取消掩盖
h2=plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'.');
axis(v),grid
subplot(2,2,4),
h4=plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'o');
set(h4,'markersize',2),
hold on,
axis("vis3d"),
grid on
end
%}