数学建模社区-数学中国

标题: [公告]把你的问题提出来 [打印本页]

作者: liwenhui    时间: 2008-7-23 20:49
标题: [公告]把你的问题提出来
大家有什么问题都发帖子提出来,也可以到“数学问吧”里提问(推荐)。希望大家支持工作!
作者: 小水一把    时间: 2008-7-23 21:21
[BR]刚跑完步回来···一身大汗就来捧场···
作者: ljm1989    时间: 2008-8-6 21:40
标题: 求助啊!!

 

$ n& V$ b* X! {$ W1 ^& T

我想画一个四维的图形??例如

1 a, @/ A& [) G; K/ D' J

x=[6744 4536 4893.6 8860.8 7507.2 6782.4 9494.4 6216 6432 8380.8 7375.2]';
z=[6.1 8.5 6.4 5.4 6.3 6.7 5.2 6.6 6.9 5.7 5.9]';
w=[20.95 29.6 21.77 21.38 22.19 24.24 25.93 19.52 19.39 21.54 24.44]';
y=[9460.9 19076.5 13846.7 8342.2 10685.1 11403.1 7384.6 12873.5 10524.5 10446.0 12280.6]';

, {5 }2 u( [: t( l, K; F- g

找出y与x,w,z的多元曲面关系

& Q: `$ |" N& u1 {. N) D

怎么办??

/ i; x) y( G. C8 C( W

怎么用插值法或多元回归做呢??

6 P, V1 T' H0 b% N8 \

做出请发邮件给我:

# f( d. Y. ?1 n. X/ |

363167932@qq。com

3 H, g7 L$ r) N; Y* n7 T

谢谢!!!

& Y- T6 ]3 o- \9 Z; R
[此贴子已经被作者于2008-8-6 21:41:28编辑过]

作者: liwenhui    时间: 2008-8-7 20:10

(你应该发一个新帖让网友们帮你解答,不要跟我的帖)

5 g3 g8 |/ s2 m" w: J2 C7 E

多元回归的代码如下:
x=[6744 4536 4893.6 8860.8 7507.2 6782.4 9494.4 6216 6432 8380.8 7375.2]';
z=[6.1 8.5 6.4 5.4 6.3 6.7 5.2 6.6 6.9 5.7 5.9]';
w=[20.95 29.6 21.77 21.38 22.19 24.24 25.93 19.52 19.39 21.54 24.44]';
y=[9460.9 19076.5 13846.7 8342.2 10685.1 11403.1 7384.6 12873.5 10524.5 10446.0 12280.6]';
X=[ones(11,1),x,y,w];
[b,bint,r,rint,s]=regress(y,X,0.5)
得到回归结果:
b =

8 A! d4 s7 g* Z/ K d* W

  1.0e+003 *

, D; W+ Y& M, ~( C1 B, P1 T

   7.63365814703216
  -0.00125510669136
   0.83801013875319
   0.32224150099284

/ L& _1 f _# D A& e2 D, r


bint =

) V, K2 Z6 p9 e! P9 m. ]6 L% P

  1.0e+004 *

5 V, z4 A9 z! E# b

  -1.30439270279949   2.83112433220592
  -0.00025564588870   0.00000462455043
  -0.15661574000751   0.32421776775815
  -0.00588755516686   0.07033585536543

( U1 u# U& P6 _% y5 L6 O, k/ T


r =

& k2 ]( D$ n. h; c

  1.0e+003 *

) K, ~( p/ G4 e+ @+ d' S9 V

  -1.57113991271872
   0.47457119617056
  -0.02343040684476
   0.41501318324491
   0.04377602514276
  -1.14372443728898
  -1.04594801887843
   1.22056403128825
  -1.06684456987565
   1.61340028960900
   1.08376262015110

) w7 n+ |" N- f5 E. ?* C1 z


rint =

2 H- {0 J6 q, U# A

  1.0e+003 *

5 \ f# S# ?( _4 f. x

  -4.18060826996541   1.03832844452798
  -0.58713163609601   1.53627402843712
  -1.63265383009492   1.58579301640540
  -2.41980607307114   3.24983243956095
  -3.01244690932543   3.09999895961095
  -4.02599379571372   1.73854492113576
  -3.03725872250225   0.94536268474539
  -1.33675472238316   3.77788278495966
  -3.27677184264986   1.14308270289856
  -0.92047402141587   4.14727460063388
  -1.61182764620170   3.77935288650389

- y7 t M' h7 x


s =

3 f; |2 ~( k# C

  1.0e+006 *

6 {* k2 x' ]1 v2 w# C6 l, ~4 o" {# P

   0.00000088230440   0.00001749182053   0.00000000124110   1.66800690870350
也就是说,方程为:
y=7633.65814703216x-1.25510669136z+838.01013875319w+322.24150099284

& c1 V" X1 `& ^6 [$ K' l/ d

都不是很显著,因为置信区间都过了0点!

- w6 q9 d* X2 S, f, g" U% \
[此贴子已经被作者于2008-8-9 13:02:27编辑过]

作者: bitwt    时间: 2008-10-27 02:01
有没有接触过线性矩阵不等式(LMI)的。我现在的问题是 一般的lmi可以用matlab来解的,但是我现在解的这个维数比较大,所以在运算的时候总是出现 out of memery的提示,电脑的配置应该很高了。试过了好多办法都不行。请问,有没有好的办法解决这个问题。如果没有的话,还有没有别的软件可以解lmi~~~谢谢各位老师·~~
作者: trisol    时间: 2008-12-20 10:22
提示: 作者被禁止或删除 内容自动屏蔽




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5