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

Solemn胜威

6 ^& O( C. M- y  d$ N4 D3 f
! X( d0 W9 C& q3 F, G

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

griddata（），，，x是换算转速，y是压比函数，z就是相应的Prc，nc, Wc
; D/ S. `/ p* g2 [% B; |" m- S
试问，答案给出的风扇和CDFS的结果，是怎么算出来的，？

! T6 _; O; W, ?9 `- {, p. [我的程序代码，
2 g" l  A% b) x: j& d% 线性插值法计算出压气机的换算转速为ncor、压比函数值为zz 时的特性图上的增压比prc、效率nc 和换算流量Wc
. }( r3 c$ i/ m, R; ?nint1 = 0:0.01:2;  % 转速
zz1 = 0:0.005:1;   % 压比函数值
[nint2,zz2] = meshgrid(nint1,zz1);
ncor_data =[];
ncor_data(:,:,1) =  griddata(nint,zz,prc,nint2,zz2,'v4'); % 样条函数内插  % 增压比
* E. f* U  j$ _' X0 }+ o+ K* T+ Bncor_data(:,:,2) =  griddata(nint,zz,nc,nint2,zz2,'v4');  % 样条函数内插  % 效率
' s$ V/ B3 s5 l/ j+ {! Z( I+ o1 c& Qncor_data(:,:,3) =  griddata(nint,zz,Wc,nint2,zz2,'v4');  % 样条函数内插  % 流量
7 ^3 |2 N7 ~7 i  [7 b6 G4 \
1 t8 y  ^2 |1 e. s% S( M$ K算出来的结果：

$ S/ u, y# [: U, h( [* f6 @2 q& V4 l进风管输出 总温 = 244.381200
3 Z8 Z( A7 u% @# W9 ?. y* [进风管输出 总压 = 0.344739
1 J' Y  i% X5 u7 y9 B' o
pnW =
1 r+ c# x) e& W   prc            效率nc     流量Wc
: F6 x: B2 U- j4 ?    3.7834    0.3920     111.7823
! I) g) T5 o' k' f
风扇输出 总温 = 530.341391
3 H- B0 j- Z3 a+ ~5 I; b风扇输出 总压 = 1.304293
9 M! s$ A( @3 d0 ^4 y风扇输出 流量 = 41.297327
% Q0 \! Y; L+ N$ E风扇输出 功 = 289881.677719
风扇输出 功率 = 11971338.445019

- i  D6 ?* h# x1 R9 a* G3 LCDFS输出 总温 = 538.744399
3 J% c" q. P  T) q" r8 F6 JCDFS输出 总压 = 1.385129
* B9 x3 K& I2 J( OCDFS输出 流量 = 41.297327
CDFS输出 功 = 8722.346390
8 W) }; E$ g- aCDFS输出 功率 = 360209.591341

+ U& u0 c7 @/ T2 X和官方公布的答案，不一样，请大家解释一下，我试 了各种方法，均不得行。

2 k, c/ e& s! X% @谢谢
" z" m2 k! }: ^
附上参考答案。
. a  r8 V1 y. M( c
三等奖啊，不甘心啊，请大家帮忙指正。

Solemn胜威

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

kendy1992

楼主有B题的吗，谢谢啦

sleepstealer

看不懂。。。帮你顶帖吧

yuanyuanxiang

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

B题参考答案在这里。

17263249

Lz的griddata采用的是v4插值法，参考答案给的是线性的结果。你可以把v4改为linear试试。
1 ]$ ^- H1 q4 U* r! n$ x, D
本人没做过本题，若说的不对请见谅。

Solemn胜威

17263249 发表于 2013-11-23 21:53
Lz的griddata采用的是v4插值法，参考答案给的是线性的结果。你可以把v4改为linear试试。

6 |- f% e5 {$ F$ o; ?5 A1 j% ]" ]本人没做过本题 ...

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