2017数学建模总结

迟到的建模总结，昨晚12点交完论文，打印材料，今天一天满课，刚刚把上周作业写完的我，终于有空闲时间写写这个总结了。

简单概括游戏体验非常差！！！

第一天晚上选题，看了看A，看了看B，B题分析性太强，属于应用统计，数据分析类，和队友商量着战略性放弃，然后就开始查资料，搞东西了，突然回忆起以前姜起源的书上貌似讲过CT成像这种东西,翻出来看果然有这玩意儿。然后看了看题目发现精度是4位小数，似乎可以暴力枚举答案啊(ACMer的特性)，再仔细对照姜老的模型，发现参数太多了，180个角度，还有中心，连单元距离也没告诉，，，，图片是256x256,几乎搜到建模结束也不可能搜出来，突然陷入绝望，本着第一晚看看题的想法，就这样回宿舍碎觉了。

第二天上午我突然发现了一个很好的几何特性，圆无论从那边投影，始终间距是不变的，那不就是数一数探测器格子的事吗，于是乎敲了几行代码，发现180次投影绝大多数都是29个格子，那就8/29了。然后发现由于椭圆中兴到圆心的距离知道，而中心峰值最大，那就那么通过这个就可以大致确定角度了，还想了想其实中兴也可以这样算啊，取探测器中间位置为中心线不就行了？几何关系就弄出来了，到底靠不靠谱啊，不敢想，这时大概10.30，唐老进来吹逼了，还以为它能告诉我们一个很好的思路，结果就是瞎逼逼，连CT成像的数字图像处理原理都不知道……radon变换都没提到，间距都不会求。哎感到无望的我决定先去把自己的想法实现.敲了一敲，角度还是大致精确，对着其实角度29.。。。几度一个一个往上跑，可是中兴当取不同连线的时候就不行了.又陷入绝望，2，3两问还要重建图像呢，到了下午也不知道怎么处理，这时一个老师道出了精髓，你怎么就能确定中心一定在256像素呢。瞬间茅塞顿开.再也找不到求中心的方法了，在晚上要结束的时候发现matlab有iradon变换的函数，随便放了108个角度进去，竟然出来图像了，虽然有点倾斜,然后搜了一搜，发现scipy里还是准备好这些图像处理的基本变换函数的，就接着有python试了试，果然也是可以的.

第三天就是各种方法乱上了，因为iradon没有中心这个接口，自己又不会写，反变换像素与原图像素不一致，各种乱搞，发现还原出来的图居然正好将原图缩小的一些，对比探测器长度和和原图宽度，竟然是根号2倍关系,随机明白了，开始截图，还原像素了，就这样度过了惊人的第3天.

第四天分析精度再也弄不出来什么东西了，只能xjb搞了，老师说几何方法在4问肯定没法做，那我也没办法了……可是老师说枚举180个角度确定中心去最小误差的搜索方法我至今怀疑，怎么可能一个角度就能搜出中心，而且，一次iradon变换运行时间至少为1秒啊，（可能是搞ACM的原因本能觉得是不可行的）。

在晚上的时候老唐居然在群里发出了，最恶心的嘲讽，附上截图

估计做A题的人的心得凉透了….就这样在这种嘲讽中忙活过了最后的4小时改论文时间…..

至此，数学建模经历全剧终。其实从暑假培训起就有点不想搞了，培训就xbb，什么也没讲到，原理从来不讲，就叫你用matlab命令搞一搞就行了,全程用python的我根本用不到这些东西的，而且现在python社区如此庞大，本科生搞搞小科研个人认为完全够用了，不知道学校为何还要让大家集体用盗版的matlab。（机房也是盗版的）.记得当时听见有老师说建模是他们的事情，我们所要做的就是编程实现,,,,,,这次估计打脸了，在第二天，大佬唐,(就是那个最后挖苦我们的人)说我们知道的也就那么多了，其余我们也不会了….然后一直xbb。。