编程 – 第 5 页

编程：MATLAB实现粒子二维随机运动轨迹

By: taho

On: 2012年12月26日

代码如下：

clear;
n=10;%粒子数
step=200; %粒子前进的步数
jump=0.1; %每前进一步粒子的跳跃位置，这里认为是一个自然数，甚至程序中给出小数也可以。

figure; %创建绘图窗口
hold on; %打开重绘
grid; %显示网格

%以下是为了显示两个上下极限概率的边界
x_line=(1:step);
y_line_upper=jump*x_line; %上边界的y值
y_line_lower=-jump*x_line; %下边界的y值
plot(x_line,y_line_upper); %绘制上边界
plot(x_line,y_line_lower); %绘制下边界
clear x_line;clear y_line_upper;clear y_line_lower; %清除这几个变量，节省内存

for m=1:n; %粒子数(n)个循环
move_x=zeros(100); %初始化粒子的横坐标位置数组
move_y=zeros(100); %初始化粒子的纵坐标位置数组

y=0; %初始位置的纵坐标
for x=1:step; %粒子随机前进500步
c=rand(1); %生成随机数
%以下生成的是非等概率分布
if c<=0.20
a=-jump; %概率分布为1/4
elseif c>0.10 && c<0.70
a=0; %概率分布为1/2
else
a=jump; %概率分布为1/4
end

y=y+a;%纵坐标y改变，横坐标是步长j
move_x(x)=x; %把当前步长的粒子位置横坐标给位置数组
move_y(x)=y; %把当前步长的粒子位置纵坐标给位置数组
plot(x,y,'.'); %绘当前位置点
pause(0.00001);%延时绘图，动画显示
end
plot(move_x,move_y,'r');%将点连线，形成粒子的运动轨迹
end
clear all; %清除所有变量

clear;

n=10;%粒子数

step=200; %粒子前进的步数

jump=0.1; %每前进一步粒子的跳跃位置，这里认为是一个自然数，甚至程序中给出小数也可以。

figure; %创建绘图窗口

hold on; %打开重绘

grid; %显示网格

%以下是为了显示两个上下极限概率的边界

x_line=(1:step);

y_line_upper=jump*x_line; %上边界的y值

y_line_lower=-jump*x_line; %下边界的y值

plot(x_line,y_line_upper); %绘制上边界

plot(x_line,y_line_lower); %绘制下边界

clear x_line;clear y_line_upper;clear y_line_lower; %清除这几个变量，节省内存

for m=1:n; %粒子数(n)个循环

move_x=zeros(100); %初始化粒子的横坐标位置数组

move_y=zeros(100); %初始化粒子的纵坐标位置数组

y=0; %初始位置的纵坐标

for x=1:step; %粒子随机前进500步

c=rand(1); %生成随机数

%以下生成的是非等概率分布

if c<=0.20

a=-jump; %概率分布为1/4

elseif c>0.10 && c<0.70

a=0; %概率分布为1/2

else

a=jump; %概率分布为1/4

end

y=y+a;%纵坐标y改变，横坐标是步长j

move_x(x)=x; %把当前步长的粒子位置横坐标给位置数组

move_y(x)=y; %把当前步长的粒子位置纵坐标给位置数组

plot(x,y,'.'); %绘当前位置点

pause(0.00001);%延时绘图，动画显示

end

plot(move_x,move_y,'r');%将点连线，形成粒子的运动轨迹

end

clear all; %清除所有变量

视频演示：待更新 &n阅读全文

图书：几本MATLAB函数查询手册

By: taho

On: 2012年12月25日

In: Matlab, 读书

With: 0 Comments

MATLAB函数.chm MATLAB函数速查.pdf MATLAB函数查询手册官方.pdf &nb阅读全文

编程：怎样用MATLAB制作动画

By: taho

On: 2012年12月25日

In: Matlab

With: 0 Comments

电影动画的好处就是，运行一次可以多次播放，甚至可以直接生成avi文件，直接独立与Matlab环境播放阅读全文

编程：用Matlab画圆锥的两种方法

By: taho

On: 2012年12月18日

In: Matlab

With: 0 Comments

%第一种方法，结果见图“figure1” figure; [X,Y,Z]=cylinder(0:25阅读全文

编程：MATLAB多行注释的三种方法，取整函数，eps

By: taho

On: 2012年12月11日

In: Matlab

With: 0 Comments

多行注释的三种方法 1. %{ 若干语句 %} 2. 多行注释: 选中要注释的若干语句, 编辑器菜单阅读全文

编程：MATLAB中先指定概率分布再产生随机数的方法

By: taho

On: 2012年12月10日

In: Matlab

With: 0 Comments

我的应用例子可以到本站查看一个我使用的一个代码：【编程：MATLAB实现粒子二维随机运动轨迹】产阅读全文

编程：C++输入输出流

By: taho

On: 2012年12月4日

In: VC++

With: 0 Comments

阅读全文

软件：PSPLOT——G77的绘图工具

By: taho

On: 2012年11月23日

In: Fortran, 计算物理

With: 0 Comments

Scientific Graphics using Fortran ——&#阅读全文

软件：G77的万能式安装使用和说明文件

By: taho

On: 2012年11月21日

In: Fortran, 计算物理

With: 0 Comments

我修改了一点G77，安装非常简单，可以任意安装，随你喜好！请使用这个版本的G77，下载地址：G77阅读全文

电脑技巧：如何使用CMD修改系统path环境变量、窗口和字体颜色

By: taho

On: 2012年11月21日

In: 电脑技术, 编程

With: 0 Comments

CMD获取当前目录的绝对路径 @echo off echo 当前盘符：%~d0 echo 当前盘符和阅读全文

技术：G77与MATLAB接口和FORTRAN77语法概述/简单教程

By: taho

On: 2012年11月15日

In: Fortran, Matlab, 计算物理

With: 0 Comments

MATLAB与C++、G77的混合编程《深入浅出MATLAB 7.X混合编程》董维国编著电子阅读全文

转载：技术文章——C++与C的区别

By: taho

On: 2012年11月5日

In: VC++, 转载和摘抄

With: 0 Comments

相对来说，C语言是一种简洁的语言，所涉及的概念和元素比较少，主要是：宏（macro）、指针（poin阅读全文

资料：Windows Media编程

By: taho

On: 2012年10月31日

In: VC++

With: 0 Comments

《Windows Media编程与应用》 http://book.douban.com/subjec阅读全文

编程：VC++ 类的继承，构造函数，数据初始化

By: taho

On: 2012年10月29日

In: VC++

With: 0 Comments

/*************************************************

功能：类的继承，构造函数，数据初始化

作者：TAHO

时间：2012年10月29日

**************************************************/

#include <iostream.h>

class vehicle
{
protected:
int wheels,weight;
public:
vehicle(int i,int j)
{
wheels=i;
weight=j;
}
int get2ws() //这里的return功能不能被写进构造函数vehicle(…)里，否则不被执行，构造函数只负责数据初始化，不负责返回数值等操作。
{
return wheels;
return weight;
}
};

class car:private vehicle //私有继承
{
public:
int passenger_load;
car(int i,int j,int p_l):vehicle(i,j) //构造函数，有i,j，初始化了基类的成员i,j
{
passenger_load=p_l;
cout<<“The wheels number of car is: ”<<wheels<<endl;
cout<<“The weight of car is: ”<<weight<<endl;
cout<<“The passenger load of car is ”<<passenger_load<<“nn”;
}
};

class truck:private vehicle //私有继承
{
public:
int passenger_load,payload;
truck(int i,int j,int p_l,int pay_l):vehicle(i,j) //构造函数，有i,j，又一次初始化了基类的成员i,j，旧的i,j值被替换
{
passenger_load=p_l;
payload=pay_l;
cout<<“The wheels number of truck is: ”<<wheels<<endl;
cout<<“The weight of truck is: ”<<weight<<endl;
cout<<“The passenger load of truck is ”<<passenger_load<<endl;
cout<<“The payload of truck is ”<<payload<<“nn”;
}
};

void main()
{
car car1(4,10,6);
truck truck1(6,25,2,40);
}

/*************************************************

功能：类的继承，构造函数，数据初始化

作者：TAHO

时间：2012年10月29日

**************************************************/

#include <iostream.h>

class vehicle

{

protected:

int wheels,weight;

public:

vehicle(int i,int j)

{

wheels=i;

weight=j;

}

int get2ws() //这里的return功能不能被写进构造函数vehicle(…)里，否则不被执行，构造函数只负责数据初始化，不负责返回数值等操作。

{

return wheels;

return weight;

}

};

class car:private vehicle //私有继承

{

public:

int passenger_load;

car(int i,int j,int p_l):vehicle(i,j) //构造函数，有i,j，初始化了基类的成员i,j

{

passenger_load=p_l;

cout<<“The wheels number of car is: ”<<wheels<<endl;

cout<<“The weight of car is: ”<<weight<<endl;

cout<<“The passenger load of car is ”<<passenger_load<<“nn”;

}

};

class truck:private vehicle //私有继承

{

public:

int passenger_load,payload;

truck(int i,int j,int p_l,int pay_l):vehicle(i,j) //构造函数，有i,j，又一次初始化了基类的成员i,j，旧的i,j值被替换

{

passenger_load=p_l;

payload=pay_l;

cout<<“The wheels number of truck is: ”<<wheels<<endl;

cout<<“The weight of truck is: ”<<weight<<endl;

cout<<“The passenger load of truck is ”<<passenger_load<<endl;

cout<<“The payload of truck is ”<<payload<<“nn”;

}

};

void main()

{

car car1(4,10,6);

truck truck1(6,25,2,40);

}

阅读全文