怨念的外场申请
今天，敏感日，不提，只记录一下外场申请的经历，个中感受自己参悟：
时间跨度大致为9：00-16：00，中间包含午饭、午休时间
0、分团委盖章是前一天就搞定的，导人没话说
1、9：00，保卫处，教学区门口，负责人（杨队长）说外场一律不批，他只听（副）处长的话
2、保卫处副处长办公室，西办公楼111，金副处长，口有点吃，说校领导不给批，一副“关我屁事”的架势，当被问及可否找校领导时，吞吐迟疑半天，说去找校团委副书记
3、校团委副书记办公室，知行北楼219，无人
4、校团委副书记办公室旁边的办公室，一位很和气的女老师，未问及姓名，道出真相，是保卫处不批他们的申请，他们根本不做主，但透露可以找党政办“咨询”相关情况
5、校党政办办公室1，办公楼中厅216，一位男老师，说应该找团委、保卫处，当被告之我们的经历后，去了旁边的办公室好长时间才出来，让我们进去找另一位老师
6、校党政办办公室2，一位很好的女老师，打电话联系了分团委、团委、保卫处等若干负责人，疏通了所有关系，但是提出申请表需要修改，主要因为在校团委意见一栏，盖了分团委的章
7、保卫处副处长办公室，想给那个踢皮球的副处长杀个回马枪，结果等了可能将近有一个小时，人都没有出现
8、11：00，去吃午饭，下午2点办公，期间按照党政办老师的要求重新打印了申请表
9、13：50，导办，分团委的章，导可能因为事情一直弄到了党政办，又具体的问了很多情况，不过还是很爽快的盖了章
10、保卫处副处长办公室，似乎在忙什么，一副不理不睬的样子，当提及党政办已经联系过他们处长时，又说处长没有和他说，处长在开会，然后就让我坐在旁边等
11、保卫处副处长办公室对门的办公室，处长开会回来，今副处长马上又把我当皮球踢给了刚回来坐在对面办公室的处长，陈处长了解了活动大体情况之后，说了很多官话，显得很正直，例如：“学生活动我们是一定会支持的！”。先说要团委书记签字，当被告之是分团委时，又说必须还有分团委书记签字，然后又说必须有院系盖章，至少要是副书记签字，我说分团委的章加上副书记签字可不可以，又被告之可以。
12、软件学院副书记办公室，宋老师在忙，不过很热情的接待了我，了解了一下活动情况，就答应了，说找分团委辅导员，那里有他的签字章
13、导依旧没话说，很爽快
14、保卫处处长办公室，西办公楼119，无人
15、保卫处副处长办公室，保卫科（见1）的队长还有另外几个人在里面闲聊，等了一会今副处长才来，貌似电脑出了一点小问题，于是乘势对我保持漠视态度，期间眼神呆滞的看着近乎死机的电脑或者其他聊天的几个人若干分钟，直到聊天中的一个人看着我站了半天，问我有什么事的时候，在旁人催促下，他终于开始处理我的事情。问了一堆全部被我挡回去的话，然后起身去找处长未果。回来又对了不超过100字的几张申请表深情的望了若干分钟，在旁边一个人的催促下终于签上“同意”二字，还附带了一堆注释——禁止商业性，不得使用音响设备，但*单位有专人管理（*为不认识的一个字），又去了另一个办公室盖了保卫处的章。起初连存根都没留就给我了，在旁边人的提醒下，才在我出门后叫住我，要了份作为存根。
16、至此算是搞定了，只不过花费了前后差不多7个小时的时间而已。
感谢：刘导、校团委不知名女老师、党政办不知名女老师、软件学院宋副书记（以事情发展先后顺序）和所有指路的保安
保留对保卫处相关人员评价的权利，尤其是金副处长，放心，我不会利用你口吃的弱点骂你的，因为鄙人觉得你还没有被我骂的资格
P.S.另外，水杯在漫漫外场申请路中丢失，寻找未果，留此祭奠
北航物理实验E04——Code
物理实验E04：声源定位和GPS模拟，纠结了好久，竟然还要交数据处理的源程序
程序倒不难，关键是不知道数据怎么处理的话就……于是我就杯具了……
网上竟然没有什么现成的程序，看样大家还不够团结啊～
贴在这里，完全开源，请遵守GPL协议，造福后人……
/*
*ID:38211214
*Program:声源定位
*/
/*
*in.txt
* 45.4 73.7 0 42.7
* 0 55.0 107.8 126.3
*100.2 112.0 0 25.8
*106.0 78.7 53.0 0
* 0 34.6 9.3 34.6
* 0 25.6 74.0 85.5
* 81.5 52.6 50.1 0
* 0 62.3 134.8 151.6
* 1.3 0 94.3 91.6
*/
/*
*out.txt
*( 53.29,153.49)
*( 60.26,393.04)
*( 98.92, 49.82)
*(254.10,103.05)
*( 89.57,364.71)
*( 98.04,317.42)
*(249.78,142.83)
*( 50.46,403.15)
*(151.00,405.37)
*/
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
const double sensor_x[4] = { 0, 298, 0, 298}; //传感器横坐标
const double sensor_y[4] = { 446, 446, 0, 0}; //传感器纵坐标
const double c = 2970000; //常数，声音的传播速度，单位"mm/s"
const double pi = acos(-1.0); //常数π
double det_t[4] = {0}; //时间间隔
double A, B, D, phi; //计算公式中间量A, B, D, Φ
double radius, theta; //声源坐标，极坐标表示(r,θ)
double source_x, source_y; //声源坐标，直角坐标表示
void cal(int sign, double t1, double t2, double x1, double y1, double x2, double y2, double px, double py) //计算公式的实现
{
double det1 = c * t1, det2 = c * t2;
//中间量的计算
A = x2 * (x1 * x1 + y1 * y1 - det1 * det1) - x1 * (x2 * x2 + y2 * y2 - det2 * det2);
B = y2 * (x1 * x1 + y1 * y1 - det1 * det1) - y1 * (x2 * x2 + y2 * y2 - det2 * det2);
D = det1 * (x2 * x2 + y2 * y2 - det2 * det2) - det2 * (x1 * x1 + y1 * y1 - det1 * det1);
phi = atan(B / A);
theta = phi + sign * acos(D / sqrt(A * A + B * B));
radius = (x1 * x1 + y1 * y1 - det1 * det1 ) / (2 * ( x1 * cos(theta) + y1 * sin(theta) + det1 ));
//最终结果及输出
source_x = cos (theta) * radius;
source_y = sin (theta) * radius;
cout < <fixed <<setprecision(2)
<<'(' <<setw(6) <<source_x + px << ',' <<setw(6) <<source_y + py <<')' <<endl;
}
int main()
{
freopen("in.txt", "r", stdin); //打开输入文件
freopen("out.txt", "w", stdout); //打开输出文件
for (int i = 0; i < 9; ++ i)
{
for (int j = 0; j < 4; ++ j) //数据读入
{
cin >>det_t[j];
det_t[j] /= 1e6;
}
for (int j = 0; j < 4; ++ j)
{
if (fabs(det_t[j] - 0) < 1e-6) //判断最先到达的传感器编号
{
switch (j)
{
case 0:
cal(-1, det_t[1], det_t[3], 298, 0, 298, -446, 0, 446);
break;
case 1:
cal(-1, det_t[0], det_t[2], -298, 0, -298, -446, 298, 446);
break;
case 2:
cal(1, det_t[0], det_t[1], 0, 446, 298, 446, 0, 0);
break;
case 3:
cal(1, det_t[0], det_t[1], -298, 446, 0, 446, 298, 0);
break;
}
}
}
}
return 0;
}
/*
*ID:38211214
*Program:GPS仿真
*/
/*
*in.txt
* 87.7 113.7 88.1 85.0 224.0
* 89.9 67.7 124.9 203.0 303.0
*147.5 123.5 95.9 240.0 109.0
*112.2 91.3 106.0 208.0 216.0
*112.4 140.7 56.0 49.0 137.0
* 90.3 79.9 112.1 179.0 269.0
*128.8 88.7 121.0 272.0 209.0
* 67.2 62.5 142.3 158.0 369.0
* 94.3 111.3 98.2 128.0 228.0
* 78.5 79.8 118.5 149.0 297.0
*/
/*
*out.txt
*( -9.51,458.73)
*(313.00,460.32)
*( -8.63, -7.67)
*/
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
const double c = 2970; //常数：声音传播速度
double det_t[3][10] = {0}; //三个接收器的时间差
double trans_x[10],trans_y[10]; //发送器（transmitter）坐标;
double a1, a2, a3, a4, b1, b2, mult, D, E; //算法过程中间量
void gps(double t[], double x1, double y1, int num)
{
D = E = 0;
do
{
a1 = a2 = a3 = a4 = b1 = b2 = 0;
for (int i = 0; i < 10; ++ i)
{
double m1 = trans_x[i] * 1e-3 - x1;
double m2 = y1 - trans_y[i] * 1e-3;
double m3 = (c * t[i] * 1e-6) * (c * t[i] * 1e-6);
a1 = a1 - 3 * (trans_x[i] * 1e-3 - x1) * (trans_x[i] * 1e-3 - x1) + m2 * m2 - m3;
a2 = a2 + 2 * m1 * m2;
a3 = a2;
a4 = a4 - 3 * (y1 - trans_y[i] * 1e-3) * (y1 - trans_y[i] * 1e-3) + m1 * m1 - m3;
b1 = b1 + (m1 * m1 + m2 * m2 - m3) * (- m1);
b2 = b2 + (m1 * m1 + m2 * m2 - m3) * m2;
}
mult = a3 / a1;
a4 = a4 - a2 * mult;
b2 = b2 - b1 * mult;
E = b2 / a4;
D = (b1 - a2 * E) / a1;
x1 += D;
y1 += E;
}
while (D * D + E * E <= 1e-9); //达到精度之后停止
cout <<fixed <<setprecision(2)
<<'(' <<setw(6) <<x1 * 1e3 <<',' <<setw(6) <<y1 * 1e3 <<')' <<endl;
}
int main()
{
freopen("in.txt", "r", stdin); //打开输入文件
freopen("out.txt", "w", stdout); //打开输出文件
for (int i = 0; i < 10; ++ i) //读入数据
cin >>det_t[0][i] >>det_t[1][i] >>det_t[2][i] >>trans_x[i] >>trans_y[i];
//三次调用计算过程
gps(det_t[0], 0, 0.45, 1);
gps(det_t[1], 0.3, 0.45, 2);
gps(det_t[2], 0, 0, 3);
return 0;
}
Sleepless Night
好长时间没更新了，最近各种忙，连睡眠都成了奢望……
刷的第N个夜了，用音乐驱散倦意，终究还是耐不住寂寞，只好到这里来了……
最近，越来越觉得自己是一个Idealist，理性的外表下是深切的感性，面对世事的无常习惯于用我心看世界，坚持着自己内心深处那近乎不切实际的所谓理想却始终不愿放弃，这是狂妄还是一种超然？
最近，发现自己竟然爱上了音乐，硬盘里下满了贝多芬、莫扎特、柴可夫斯基的音乐，还有各种流行，贝多芬的《Piano Concertos Nos. 1-5》，莫扎特的《安魂曲》，苏打绿的《夏.狂热》《无与伦比的美丽》《小宇宙》《苏打绿》，这都是留下比较深刻印象的，还有各种流行，各种经典，突然有一种相见恨晚的感觉……
最近，莫名其妙的恋上了文字，看了《绿野仙踪》，看了《理想国》，看了《心理助人精要》，看了开复博士的《世界因你不同》《与未来同行》《做最好的自己》，看了《杜拉拉升职记》和《杜拉拉2华年似水》，看了《山楂树之恋》，看了痞子蔡的《鯨魚女孩·池塘男孩》，还看了张小娴的若干纯情文学《荷包里的单人床》《卖海豚的女孩》《我们都是丑小鸭》，收获颇丰～
最近，疯狂看了若干电影，《V for Vendetta》《(500) Days of Summer》《The Butterfly Effect》《Oktapodi》《Avatar》《Up》《The Hurt Locker》《It’s Tough to Be a Bird》《Body Heat》，一堆又一堆的感悟，感谢北航BT！(貌似北航校园网才可以访问)
最近，终于开始觉得自己正在走向理想的大学生活，以Coding为乐，以技术为伴，以GC为友，音乐、文字、电影成为了生活的娱乐，生活忙忙碌碌，却充实而满足，
最近，很开心，很幸福，因为小狗，因为善良的、单纯的、美丽的小狗，我明白异地、家庭给她的压力，笨笨会好好照顾自己，好好珍惜小狗……
Hello World Press！
历经艰辛，终于脱离BSP加入了Wordpress行列，值此新春佳节之际，不由得牢骚一番。
标题中的关键词来自一位朋友的误用，不过仔细琢磨一下，觉得不无道理，便斗胆在此误人子弟一回了……
觉醒……
突然睡不着了，开学一个星期了，却始终找不到开学的感觉……我和舍友调侃说：“这才是大学……”可是我明白，这不应该是我的大学……
从南京回来之后，心里一直空空的。知道自己该努力了，知道大学不是用来苦学如何荒废人生的地方。近期太多的经历让我突然明白了许多事情，在种种现实交织的生活面前，我变了许多，这是成长，可是，似乎还缺点什么……
我不缺乏动力，我学习着自己热爱的专业，我更拥有着所有爱我的关心我的人带给我的勇气。可是，我却似乎一直在独自徘徊，知道自己该前进，却连怎样迈出步伐都抛在了脑后……难道一定要失败之后才知道要努力吗？我不希望也不允许自己成为那样！机会是有限的，时间是宝贵的，我没有选择的权利！
今天（9.19）是自由软件日，北航也会举办自由软件日中国站相关活动，昨天北航开源软件俱乐部的部分成员讨论了会前筹备工作，我也混进去了。大家的热情都很高，一种前所未有的感觉。由于是开源软件的活动，没有任何经费预算，只是提供了技术支持（请牛人来演讲）和一些小礼品，大家忙前忙后，没有人为了什么，要说有，也只是为了心中的那份热情……这种感觉真的很好，我想这也是我需要的吧。
过得太现实了，便会忘记了真爱！
考虑的太多了，便会束缚了自己！
是该放下更多的东西了，对我前途负责的只有我自己！
不想再说什么了，下决心早已不是创意，做点自己想做的事情吧，这也许就足够了……
P.S.
以后我的博客里面会适当添加相关链接，感谢大家关注。
P.S.P.S.
饭否流民正式定居twitter！（这个链接需要翻墙才能访问）
认识的人都来吼声啊！
