今天是六一儿童节，昨天一天都在折腾体育缓考，简单的把流程记录下，以飨后人。
注：学校为北京航空航天大学，院系为软件学院，其他院系亦可参照
1、到指定医院开具“诊断证明书”，重要的就是那个公章
注：校医院只适合证明，不适合诊断，你懂的
2、到系教务处（6号楼312）找本科教务老师填写“缓考申请表”，填写完之后可能需要等待半天，因为还要院系领导签字
注：院系领导签字就是盖个私章，而且章就在教务老师那，不知为什么一定非要等半天，可能这就是中国特色吧
3、到校教务处（办公楼东配楼314）找教务处领导签字，至此表格完整
4、去打印店将表格复印两份，诊断证明书需复印在背面
注：最近的复印店为5号学生公寓地下室的打印店
5、原件交至校教务处、复印件交至系教务处和体育部学分管理处（体育场看台背面的大门进去2楼右拐左手第一间），至此全部搞定
P.S.前后要跑很多路，如果行动不方便的，最好找个好朋友帮忙跑下腿
今天下午去了学校的“心理咨询与发展辅导研究中心”。虽然很早觉得有必要去，但因为各种原因一直拖到今天。老师很和蔼，帮我整理思绪，聊得很愉快，心情平静了许多，也想明白了很多事情。
感谢心理老师马喜亭教授，以前选修过你的《大学生心理学》，当时就被马老师的魅力吸引了，虽然还是逃了很多次课，但那些却变成了遗憾，以后有机会一定去旁听。
感谢所有关心、帮助过我的朋友，你们为我做过的一切我都不会忘记，大家都要好好努力，未来在我们手里！
感谢爸爸妈妈，请原谅儿子的任性，一直让你们很担心。请原谅他的年轻倔强，他也在成长，他更在努力成为你们的骄傲！
感谢你，我的小狗，陪我走过那段美好的回忆，我会努力的向前，至少证明你当初没有看错。你也要好好加油，我会等你的，相信你一定可以……
虽然可能还是会痛苦，虽然可能还是会纠结，但我相信，他们的存在是有意义的，这就是生活的体验吧，这不正是我一直以来追求的吗？是的……
好久没发日志了，Calendar排的满满的，忙碌的学期就这样开始了……
1、ICST：实习机会很宝贵，想好好珍惜的，难得老师赏识，我也想踏踏实实真正做点东西
2、必修课：虽然实验室导师给予了保研的承诺，但还是以我有保研资格为前提的，这学期好多理论课，要努力静下心好好学学
3、Google Camp：说实话，现在对社团越来越迷惘了，或许我总是避免不了理想主义的禁锢，到头来总让自己不是滋味，决定逐渐放下了。小白，你要加油！
4、奶奶：苦了一辈子的老人家，到头来还要经受病痛的折磨。异地求学让陪伴都变成了奢望，然而回去又会陷入无能为力的哀伤，真不知道该怎样才好……
5、爸爸妈妈：如果世界上没有钱这种东西，一定会少了很多烦恼。你们辛苦了，你们要好好的，钱根本没有你们重要……
6、小狗：笨笨知道，你一个人更不容易。笨笨真的想可以和你能有结果，但是无论怎样，你一定要幸福……不然我不会饶了你！
夜深了，最近其实效率挺低的，还是有点心浮气躁吧，回去休息了，新的一天已经开始了……
今天，敏感日，不提，只记录一下外场申请的经历，个中感受自己参悟：
时间跨度大致为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：声源定位和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;
}