怨念的外场申请
今天，敏感日，不提，只记录一下外场申请的经历，个中感受自己参悟：
时间跨度大致为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;
}
国庆，真的来了……
本来说好早起去学校的如心报告厅看国庆直播的，结果昨晚睡得太晚直接导致早上7点才起床。急急忙忙赶去了如心，不出所料前面的好座位都被占了，而且还发现了帮忙占座的。我就奇怪了，他们是来自习的吗？
不过虽然没坐在前面，我们还是很幸运的，因为最后报告厅里面坐满了人，最后其他人都是站着或者坐在地上的，我终于知道报告厅铺地毯的真正用途了……
话说中国最大的优势就是人口优势，今天也得到了充分的体现，来自全国各地的北航爱国学子齐聚在了这里，一同见证了伟大祖国60年华诞的历史。整个国庆汇演过程中，一阵阵的掌声在耳边回响，将报告厅的气氛一次次引入高潮。每当镜头对准温家宝总理的时候，会场内都会响起一阵掌声。每当那种情形浮现在眼前，心中都对温总理深深的人格魅力所折服，我想，这也许才是人生的意义之所在吧。在各省市直辖市的花车方阵经过时，会场中又是一阵一阵的掌声，虽然每一辆花车带来了不同地域的掌声，但是我们都明白，中华民族的心是永远紧紧相连的。当然，凡事总有例外，当四川、西藏还有台湾的花车经过时，大家不约而同的都鼓起了掌，不用说大家也都应该明白其中的意义。
一直到12点，汇演结束后，吃午饭，宅宿舍。玩了一会电脑，一觉睡到了6点，然后又在网上闲荡了一会，就看起了国庆晚会的现场直播。很精彩，很震撼，很自豪！
国庆真的到了，那些激动和感慨，并未因为午夜的来临而又丝毫的褪色。国庆八天假期，没什么太大的活动，外面人多，也不想到处乱转，正好可以散散心、看看书、调整调整状态，之前都太浮躁、太低效了，希望借此能有所改善。
最后祝所有朋友国庆快乐！华夏儿女，普天同庆！
开学一周……
开学一周……
选课结束……
繁重的课表……
悲惨的学期……
开始了……
P.S. 省略的课程还有
每两周一次的物理实验（每次要半天！）
每周五下午还有党支部组织的支教（入党积极分子必须参加……囧）
P.S.P.S.课表如下：
上面的两个不算才48课时……48课时！！！
（配图不太清楚，可以点击查看原图）
觉醒……
突然睡不着了，开学一个星期了，却始终找不到开学的感觉……我和舍友调侃说：“这才是大学……”可是我明白，这不应该是我的大学……
从南京回来之后，心里一直空空的。知道自己该努力了，知道大学不是用来苦学如何荒废人生的地方。近期太多的经历让我突然明白了许多事情，在种种现实交织的生活面前，我变了许多，这是成长，可是，似乎还缺点什么……
我不缺乏动力，我学习着自己热爱的专业，我更拥有着所有爱我的关心我的人带给我的勇气。可是，我却似乎一直在独自徘徊，知道自己该前进，却连怎样迈出步伐都抛在了脑后……难道一定要失败之后才知道要努力吗？我不希望也不允许自己成为那样！机会是有限的，时间是宝贵的，我没有选择的权利！
今天（9.19）是自由软件日，北航也会举办自由软件日中国站相关活动，昨天北航开源软件俱乐部的部分成员讨论了会前筹备工作，我也混进去了。大家的热情都很高，一种前所未有的感觉。由于是开源软件的活动，没有任何经费预算，只是提供了技术支持（请牛人来演讲）和一些小礼品，大家忙前忙后，没有人为了什么，要说有，也只是为了心中的那份热情……这种感觉真的很好，我想这也是我需要的吧。
过得太现实了，便会忘记了真爱！
考虑的太多了，便会束缚了自己！
是该放下更多的东西了，对我前途负责的只有我自己！
不想再说什么了，下决心早已不是创意，做点自己想做的事情吧，这也许就足够了……
P.S.
以后我的博客里面会适当添加相关链接，感谢大家关注。
P.S.P.S.
饭否流民正式定居twitter！（这个链接需要翻墙才能访问）
认识的人都来吼声啊！
