WHOOL! fantastic world 工作,学习,生活

28Dec/090

King and Queens

3年前的一个深夜,看着queens的温布利演唱会,以为那是绝唱。Freddie
Mercury
的声音充满着奇幻的色彩。
昨晚,深夜,看着MJ的That is
it。再一次被震撼,这些年歌听的很多,不过给自己震撼的感觉没有比这两个更强悍了。
The King of POP and queens. I love you,always.
Tags:
Filed under: 碎碎念 No Comments
13Dec/090

HALTON序列

function result = halton( m,base,seeder )

%生成HALTON序列

% Check inputs
if nargin < 3
seeder = 0;
if nargin < 2

error(‘hATLAB:Halton:NotEnoughInputs’,…
‘Not enough input arguments. See Halton.’);

end
end
res=0;
n=length(base);
for i=1:m
for j=1:n
element=0;
temp=seeder+i;
k=1;
while
temp>0

element(k)=rem(temp,base(j));

Tags: ,
12Dec/090

Generating Correlated Random Numbers

This article describes common methods that are used in
generating correlated random numbers.
Contents
Generating two sequences of correlated random numbers
Generating two sequences of random numbers with a given
correlation is done in two simple steps:
Generate two sequences of uncorrelated normal distributed
random numbers X_1, X_2 ,

Define a new sequence Y_1 = rho X_1 + sqrt{1-rho^2} X_2

This new Y_1 ,
sequence will have a correlation of rho,
with the X_1 ,
sequence.
[edit]
Before and after correlating
Tags: ,
9Dec/090

MATLAB多线程实时工具——Timer对象

MATLAB先天的多线程残疾,给编写实时程序带来了致命的打击。一般情况只能按堆栈的形式,排队执行Matlab命令,可是无巧不成书,有时我们偏偏需要进行多线程和实时编程
比如:

1.我想实时更新一个电子表的时间,也就是每1s读取一次系统时间或者自动加1s

2.我需要这个程序在10s后终止

明显这些需要,使用一般的Matlab编程是没法或者说很难实现的。

还好Matlab自知在这方面亏欠大家太多,提供了唯一一个多线程的实现方法——Timer对象
一、Timer对象的创建
T = timer

T = timer(‘PropertyName1′, PropertyValue1, ‘PropertyName2′,
PropertyValue2,…)

括号中的参数可以在创建Timer是给出,也可以通过以后的set/get命令设置和读取
二、Timer属性参数

1.属性设置

(1)创建的时候直接赋值

%创建一个Timer对象变量TaskTimer,并给出部分参数

TaskTimer=timer(…

‘Name’,'Matlabsky’,…

‘TimerFcn’,@ExecuteTask,…

‘ErrorFcn’,@ExecuteError,…

‘Period’,1,…

‘ExecutionMode’,'fixedrate’);

(2)使用set/get设置和读取

%创建一个默认的Timer对象

TaskTimer=timer

set(TaskTimer,’Name’,'Matlabsky’,'TimerFcn’,@ExecutTask)

get(TaskTimer,’Name’)

其实set/get的使用方法不只这个,还其他形式,特别是在同时设置多个句柄的多个属性时,特别方便,感兴趣的网友可以直接参考Matlab的帮助文档

(3)使用结构体设置

TaskTimer=timer

TaskTimer.Name=’Matlabsky’

TaskTimer.TimerFcn=@ExecutTask
2.常用属性(下面是几个比较重要的,希望大家掌握)

(1)ExecutionMode   执行的模式

就是确定TimerFcn执行周期计算的始末点

’singleShot’
:只能执行一次,故Period属性不起作用,其他模式都可以执行多次

‘fixedSpacing’
:上一次TimerFcn执行完毕时刻到下一次TimerFcn被加入队列时刻之间的间隔

‘fixedDelay’
:上一次开始执行到下一次被加入队列之间的间隔

‘fixedRate’     :前后两次被加入到执行语句队列时刻之间的间隔

(2)Period  时间间隔(执行周期)

也就是每个Period执行一次TimerFcn,执行周期到底指的是哪两个时刻之间的间隔由ExecutionMode决定

(3)StartDelay   启动时延

从启动Timer开始到第一次把TimerFcn的加入到Matlab的执行语句队列中去的时延,
默认值为0s

(4)TasksToExecute    执行次数

TimerFcn被执行的次数,默认为1次,设置多次时需要设置执行周期Period

(5)TimerFcn    Timer函数

是Timer对象的核心,Timer的多线就是通过执行这个TimerFcn来实现的,也是Timer对象的必选参数

当然还有其它属性,但是相对来说比较次要,更详细的说明,感兴趣网友可以自己查看MATLAB帮助文档
三、回调函数callback
1.callback的建立

Timer对象的回调函数有TimerFcn、ErrorFcn、StartFcn和StopFcn,其中TimerFcn是必选的

Timer的所有回调函数callback,至少有两个输入参数obj和event,其经典格式为

  1. function my_callback_fcn(obj,event,p1, p2)
  2. %by dynamic
  3. %see also http://www.matlabsky.com
  4. %20092.15
  5. %
  6. %obj和event为必选输入参数
  7. %
  8. %obj就是前面创建的Timer对象,它包含着Timer的所有参数
  9. %
  10. %event是一个结构体,包含Type和Data两个字段,Type保存的是当前的Timer对象执行的时间,比如StartFcn、StopFcn等
  11. %
  12. %而Data还是一个结构体,包含time一个字段,保存的是执行该事件的系统时间
  13. %
  14. %所以event的使用是event.Type获取Timer到底在执行哪个事件,event.Data.time或者执行该时间的系统时间
  15. %
  16. %p1,p2…等是其它自定义的输入参数,根据需要在自己添加
复制代码
2.回调函数的属性值
  1. 回调函数
    属性值
  2. function
    myfcn(obj,event)
    set(T,’StartFcn’,'myfcn’)或者T.StartFcn=’myfcn’(下同)
  3. function
    myfcn(obj,event)
    set(T,’StartFcn’,@myfcn)(同上)
  4. function myfcn(obj,event,p1,p2)
    set(T,’StartFcn’,{‘myfcn’,p1,p2})或者T.StartFcn={‘myfcn’,p1,p2}(下同)
  5. function myfcn(obj,event,p1,p2)
    set(T,’StartFcn’,{@myfcn,p1,p2})(同上)
复制代码
3.回调函数实例
  1. function TimerFcn_callback(obj,event,var)
  2. txt=’ event occurred at ‘;
  3. event_type = event.Type;%获取当前事件
  4. event_time = datestr(event.Data.time);
  5. msg = [event_type txt event_time];
  6. disp(‘Hello it’s Matlabsky Timer Object Demo’)
  7. disp(‘Now! Begin…’)
  8. disp(msg)
  9. x=2*pi*[-1:0.01:1];
  10. y=var*sin(var*x)+var
  11. plot(x,y)
复制代码
四、Timer对象应用实例
  1. function matlabsky_timer_demo
  2. %
  3. %by dynamic
  4. %see also http://www.matlabsky.com
  5. %2009.1.23
  6. %
  7. clc
  8. %创建Timer对象
  9. T = timer(‘TimerFcn’,@TimerFcn,
    ‘StartDelay’,10,’StartFcn’,@StartFcn);
  10. %启动Timer
  11. start(T);
  12. %开始循环
  13. while T.userdata==’f'
  14. disp(‘Love Matlab Love Matlabsky——Timer
    Object Demo’);
  15. pause(1)
  16. end
  17. %终止Timer
  18. stop(T);
  19. function TimerFcn(obj, event)
  20. %10s后终止程序
  21. obj.userdata=’t';
  22. disp(['the BeiJing time is ' datestr(event.Data.time)])
  23. disp(’10s Has Escaped, Following The Orde of TimerFcn, Now I
    Have To Stop This Progam!’)
  24. function StartFcn(obj,event)
  25. obj.userdata=’f';
  26. disp(['the BeiJing time is ' datestr(event.Data.time)])
  27. disp(‘The Timer Object Demo Begin Now!’)
复制代码
Tags: , ,