2013年研究生建模A题答案质疑

Solemn胜威

& x/ ]# T9 k& X/ z6 z, j5 A% F
' w. O* ^9 n; q3 p) Y/ t. y

（1）  附录4分别给出了风扇、CDFS,高压压气机的特性数据。利用线性插值法计算出压气机的换算转速prc、压比函数值为zz时的特性图上的增压比prc、效率nc和换算流量Wc。

griddata（），，，x是换算转速，y是压比函数，z就是相应的Prc，nc, Wc

试问，答案给出的风扇和CDFS的结果，是怎么算出来的，？
' d4 D7 v3 i: b6 p0 w4 p
我的程序代码，
) k8 k4 Z7 n3 a2 d2 M7 g, m# W% 线性插值法计算出压气机的换算转速为ncor、压比函数值为zz 时的特性图上的增压比prc、效率nc 和换算流量Wc
2 A& m; |7 Y- s1 m0 v( G! T1 Inint1 = 0:0.01:2;  % 转速
zz1 = 0:0.005:1;   % 压比函数值
[nint2,zz2] = meshgrid(nint1,zz1);
9 k/ |' T  T# k3 [2 ~7 C2 F  h* Xncor_data =[];
ncor_data(:,:,1) =  griddata(nint,zz,prc,nint2,zz2,'v4'); % 样条函数内插  % 增压比
ncor_data(:,:,2) =  griddata(nint,zz,nc,nint2,zz2,'v4');  % 样条函数内插  % 效率
ncor_data(:,:,3) =  griddata(nint,zz,Wc,nint2,zz2,'v4');  % 样条函数内插  % 流量
" V% Y% M2 t, O0 l$ `
2 k" J8 j( G7 k% }8 G/ Y算出来的结果：
% i9 \) z& \) _/ i3 L5 }
进风管输出 总温 = 244.381200
进风管输出 总压 = 0.344739
' e' s# B: ^. [) D0 @
pnW =
/ V5 \5 x3 V: A% V! t: G7 T   prc            效率nc     流量Wc
    3.7834    0.3920     111.7823
" @1 B. q! f  X' Q4 f
风扇输出 总温 = 530.341391
0 A" Q9 v3 J/ \  Y$ l- `风扇输出 总压 = 1.304293
风扇输出 流量 = 41.297327
5 ?& G7 i6 f! w& Z# G; r风扇输出 功 = 289881.677719
# O, P! j* a9 e" C风扇输出 功率 = 11971338.445019

8 Z2 v1 T; S8 ICDFS输出 总温 = 538.744399
4 v' s8 W* t; K1 MCDFS输出 总压 = 1.385129
CDFS输出 流量 = 41.297327
; y+ b* h* j) KCDFS输出 功 = 8722.346390
CDFS输出 功率 = 360209.591341
2 S2 c) w1 w8 I: R/ v' K' ~; d* ^
和官方公布的答案，不一样，请大家解释一下，我试 了各种方法，均不得行。

+ v0 U  O4 J5 F谢谢
7 k( q* i5 E: o" X1 N" v& ]  _
* A  d6 {5 ?+ E3 H! R附上参考答案。

3 r6 i" M$ A" p2 U  w7 z$ I) ^三等奖啊，不甘心啊，请大家帮忙指正。

& e4 _; D2 Z% b) N3 b9 U

Solemn胜威

请大家踊跃指正，提出自己的见解

kendy1992

楼主有B题的吗，谢谢啦

sleepstealer

看不懂。。。帮你顶帖吧

yuanyuanxiang

kendy1992 发表于 2013-11-23 09:13
楼主有B题的吗，谢谢啦

B题参考答案在这里。

17263249

Lz的griddata采用的是v4插值法，参考答案给的是线性的结果。你可以把v4改为linear试试。
7 T0 `& z* {& f' |, Q6 g7 `
本人没做过本题，若说的不对请见谅。

Solemn胜威

17263249 发表于 2013-11-23 21:53
! n* \0 @' X# s" GLz的griddata采用的是v4插值法，参考答案给的是线性的结果。你可以把v4改为linear试试。

! \' u2 j6 n/ n( c! i' c0 B' w! Y本人没做过本题 ...

哥们你这个回答，可以一试，谢谢哈