三维曲线的绘制方法
1.plot3函数
plot3(x，y，z) 其中，参数x、y、z组成一组曲线的坐标。
例题：绘制一条空间折线。
x=[0.2,1.8,2.5];
y=[1.3,2.8,1.1];
z=[0.4,1.2,1.6];
plot3(x,y,z)
grid on;
axis([0,3,1,3,0,2]);
例题：绘制螺旋线
t=linspace(0,10*pi,200);
x=sin(t)+t.*cos(t);
y=cos(t)-t.*sin(t);
z=t;
subplot(1,2,1);
plot3(x,y,z)
grid on
subplot(1,2,2);
plot3(x(1:4:200),y(1:4:200),z(1:4:200))
grid on
 2.plot3函数参数的变化形式
plot3(x，y，z) 其中，参数x，y，z是同型矩阵。
方法一：
t=0:0.01:2*pi;
t=t';
x=[t,t,t,t,t];
y=[sin(t),sin(t)+1,sin(t)+2,sin(t)+3,sin(t)+4];
z=[t,t,t,t,t];
plot3(x,y,z)
方法二：
t=0:0.01:2*pi;
x=t';
y=[sin(t);sin(t)+1;sin(t)+2;sin(t)+3;sin(t)+4];
z=t';
plot3(x,y,z)
3.含多组输入参数的plot3函数
plot3(x1,y1,z1,x2,y2,……,xn,yn,zn)
每一组x，y，z向量构成一组数据点的坐标，绘制一条曲线。
例题：绘制三条不同长度的正弦曲线：
t1=0:0.01:1.5*pi;
t2=0:0.01:2*pi;
t3=0:0.01:3*pi;
plot3(t1,sin(t1),t1,t2,sin(t2)+1,t2,t3,sin(t3)+3,t3)
4.含选项的plot3函数
Plot3(x，y，z，选项) 选项用于说明曲线的线型、颜色和数据点标记。
例题：绘制空间曲线：
解答：
t=0:pi/50:6*pi;
x=cos(t);
y=sin(t);
z=2*t;
plot3(x,y,z,'p')
xlabel('X');ylabel('Y'),zlabel('Z');
grid on
5.fplot(funx，funy，funz，tlims)
其中，funx，funy，funz，tlims代表定义曲线x，y，z坐标的函数，通常采用函数句柄的形式。Lims为参数函数自变量的取值范围，用二元向量[tmin,tmax]描述，默认为[-5,5]
例题：绘制墨西哥帽顶曲线，曲线的参数方程如下：
解答：
xt=@(t)exp(-t/10).*sin(5*t);
yt=@(t)exp(-t/10).*cos(5*t); 
zt=@(t)t;
fplot3(xt,yt,zt,[-12,12])
在fplot3函数中，可以指定曲线的线型、颜色和数据点标记。
xt=@(t)exp(-t/10).*sin(5*t);
yt=@(t)exp(-t/10).*cos(5*t);
zt=@(t)t;
fplot3(xt,yt,zt,[-12,12],'r-.')

Java绘制不规则几何图形，比如划曲线，写字，线条随意画，如截图所示，甚至可以写出文字：
不规则图形的绘制代码：
public class IrregularShapeDemo extends JFrame {
GeneralPath gPath= new GeneralPath(); //GeneralPath对象实例
//构造函数
public IrregularShapeDemo() {
super("不规则图形的绘制"); //调用父类构造函数
enableEvents(AWTEvent.MOUSE_EVENT_MASK|AWTEvent.MOUSE_MOTION_EVENT_MASK); //允许事件
setSize(300, 200); //设置窗口尺寸
setVisible(true); //设置窗口可视
setDefaultCloseOperation(JFrame.EXIT_ON_CLOSE); //关闭窗口时退出程序
}
public void paint(Graphics g) { //重载窗口组件的paint()方法
Graphics2D g2D = (Graphics2D)g;//获取图形环境
g2D.draw(gPath); //绘制路径
}
public static void main(String[] args) {
new IrregularShapeDemo();
}
protected void processMouseEvent(MouseEvent e) { //鼠标事件处理
if(e.getID() == MouseEvent.MOUSE_PRESSED) {
aPoint = e.getPoint(); //得到当前鼠标点
gPath = new GeneralPath(); //重新实例化GeneralPath对象
gPath.moveTo(aPoint.x,aPoint.y); //设置路径点
}
}
protected void processMouseMotionEvent(MouseEvent e) { //鼠标运动事件处理
if(e.getID() == MouseEvent.MOUSE_DRAGGED) {
aPoint = e.getPoint(); //得到当前鼠标点
gPath.lineTo(aPoint.x, aPoint.y); //设置路径
gPath.moveTo(aPoint.x, aPoint.y);
repaint(); //重绘组件
}
}
}

一、模型的常微分方程及初值
二、ode45求解函数
function r=hudie3(t,x)
global a;
global b;
global c;
a=8/3;b=10;c=28;
r=zeros(3,1);
r(1)=-a*x(1)+x(2)*x(3);
r(2)=-b*x(2)+b*x(3);
r(3)=-x(1)*x(2)+c*x(2)-x(3);
end
三、绘制混沌系统三维曲线图
a=8/3;b=10;c=28;
t0=[0,100];f0=[0,0,1e-10];
[t,x]=ode45('hudie3',t0,f0);
plot3(x(:,1),x(:,2),x(:,3));
title('Lorenz模型');xlabel('x');ylabel('y');zlabel('z');
grid on
四、绘制混沌系统三个平面上的投影图
a=8/3;b=10;c=28;
t0=[0,100];
f0=[0,0,1e-10];
[t,x]=ode45('hudie3',t0,f0);
subplot(2,2,1)
plot3(x(:,1),x(:,2),x(:,3));
title('Lorenz模型');xlabel('x');ylabel('y');zlabel('z');
grid on
subplot(2,2,2)
plot(x(:,1),x(:,2));
title('Lorenz模型X-Y平面图');xlabel('x');ylabel('y');
grid on
subplot(2,2,3)
plot(x(:,1),x(:,3));
title('Lorenz模型X-Z平面图');xlabel('x');ylabel('z');
grid on
subplot(2,2,4)
plot(x(:,2),x(:,3));
title('Lorenz模型Y-Z平面图');xlabel('y');ylabel('z');
grid on
五、修改初值观察混沌系统图像的变化
以下是给出的几组任意的初值：
其图像如下：
六、探究各初值对图像的影响的参考代码
a=8/3;b=10;c=28;
for k=0:50:400
t0=[0,100];f0=[k,0,1e-10];
[t,x]=ode45('hudie3',t0,f0);
subplot(3,3,(k+50)/50)
plot3(x(:,1),x(:,2),x(:,3));
title('Lorenz模型');xlabel('x');ylabel('y');zlabel('z');
grid on
end
上面程序的图像如下:
探究第一个初值的变化对图像的影响的时候将变量K放在f0向量的第一个参数上，如上程序所示，探究第二个初值的变化对图像的影响的时候将变量K放在f0向量的第二个参数上，以此类推。
参数的取值范围通过修改 for k=0:50:400语句实现，但是注意在修改该语句时对应的 subplot(3,3,(k+50)/50)语句中(k+50)/50也要修改，比如将for k=0:50:400改为for k=0:500:4000，则subplot(3,3,(k+50)/50)应该相应的改为subplot(3,3,(k+500)/500)

前面介绍的余弦曲线的绘制，我们看到的是一个完整的静态图形，能否动态地展现绘制的过程？答案是肯定的，我们可以采用线程的方式来实现，参考代码如下：
import java.applet.Applet;
import java.awt.Color;
import java.awt.Graphics;
public class donghua_cos extends Applet implements Runnable
//通过实现Runnable接口实现线程操作
{
int x,y;
double a;
int xpos=0;
Thread runner;
boolean painted=false;
public void init() //Applet创建即启动执行，坐标初始化
{
// TODO Auto-generated method stub
Graphics g=getGraphics();  //画画之前，必须先取得画笔
for(x=0;x<=750;x+=1)    //画x轴
{
g.drawString("·",x,200);
if(x<=385) g.drawString("·",360,x);
}
g.drawString("Y",330,20);   //画y轴
for(x=360;x<=370;x+=1)   //画y轴箭头
{
g.drawString("·",x-10,375-x);
g.drawString("·",x,x-355);
}
g.drawString("X",735,230);
for(x=740;x<=750;x+=1)   //画x轴箭头
{
g.drawString("·",x,x-550);
g.drawString("·",x,950-x);
}
}
public void start() //Applet创建后自启动方法
{
// TODO Auto-generated method stub
if(runner==null){
runner=new Thread(this);   //通过Thread类来启动Runnable
runner.start();  //线程启动
}
}
public void stop() //Applet生命周期结束后自启动方法
{
// TODO Auto-generated method stub
if(runner!=null){
runner=null;//结束线程
}
}
public void run() //线程运行方法
{
// TODO Auto-generated method stub
while(true){
for(xpos=0;xpos<900-90;xpos+=3)
//循环设置曲线x轴坐标边界
{
repaint();//调用paint()方法
try{
Thread.sleep(100);//线程休息100毫秒
}catch(InterruptedException e){}
if(painted)
{
painted=false;
}
}
}
}
public void paint(Graphics g)//画图方法
{
for(x=0;x<=xpos;x+=1)//循环画曲线
{
a=Math.cos(x*Math. PI/180+Math.PI);
y=(int)(200+80*a);//放大80倍并向下平移200个像素
g.drawString("·",x,y);
}
painted=true;
}