在返回上海的途中，我们顺便游玩了一下新安江，号称江南小三峡，水质不行，混浊不堪，水流不急，真不知道为什么被称为小三峡的，也许是古时候，这里的水流较急吧，不然当年伍子胥也不需要那个老渔夫送他过江了，今天也就没有“胥江野渡”这个景点了。两岸都是郁郁葱葱的树木，这天天色不太好，没有白天蓝云，否则的话，景色会更好看点。

   上岸后，爬了个小山丘，看看那个所谓的“葫芦飞瀑”，这些所谓的景点都是名字取得玄乎又玄，引人入胜，要真到那里了，也就那么回事，没有啥好看的。叫毛毛拍了几张照片做留念，证明老子曾经到此一游。

    回上海的车上，还是和过来时候一样——打牌消遣时光，不过这次打的是固定队友的，我和戴靓仔对家，不过把人家高人给拖累坏了，我们被裴氏传人打得体无完肤，惶惶若丧假之犬，我们想投降都no way，他们一定要慢慢地把我们折磨到死，简直就是没人性啊，现在我终于明白为什么美国老说我国没有人权了。

   嘉兴小便期间，杨老师给大家买了些粽子，每人五个，很赞啊！这些粽子放在lab402的冰箱里面，是后来几天我的夜宵。

   回到学校后，大家去国定365搓了一顿再回宿舍休息。

（四）建德激战

    周一早上六点二十分我就醒了，其实由于紧张过度，基本是彻夜未眠。起床后，洗个澡，然后再整理整理就下去吃早餐，八点左右开会去了。

    Coffee break之前我是听别人讲，不过也是身在曹营心在汉了，哪能听得进去啊，都在想着自己的讲稿。本来还以为report时可以拿个讲稿看看呢，后来发现情况不允许，因为一个手得拿话筒，一个手得拿激光笔。Coffee break之后，第三个是我，没办法，硬着头皮上了，死就死，怕什么啊。说实话，开讲之前很紧张的，”good morning……”之后，就一点都不怕了，感觉就跟下面没有人似的，我就是不停的blablabla…….,主持人提醒我还有一分钟时候，我刚好结束，时间还挺凑巧的，不过就是英语太烂，说起来不连贯，发音带广西口音。可能是下面的人也没有听懂，所以也没有什么问题，就那个最爱提问题的韩国人问了我个问题，后来一个问题也被这个家伙帮我回答了，Over了！解放了！之后就是认认真真的听别人讲了。

    说真的，参加这种会议还是挺有好处的，可以了解到当前光学的研究领域，和大家的研究兴趣所在，开阔了视野，不仅仅局限于自己实验室那几个课题。

   下午有个poster。做poster的同学终于出场了，我也就转了一圈而已，没怎么认真看。吃点点心，喝喝coffee，跟导游小姐瞎侃了一下。

（三）征战途中

    之前大多数时间都是在净化楼做实验，很少呆在科学楼，因此很多消息都不知道，就如同驻外大使，远离中央，既无法参与中央的决定，也无法及时获得中央的决议，导致出发开会前的三天我才知道我要在会议上作report，时间紧，任务急，只能没日没夜加班加点地prepare了，虽然知道自己很烂，但也应该尽一切可能地准备好一点了，好歹也是个国际会议，总不能丢国家的脸吧。

    一般情况下，你没事的话，就会一直都没什么事情的，当你忙的时候，事情会一件又一件的给你冒出来的。本来就为这个report搞死了，ZXGJ的笔试面试恰好在这个时候过来，在本来就不够的三天时间中还得抽出一天去张江笔试和面试，这不是火上加油，雪中送冰吗！面试是没法准备的了，只好匆忙应战了，幸好有逃之夭夭姐姐帮忙。

    周五早上6点起床，七点坐校车到张江校区，然后打的到ZXGJ，逃之夭夭姐姐请我吃公司的早餐，不过一心想着笔试面试的事情，没什么胃口，吃了两口就吃不下了。九点半hr的小姐先给我们几个人介绍一下ZXGJ的概况，十点开始笔试，内容就是英语和逻辑题，英语是四段短文阅读，和CET6差不多，逻辑题嘛，本来就是我的强项，一点也不怕，不过太紧张了，我最拿手的那两道数列题居然没做出来，SHIT!,十一点开始是性格测试题。

    考完试后，逃之夭夭姐姐请我去公司的饭堂吃饭，感觉还挺不错，和学校的差不多。

    下午就是面试了，先是hr部门的chen小姐面试我，hr嘛，半桶水一堆，就是问些漫无边际抽象无比的SB问题，这瞎扯蛋也是我的强项，对于这种问题，答案就是瞎扯蛋，问到我最想感谢的三个人时候，我用最煽情最感人的方式回答了她我要感谢的人以及原因，当我回答完之后，发现她眼泪都快要流出来了，哎！太投入了。至于其他问题嘛，信手捏来，从容不迫，侃侃而谈，谈笑风生，时光在不知不觉中流逝，刹那间已面了三十多分钟，面我的时间比其它几个人的加起来都长，唉！我这人就有这毛病，侃起来就没法stop，这下侃多了。之后是部门主管面试，有三个主管，其中两个chinese，一个korean。对于Chinese，我当然不怕了，这辈子咱们怕过谁啊，不过面试那个高丽人还真是有点紧张，因为要用到英语，这English啊，我是烂到没法再烂了。第一个面试的是chinese，我靠！碰上狗屎运了我，居然是老乡，老乡嘛，当然是聊聊家乡的情况了，再跟他套套公司的前景和前途如何，聊得不亦乐乎！开心啊，有个老乡助阵，成功的概率突然从0提高到了50%。Next one, 高丽人，据逃之夭夭姐姐说这个人还挺好的，不过我看他横看竖看正看侧看都像个十足的大色狼；忘了带个英文版简历过去，那个中文版的他又看不懂，所以剩下的就是还是瞎扯蛋，他夸我英文好，我夸他中文好，娘的！这家伙居然问我专业问题，没办法了，只能使出我的杀手锏了“I understand the questions you asked, but I can not answer you in English, because my English is so bad! I can answer you in Chinese if you can understand Chinese…..”，他就无语了，转而聊其它非技术类的问题了，比如爱好等。OK! The last one, a Chinese guy，人挺好的，头有点秃，一看就是高技术类，他就问了些比较实际的问题，比如，学过什么半导体相关的知识，对这个部门有兴趣吗，以前在学生会做部长时主要是做什么东西的啊，公司的薪水这么低你为什么还愿意过来啊……。面试结束时，恰好逃之夭夭姐姐过来，跟他们聊了一会儿，说这是我是她师弟，希望大家多多关照之类的；然后带我到她那里坐了一会儿聊了聊。之后我就匆匆赶回学校了，因为还得准备中韩会议的report呢。

   这段时间确实太忙了，忙到我根本找不出个合适的词来形容了，不过总算忙过去了。回首这段时间真是不胜唏嘘，从中韩光学会议前的没日没夜的准备，到大会上的演讲，再到新安江上赏景，最后回到学校修改文章……

两点到了北馆，远远就看见那个师傅跟我笑眯眯地打招呼了
“你是那个***同学吧”
“是啊”，心一愣，妈的，他怎么认识我的
“你定的是两点到三点的主场？”
“yes!”
“不好意思啊，跟你说啊，那个场地现在有学校领导来玩，你去玩玩乒乓好吗，免费任你玩......”
......
不好也只能好了，有什么办法呢，谁叫人家是领导
顶你个肺! 老子大半年不来玩一次羽毛球了，怎么就那么倒霉呢，遇上这些
所谓的领导，想起了唐伯虎点秋香里面的："蛮横不留情，占我大屋夺我田....."

“那我定的那个三点到四点的边场呢？能玩吗”
“你还定了个边场啊？可以，可以...”
“那我买票了”
“这个不用了，免费给你玩，还真不好意思，把你的主场让给那些领导了，谢谢你的体谅啊...”
......

就这样拣了个场地玩玩，省了10块钱
.....
.....

国庆节这几天哪里都没去，
就整天呆在净化楼大炉子那里研究停判算法
本来雄心壮志的——“匈奴未破，何以为家”，“不破楼兰终不还”
几天下来，机关算尽，天不遂人愿
折戟沉沙，翦羽而归
被打击坏了我，能想得到的算法都用过了，还是不行
现在已是心疲力竭，无可奈何了
要么rp用完了，要么就是江郎才尽了
认命了！
妈的！怪不得当年就没有人想来搞这个项目

我的青春和心血就这样无私的奉献给了这个实验室
现在连找工作投简历的时间都没有
搞不好到最后，毕业就失业
想想就觉得不寒而栗，毛骨悚然

等这项目验收了，打死都不干了
这冤大头不能再做了
也许也不能说是冤大头
应该说是这种不是自己能胜任的活绝对不能干了
管它学位不学位的
能找到好工作才是王道啊
找不到好工作，给个博士学位都没意思

===============================================
小整理：MATLAB常用的基本数学函数
abs(x)：纯量的绝对值或向量的长度
angle(z)：复数z的相角(Phase angle)
sqrt(x)：开平方
real(z)：复数z的实部
imag(z)：复数z的虚部
conj(z)：复数z的共轭复数
round(x)：四舍五入至最近整数
fix(x)：无论正负，舍去小数至最近整数
floor(x)：地板函数，即舍去正小数至最近整数
ceil(x)：天花板函数，即加入正小数至最近整数
rat(x)：将实数x化为分数表示
rats(x)：将实数x化为多项分数展开
sign(x)：符号函数 (Signum function)。

当x<0时，sign(x)=-1；

当x=0时，sign(x)=0;

当x>0时，sign(x)=1。
rem(x,y)：求x除以y的馀数
gcd(x,y)：整数x和y的最大公因数
lcm(x,y)：整数x和y的最小公倍数
exp(x)：自然指数
pow2(x)：2的指数
log(x)：以e为底的对数，即自然对数或
log2(x)：以2为底的对数
log10(x)：以10为底的对数
===============================================
小整理：MATLAB常用的三角函数
sin(x)：正弦函数
cos(x)：馀弦函数
tan(x)：正切函数
asin(x)：反正弦函数
acos(x)：反馀弦函数
atan(x)：反正切函数
atan2(x,y)：四象限的反正切函数
sinh(x)：超越正弦函数
cosh(x)：超越馀弦函数
tanh(x)：超越正切函数
asinh(x)：反超越正弦函数
acosh(x)：反超越馀弦函数
atanh(x)：反超越正切函数
====================================================
小整理：适用於向量的常用函数有：
min(x): 向量x的元素的最小值
max(x): 向量x的元素的最大值
mean(x): 向量x的元素的平均值
median(x): 向量x的元素的中位数
std(x): 向量x的元素的标准差
diff(x): 向量x的相邻元素的差
sort(x): 对向量x的元素进行排序（Sorting）
length(x): 向量x的元素个数
norm(x): 向量x的欧氏（Euclidean）长度
sum(x): 向量x的元素总和
prod(x): 向量x的元素总乘积
cumsum(x): 向量x的累计元素总和
cumprod(x): 向量x的累计元素总乘积
dot(x, y): 向量x和y的内积
cross(x, y): 向量x和y的外积
（大部份的向量函数也可适用於矩阵，详见下述。）
====================================================
小整理：MATLAB的查询命令
help：用来查询已知命令的用法。例如已知inv是用来计算反矩阵，键入
help inv即可得知有关inv命令的用法。（键入help help则显示help的用
法，请试看看！）
lookfor：用来寻找未知的命令。例如要寻找计算反矩阵的命令，可键入
lookfor
inverse，MATLAB即会列出所有和关键字inverse相关的指令。找到所需的
命令後，即可用help进一步找出其用法。（lookfor事实上是对所有在搜寻
路径下的M档案进行关键字对第一注解行的比对，详见後叙
===============================================
下表即为MATLAB常用到的永久常数。
小整理：MATLAB的永久常数
i或j：基本虚数单位（即）
eps：系统的浮点（Floating-point）精确度
inf：无限大， 例如1/0
nan或NaN：非数值（Not a number），例如0/0
pi：圆周率 p（= 3.1415926...）
realmax：系统所能表示的最大数值
realmin：系统所能表示的最小数值
nargin: 函数的输入引数个数
nargin: 函数的输出引数个数
====================================================
1-7、结束MATLAB

有三种方法可以结束MATLAB：

   1.键入exit
   2.键入quit
   3.直接关闭MATLAB的命令视窗（Command window）
====================================================
小整理：MATLAB基本绘图函数
plot: x轴和y轴均为线性刻度（Linear scale）
loglog: x轴和y轴均为对数刻度（Logarithmic scale）
semilogx: x轴为对数刻度，y轴为线性刻度
semilogy: x轴为线性刻度，y轴为对数刻度
====================================================
若要画出多条曲线，只需将座标对依次放入plot函数即可：
plot(x, sin(x), x, cos(x));
若要改变颜色，在座标对後面加上相关字串即可：
plot(x, sin(x), 'c', x, cos(x), 'g');
若要同时改变颜色及图线型态（Line style），也是在座标对後面加上相
关字串即可：
plot(x, sin(x), 'co', x, cos(x), 'g*');
====================================================
小整理：plot绘图函数的叁数
字元  颜色 字元  图线型态
y  黄色 .  点
k  黑色 o  圆
w  白色 x  x
b  蓝色 +  +
g  绿色 *  *
r  红色 -  实线
c  亮青色 :  点线
m  锰紫色 -.  点虚线
   --  虚线
====================================================
图形完成後，我们可用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范
围：

axis([0, 6, -1.2, 1.2]);

此外，MATLAB也可对图形加上各种注解与处理：

xlabel('Input Value'); % x轴注解
ylabel('Function Value'); % y轴注解
title('Two Trigonometric Functions'); % 图形标题
legend('y = sin(x)','y = cos(x)'); % 图形注解
grid on; % 显示格线

我们可用subplot来同时画出数个小图形於同一个视窗之中：

subplot(2,2,1); plot(x, sin(x));
subplot(2,2,2); plot(x, cos(x));
subplot(2,2,3); plot(x, sinh(x));
subplot(2,2,4); plot(x, cosh(x));

MATLAB还有其他各种二维绘图函数，以适合不同的应用，详见下表。

====================================================
小整理：其他各种二维绘图函数
bar  长条图
errorbar  图形加上误差范围
fplot  较精确的函数图形
polar  极座标图
hist  累计图
rose  极座标累计图
stairs  阶梯图
stem  针状图
fill  实心图
feather  羽毛图
compass  罗盘图
quiver  向量场图
====================================================
以下我们针对每个函数举例。

当资料点数量不多时，长条图是很适合的表示方式：

close all; % 关闭所有的图形视窗
x=1:10;
y=rand(size(x));
bar(x,y);

如果已知资料的误差量，就可用errorbar来表示。下例以单位标准差来做
资料的误差量：

x = linspace(0,2*pi,30);
y = sin(x);
e = std(y)*ones(size(x));
errorbar(x,y,e)

对於变化剧烈的函数，可用fplot来进行较精确的绘图，会对剧烈变化处进
行较密集的取样，如下例：

fplot('sin(1/x)', [0.02 0.2]); % [0.02 0.2]是绘图范围

若要产生极座标图形，可用polar：

theta=linspace(0, 2*pi);
r=cos(4*theta);
polar(theta, r);

对於大量的资料，我们可用hist来显示资料的分　情况和统计特性。下面
几个命令可用来验证randn产生的高斯乱数分　：

x=randn(5000, 1); % 产生5000个 ?=0，?=1 的高斯乱数
hist(x,20); % 20代表长条的个数

rose和hist很接近，只不过是将资料大小视为角度，资料个数视为距离，?⒂眉座标绘制表示?
x=randn(1000, 1);
rose(x);

stairs可画出阶梯图：

x=linspace(0,10,50);
y=sin(x).*exp(-x/3);
stairs(x,y);

stems可产生针状图，常被用来绘制数位讯号：

x=linspace(0,10,50);
y=sin(x).*exp(-x/3);
stem(x,y);

stairs将资料点视为多边行顶点，并将此多边行涂上颜色：

x=linspace(0,10,50);
y=sin(x).*exp(-x/3);
fill(x,y,'b'); % 'b'为蓝色

feather将每一个资料点视复数，并以箭号画出：

theta=linspace(0, 2*pi, 20);
z = cos(theta)+i*sin(theta);
feather(z);

compass和feather很接近，只是每个箭号的起点都在圆点：

theta=linspace(0, 2*pi, 20);
z = cos(theta)+i*sin(theta);
compass(z);
--
3.基本XYZ立体绘图命令

在科学目视表示（Scientific visualization）中，三度空间的立体图是
一个非常重要的技巧。本章将介绍MATLAB基本XYZ三度空间的各项绘图命
令。

mesh和plot是三度空间立体绘图的基本命令，mesh可画出立体网状图，
plot则可画出立体曲面图，两者产生的图形都会依高度而有不同颜色。下
列命令可画出由函数 形成的立体网状图:
x=linspace(-2, 2, 25); % 在x轴上取25点
y=linspace(-2, 2, 25); % 在y轴上取25点
[xx,yy]=meshgrid(x, y); % xx和yy都是21x21的矩阵
zz=xx.*exp(-xx.^2-yy.^2); % 计算函数值，zz也是21x21的矩阵
mesh(xx, yy, zz); % 画出立体网状图

surf和mesh的用法类似：

x=linspace(-2, 2, 25); % 在x轴上取25点
y=linspace(-2, 2, 25); % 在y轴上取25点
[xx,yy]=meshgrid(x, y); % xx和yy都是21x21的矩阵
zz=xx.*exp(-xx.^2-yy.^2); % 计算函数值，zz也是21x21的矩阵
surf(xx, yy, zz); % 画出立体曲面图

为了方便测试立体绘图，MATLAB提供了一个peaks函数，可产生一个凹凸有
致的曲面，包含了三个局部极大点及三个局部极小点，其方程式为：

要画出此函数的最快方法即是直接键入peaks：

peaks

z = 3*(1-x).^2.*exp(-(x.^2) - (y+1).^2) ...
- 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2) ...
- 1/3*exp(-(x+1).^2 - y.^2)

我们亦可对peaks函数取点，再以各种不同方法进行绘图。meshz可将曲面
加上围裙：

[x,y,z]=peaks;
meshz(x,y,z);
axis([-inf inf -inf inf -inf inf]);

waterfall可在x方向或y方向产生水流效果：

[x,y,z]=peaks;
waterfall(x,y,z);
axis([-inf inf -inf inf -inf inf]);

下列命令产生在y方向的水流效果：

[x,y,z]=peaks;
waterfall(x',y',z');
axis([-inf inf -inf inf -inf inf]);

meshc同时画出网状图与等高线：

[x,y,z]=peaks;
meshc(x,y,z);
axis([-inf inf -inf inf -inf inf]);

surfc同时画出曲面图与等高线：

[x,y,z]=peaks;
surfc(x,y,z);
axis([-inf inf -inf inf -inf inf]);

contour3画出曲面在三度空间中的等高线：

contour3(peaks, 20);
axis([-inf inf -inf inf -inf inf]);

contour画出曲面等高线在XY平面的投影：

contour(peaks, 20);

plot3可画出三度空间中的曲线：

t=linspace(0,20*pi, 501);
plot3(t.*sin(t), t.*cos(t), t);

亦可同时画出两条三度空间中的曲线：

t=linspace(0, 10*pi, 501);
plot3(t.*sin(t), t.*cos(t), t, t.*sin(t), t.*cos(t), -t);

y(2:4)-1 % 取出y的第二至第四个元素来做运算

同样的方法可用於产生公差为1的等差数列：x = 7:16

x = 7:3:16 % 公差为3的等差数列

x = linspace(4, 10, 6) % 等差数列：首项为4,末项为10,项数为6

若要重新安排矩阵的形状，可用reshape命令：B = reshape(A, 4, 2) % 4是新矩阵的列数，2是新矩阵的行数

举例来说，下列命令会产生一个长度为6的调和数列（Harmonic
sequence）：
x = zeros(1,6); % x是一个16的零矩阵
for i = 1:6,
x(i) = 1/i;
end

for咳可以是多层的，下例产生一个16的Hilbert矩阵h，其中为於第i
列、第j行的元素为：
h = zeros(6);
for i = 1:6,
for j = 1:6,
h(i,j) = 1/(i+j-1);
end
end

format rat % 使用分数来表示数值
>>disp(x)
 1 1/2 1/3 1/4 1/5 1/6

function output = fact(n)
% FACT Calculate factorial of a given positive integer.
output = 1;
for i = 1:n,
output = output*i;
end
其中fact是函数名，n是输入引数，output是输出引数，而i则是此函数用
到的暂时变数。要使用此函数，直接键入函数名及适当输入引数值即可：

MATLAB的函数也可以是递渴降模≧ecursive），也就是说，一个函数可以
呼叫它本身。举例来说，n! =
n*(n-1)!，因此前面的阶乘函数可以改成递渴降男捶ǎ
function output = fact(n)% FACT Calculate factorial of a given positive integer recursively.
if n == 1, % Terminating condition
output = 1;
return;
end
output = n*fact(n-1);
在写一个递亢数时，一定要包含结束条件（Terminating
condition），否则此函数将会一再呼叫自己，永远不会停止，直到电脑的
记忆体被耗尽为止。以上例而言，n==1即满足结束条件，此时我们直接将
output设为1，而不再呼叫此函数本身。