原文地址：matlab分段函数绘图
1.一元分段函数绘图
例如：
把下面的函数保存为Piecewise_x.m文件
function F=Piecewise_x(x)
F=x.^2.*(x>=0 & x<1)+cos(pi*(x-1)).*(x>=1 & x<2)+(-x.^2./(x+2)).*(x>=2 & x<=4);
end
新浪博客
1.一元分段函数绘图
例如：
把下面的函数保存为Piecewise_x.m文件
function F=Piecewise_x(x)
F=x.^2.*(x>=0 & x<1)+cos(pi*(x-1)).*(x>=1 & x<2)+(-x.^2./(x+2)).*(x>=2 & x<=4);
end
运行：
x=linspace(0,4);
F=Piecewise_x(x);%计算相应函数值
plot(x,F);%绘制曲线
hold on;
plot(1*on
plot(2*on
2.二元分段函数绘图
例如：
把下面的函数保存为Piecewise_xy.m文件
function Pxy=Piecewise_xy(x,y)
Pxy=0.5457*exp(-0.75*y^2 - 3.75*x^2 - 1.5*x).*(x+y>1)+...
end
运行：
[x,y]=meshgrid(-3:0.1:3);
s=size(x)
Pxy=zeros(s(1),s(2));
for i=1:s(1)
end
mesh(x,y,Pxy)
3
举例如下：
在matlab中，>> x1=0:0.01:1;
4例2
0<t<1 u(t)=0
其他时候u(t)=1
这个用matlab语言怎么定义函数u
如下：
function myfun(t) if 0<t && t<1 u=0 else u=1 end 输入myfun(8),可以得到结果
例3
比如说y=5/cos(x),我想绘制此曲线.但是x的定义域中不能包括1/2*pi和3/2*pi,怎样才能定义x,可以在除x=1/2*pi和3/2*pi以外绘制区间为[0,2*pi]的曲线. 谢谢各位厉害的大虾帮忙. 尽量写得详细一些.
如下：
>> x=0:.001:2*pi; >> c=.1; >> y=5./cos(x); >> y(abs(x-pi/2)<c | abs(x-3*pi/2)<c)=NaN; >> plot(x,y);
例4 比较经典的一个分段函数绘图例子
x=-3:0.01:3;
y1=zeros(size(x));
y2=zeros(size(x));
y3=zeros(size(x));
N=length(x);
for k=1:N
if x(k)<-1&x(k)>=-3;
y1(k)=(-x(k).^2-4*x(k)-3)/2;
elseif x(k)>=-1&x(k)<1 ;
y2(k)=-x(k).^2+1;
else x(k)<=3&x(k)>=1 ;
y3(k)=(-x(k).^2+4*x(k)-3)/2;
end
end
y=y1+y2+y3;
plot(x,y)
x=linspace(0,4);
F=Piecewise_x(x);%计算相应函数值
plot(x,F);%绘制曲线
hold on;
plot(1*on
plot(2*on
2.二元分段函数绘图
例如：
把下面的函数保存为Piecewise_xy.m文件
function Pxy=Piecewise_xy(x,y)
Pxy=0.5457*exp(-0.75*y^2 - 3.75*x^2 - 1.5*x).*(x+y>1)+...
end
运行：
[x,y]=meshgrid(-3:0.1:3);
s=size(x)
Pxy=zeros(s(1),s(2));
for i=1:s(1)
end
mesh(x,y,Pxy)
3
举例如下：
在matlab中，>> x1=0:0.01:1;
4例2
0<t<1 u(t)=0
其他时候u(t)=1
这个用matlab语言怎么定义函数u
如下：
function myfun(t) if 0<t && t<1 u=0 else u=1 end 输入myfun(8),可以得到结果
例3
比如说y=5/cos(x),我想绘制此曲线.但是x的定义域中不能包括1/2*pi和3/2*pi,怎样才能定义x,可以在除x=1/2*pi和3/2*pi以外绘制区间为[0,2*pi]的曲线. 谢谢各位厉害的大虾帮忙. 尽量写得详细一些.
如下：
>> x=0:.001:2*pi; >> c=.1; >> y=5./cos(x); >> y(abs(x-pi/2)<c | abs(x-3*pi/2)<c)=NaN; >> plot(x,y);
例4 比较经典的一个分段函数绘图例子
x=-3:0.01:3;
y1=zeros(size(x));
y2=zeros(size(x));
y3=zeros(size(x));
N=length(x);
for k=1:N
if x(k)<-1&x(k)>=-3;
y1(k)=(-x(k).^2-4*x(k)-3)/2;
elseif x(k)>=-1&x(k)<1 ;
y2(k)=-x(k).^2+1;
else x(k)<=3&x(k)>=1 ;
y3(k)=(-x(k).^2+4*x(k)-3)/2;
end
end
y=y1+y2+y3;
plot(x,y)