版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

当前主题已经存档。

zjl0421

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 485.5
帖子: 19
在线: 1.8小时
虫号: 543674
注册: 2008-04-11
专业: 工业催化

[交流] 【求助】高难度的拟合！！请高手帮忙，急！急！

如果哪位高手肯帮忙的话，我愿出50-100个金币。
请高手帮忙：拟合式中的K,E,p.参量比较多，这个拟合的难度太大，我matlab刚入门，实在是解决不了这样的问题，而这个方程我着急使用，哪位大虾能够帮忙，小弟感激不尽。谢谢！！总共有四十组数据。方程和数据如下所示：
r1 = K1 * exp ( -E1 / R T ) * C1 ^ P1 *C2 ^ P2 - K2 * exp( -E2 / R T ) *C4 ^ P3 * C7 ^ P4
r2 = K3 * exp ( -E3 / R T ) * C1 ^ P5 * C4 ^ P6 - K4 * exp ( -E4 / R T ) * C3 ^ P7 * C7 ^ P8
r3 = K5 * exp ( -E5 / R T ) * C4 ^ P9 * C9 ^ P10 - K6 * exp ( -E6 / R T ) * C5 ^ P11 * C1 ^ P12 * C6 ^ P13
r4 = K7 * exp ( -E7 / R T ) * C3 ^ P14 - K8 * exp ( -E8 / R T ) * C8 ^ P15 * C6 ^ P16
C1 = 1 - C2 - C3 - C4 - C5 - C6 - C7
dC2/dt = -r1
dC3/dt = r2-r3-r4
dC4/dt = r1-r2-r3
dC5/dt = r3
dC6/dt = r3 + r4
dC7/dt = r1 + r2
dC8/dt = r4
其中： t：时间  T：温度 R:8.314
   所有p为 0.2 - 2 ，初值取1；  k和E不能确定；
总共取了40组数据。数据如下：

时间 t = 0 时

         C1       C2       C3       C4       C5       C6       C7       C8
1组       0.95       0.05       0       0       0       0       0       0
2组       0.95       0.05       0       0       0       0       0       0
3组       0.95       0.05       0       0       0       0       0       0
4组       0.95       0.05       0       0       0       0       0       0
5组       0.95       0.05       0       0       0       0       0       0
6组       0.95       0.05       0       0       0       0       0       0
7组       0.95       0.05       0       0       0       0       0       0
8组       0.95       0.05       0       0       0       0       0       0
9组       0.95       0.05       0       0       0       0       0       0
10组       0.95       0.05       0       0       0       0       0       0
11组       0.95       0.05       0       0       0       0       0       0
12组       0.95       0.05       0       0       0       0       0       0
13组       0.95       0.05       0       0       0       0       0       0
14组       0.95       0.05       0       0       0       0       0       0
15组       0.95       0.05       0       0       0       0       0       0
16组       0.95       0.05       0       0       0       0       0       0
17组       0.95       0.05       0       0       0       0       0       0
18组       0.95       0.05       0       0       0       0       0       0
19组       0.95       0.05       0       0       0       0       0       0
20组       0.95       0.05       0       0       0       0       0       0
21组       0.93       0.07       0       0       0       0       0       0
22组       0.93       0.07       0       0       0       0       0       0
23组       0.93       0.07       0       0       0       0       0       0
24组       0.93       0.07       0       0       0       0       0       0
25组       0.93       0.07       0       0       0       0       0       0
26组       0.93       0.07       0       0       0       0       0       0
27组       0.93       0.07       0       0       0       0       0       0
28组       0.93       0.07       0       0       0       0       0       0
29组       0.93       0.07       0       0       0       0       0       0
30组       0.93       0.07       0       0       0       0       0       0
31组       0.93       0.07       0       0       0       0       0       0
32组       0.93       0.07       0       0       0       0       0       0
33组       0.93       0.07       0       0       0       0       0       0
34组       0.93       0.07       0       0       0       0       0       0
35组       0.93       0.07       0       0       0       0       0       0
36组       0.93       0.07       0       0       0       0       0       0
37组       0.93       0.07       0       0       0       0       0       0
38组       0.93       0.07       0       0       0       0       0       0
39组       0.93       0.07       0       0       0       0       0       0
40组       0.93       0.07       0       0       0       0       0       0

时间 t = t 时

      C1       C2       C3       C4       C5       C6       C7       C8
1组    0.9237       0.0017       0.0016       0.0568       0    0.0013       0.0137       0.0013
2组       0.9231       0.0017       0.0035       0.0548       0.0002       0.0014       0.0139       0.0014
3组       0.9226       0.0015       0.0047       0.0536       0.0005       0.0016       0.0142       0.0014
4组       0.9217       0.0016       0.0074       0.0507       0.0007       0.0017       0.0147       0.0015
5组       0.9194       0.0014       0.0138       0.0445       0.0012       0.0021       0.016       0.0015
6组       0.9229       0.0021       0.0023       0.0547       0.0005       0.0016       0.0145       0.0014
7组       0.9209       0.0016       0.0094       0.0489       0.0007       0.0018       0.0152       0.0015
8组       0.9199       0.0015       0.011       0.0474       0.0009       0.0019       0.0159       0.0015
9组       0.9186       0.0013       0.0144       0.0442       0.0011       0.0021       0.0166       0.0017
10组       0.9164       0.001       0.022       0.0376       0.0015       0.0024       0.0174       0.0018
11组       0.9214       0.0032       0.0096       0.046       0.0009       0.002       0.0152       0.0016
12组       0.9207       0.0031       0.0105       0.0449       0.0012       0.0022       0.0156       0.0017
13组       0.9199       0.0029       0.0122       0.0434       0.0013       0.0023       0.0162       0.0018
14组       0.9183       0.0027       0.0159       0.04       0.0015       0.0025       0.0171       0.0019
15组       0.9162       0.0021       0.0235       0.0336       0.0018       0.0029       0.0177       0.0021
16组       0.9205       0.0029       0.0099       0.0454       0.0011       0.0025       0.0158       0.002
17组       0.9195       0.0026       0.0119       0.0434       0.0014       0.0027       0.0164       0.0022
18组       0.9188       0.0025       0.0133       0.0418       0.0016       0.0029       0.0168       0.0022
19组       0.9171       0.0022       0.0173       0.0381       0.0019       0.0032       0.0178       0.0023
20组       0.9157       0.0021       0.0231       0.0323       0.0026       0.0035       0.0184       0.0024
21组       0.8932       0.0025       0.0032       0.0779       0.0005       0.0019       0.0191       0.0017
22组       0.8919       0.0019       0.0056       0.0761       0.0007       0.0021       0.0198       0.0018
23组       0.891       0.0017       0.0085       0.0734       0.0008       0.0023       0.0202       0.002
24组       0.8888       0.0016       0.0135       0.0689       0.001       0.0027       0.0213       0.0022
25组       0.8851       0.0013       0.025       0.0574       0.0023       0.0034       0.0233       0.0023
26组       0.8924       0.0024       0.0054       0.0756       0.0006       0.0021       0.0195       0.002
27组       0.8911       0.0021       0.0076       0.0736       0.0008       0.0024       0.0203       0.002
28组       0.8908       0.002       0.0088       0.0722       0.0011       0.0026       0.0204       0.0021
29组       0.8886       0.0019       0.0146       0.0665       0.0015       0.003       0.0215       0.0023
30组       0.8848       0.0016       0.0257       0.0555       0.0027       0.0037       0.0236       0.0024
31组       0.8916       0.0029       0.006       0.0739       0.0006       0.0025       0.0203       0.0022
32组       0.8911       0.0028       0.0077       0.072       0.0009       0.0027       0.0204       0.0023
33组       0.8905       0.0025       0.0101       0.0696       0.0012       0.003       0.0205       0.0025
34组       0.8882       0.0021       0.016       0.064       0.0018       0.0035       0.0217       0.0027
35组       0.8847       0.0016       0.025       0.0545       0.0031       0.0043       0.0239       0.0029
36组       0.8911       0.0029       0.0085       0.0711       0.0007       0.0029       0.0203       0.0026
37组       0.8904       0.0025       0.01       0.0692       0.0013       0.0033       0.0206       0.0027
38组       0.8898       0.0021       0.0132       0.0654       0.0022       0.0038       0.0207       0.0028
39组       0.8881       0.0021       0.0168       0.0612       0.0029       0.0042       0.0219       0.0028
40组       0.8846       0.0016       0.0249       0.0536       0.0035       0.0047       0.024       0.0031

各组数据所对应的时间 t 和温度 T

      时间t       温度T
1组       600       443.15
2组       750       443.15
3组       1000       443.15
4组       1500       443.15
5组       3000       443.15
6组       600       448.15
7组       750       448.15
8组       1000       448.15
9组       1500       448.15
10组       3000       448.15
11组       600       453.15
12组       750       453.15
13组       1000       453.15
14组       1500       453.15
15组       3000       453.15
16组       600       458.15
17组       750       458.15
18组       1000       458.15
19组       1500       458.15
20组       3000       458.15
21组       600       443.15
22组       750       443.15
23组       1000       443.15
24组       1500       443.15
25组       3000       443.15
26组       600       448.15
27组       750       448.15
28组       1000       448.15
29组       1500       448.15
30组       3000       448.15
31组       600       453.15
32组       750       453.15
33组       1000       453.15
34组       1500       453.15
35组       3000       453.15
36组       600       458.15
37组       750       458.15
38组       1000       458.15
39组       1500       458.15
40组       3000       458.15

[ Last edited by kuhailangyu on 2009-4-15 at 15:34 ]

回复此楼

» 猜你喜欢

266求调剂已经有7人回复
材料调剂已经有16人回复
求调剂：085600材料与化工，考材科基，总分319 已经有22人回复
085600，320分求调剂已经有4人回复
085701求调剂已经有3人回复
考研调剂已经有3人回复
0710生物学336分求调剂已经有6人回复
电子信息调剂交叉学科有推荐吗已经有7人回复
一志愿北京化工大学，初试成绩350求调剂已经有3人回复
280求调剂已经有20人回复

1楼 2009-04-15 15:03:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

应助: 43 (小学生)
贵宾: 10.33
金币: 21115.1
散金: 737
红花: 22
帖子: 3777
在线: 639.7小时
虫号: 293038
注册: 2006-11-04
专业: 计算机应用技术
管辖: 仿真模拟

哇，兄弟的数据够庞大，期待高手来帮忙

赞一下

回复此楼

[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&amp;amp;say=各位版主辛苦了！[/img][/url][/center]

2楼2009-04-15 23:10:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fspdlh

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 345.2
红花: 3
帖子: 522
在线: 4小时
虫号: 578902
注册: 2008-07-04

★ ★ ★
sunxiao(金币+3,VIP+0):谢谢参与，欢迎常来仿真编程版 4-16 01:10

看了一下，对于拟合来说，前四个方程就够了，MATLAB里拟合功能是很强大的，不过我只搞过拟合方程，没拟合过方程组。

我只把每组数据的r求出来了，您的数据是每五个一组的，加上初值共六个，可以用数值求导求出各dC/dt，方法用五点法(我也发过帖子的，http://muchong.com/bbs/viewthread.php?tid=1228217 )，然后用线性代数的方法求r。

至于拟合，期待高手！

[ Last edited by fspdlh on 2009-4-16 at 10:46 ]

赞一下(10人)

回复此楼

3楼2009-04-16 00:13:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fspdlh

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 345.2
红花: 3
帖子: 522
在线: 4小时
虫号: 578902
注册: 2008-07-04

★ ★
sunxiao(金币+2,VIP+0):谢谢参与，欢迎常来仿真编程版 4-16 23:42

求r的程序如下

把t=0和t=t的数据分别存成data0.txt和datat.txt

运行程序

function y=readdata

%readdata
data0=load('data0.txt');
datat=load('datat.txt');
data0=[zeros(40,1) datat(:,end) data0];
datat=[datat(:,end-1:end) datat(:,1:end-2)];
n=0;
data=[];
for i=1:8
data=[data;data0(i*5,

;datat(i*5-4:i*5,

];
end
t=data(:,1)
T=data(:,2)
C=data(:,3:end)

%diff
diffc=zeros(48,8);
for i=1:8
for j=1:8
      diffc(i*6-5:i*6,j)=NumDiff(data(i*6-5:i*6,1),data(i*6-5:i*6,j+2));
end
end
Dc=diffc

%Ar=diffc
A=[ -1  0 0 0 ;... %c2
0 1 -1  -1  ;... %c3
1 -1  -1  0 ;... %c4
0 0 1 0 ;... %c5
0 0 1 1 ;... %c6
1 1 0 0 ;... %c7
0 0 0 1 ]; %c8
r=zeros(48,4);
for i=1:48
r(i,

=[A\[diffc(i,2:end)]']';
end
r

这样就可以得到拟合所要用到的48组数据，另外数据中各C值都有，我就没有用C1 = 1 - C2 - C3 - C4 - C5 - C6 - C7这个式子

[ Last edited by fspdlh on 2009-4-16 at 10:47 ]

赞一下(10人)

回复此楼

4楼2009-04-16 10:26:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fspdlh

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 345.2
红花: 3
帖子: 522
在线: 4小时
虫号: 578902
注册: 2008-07-04

★ ★
sunxiao(金币+2,VIP+0):谢谢参与，欢迎常来仿真编程版 4-16 23:42

t =

         0
      600
      750
      1000
      1500
      3000
         0
      600
      750
      1000
      1500
      3000
         0
      600
      750
      1000
      1500
      3000
         0
      600
      750
      1000
      1500
      3000
         0
      600
      750
      1000
      1500
      3000
         0
      600
      750
      1000
      1500
      3000
         0
      600
      750
      1000
      1500
      3000
         0
      600
      750
      1000
      1500
      3000

T =

  443.1500
  443.1500
  443.1500
  443.1500
  443.1500
  443.1500
  448.1500
  448.1500
  448.1500
  448.1500
  448.1500
  448.1500
  453.1500
  453.1500
  453.1500
  453.1500
  453.1500
  453.1500
  458.1500
  458.1500
  458.1500
  458.1500
  458.1500
  458.1500
  443.1500
  443.1500
  443.1500
  443.1500
  443.1500
  443.1500
  448.1500
  448.1500
  448.1500
  448.1500
  448.1500
  448.1500
  453.1500
  453.1500
  453.1500
  453.1500
  453.1500
  453.1500
  458.1500
  458.1500
  458.1500
  458.1500
  458.1500
  458.1500

C =

0.9500 0.0500       0       0       0       0       0       0
0.9237 0.0017 0.0016 0.0568       0 0.0013 0.0137 0.0013
0.9231 0.0017 0.0035 0.0548 0.0002 0.0014 0.0139 0.0014
0.9226 0.0015 0.0047 0.0536 0.0005 0.0016 0.0142 0.0014
0.9217 0.0016 0.0074 0.0507 0.0007 0.0017 0.0147 0.0015
0.9194 0.0014 0.0138 0.0445 0.0012 0.0021 0.0160 0.0015
0.9500 0.0500       0       0       0       0       0       0
0.9229 0.0021 0.0023 0.0547 0.0005 0.0016 0.0145 0.0014
0.9209 0.0016 0.0094 0.0489 0.0007 0.0018 0.0152 0.0015
0.9199 0.0015 0.0110 0.0474 0.0009 0.0019 0.0159 0.0015
0.9186 0.0013 0.0144 0.0442 0.0011 0.0021 0.0166 0.0017
0.9164 0.0010 0.0220 0.0376 0.0015 0.0024 0.0174 0.0018
0.9500 0.0500       0       0       0       0       0       0
0.9214 0.0032 0.0096 0.0460 0.0009 0.0020 0.0152 0.0016
0.9207 0.0031 0.0105 0.0449 0.0012 0.0022 0.0156 0.0017
0.9199 0.0029 0.0122 0.0434 0.0013 0.0023 0.0162 0.0018
0.9183 0.0027 0.0159 0.0400 0.0015 0.0025 0.0171 0.0019
0.9162 0.0021 0.0235 0.0336 0.0018 0.0029 0.0177 0.0021
0.9500 0.0500       0       0       0       0       0       0
0.9205 0.0029 0.0099 0.0454 0.0011 0.0025 0.0158 0.0020
0.9195 0.0026 0.0119 0.0434 0.0014 0.0027 0.0164 0.0022
0.9188 0.0025 0.0133 0.0418 0.0016 0.0029 0.0168 0.0022
0.9171 0.0022 0.0173 0.0381 0.0019 0.0032 0.0178 0.0023
0.9157 0.0021 0.0231 0.0323 0.0026 0.0035 0.0184 0.0024
0.9300 0.0700       0       0       0       0       0       0
0.8932 0.0025 0.0032 0.0779 0.0005 0.0019 0.0191 0.0017
0.8919 0.0019 0.0056 0.0761 0.0007 0.0021 0.0198 0.0018
0.8910 0.0017 0.0085 0.0734 0.0008 0.0023 0.0202 0.0020
0.8888 0.0016 0.0135 0.0689 0.0010 0.0027 0.0213 0.0022
0.8851 0.0013 0.0250 0.0574 0.0023 0.0034 0.0233 0.0023
0.9300 0.0700       0       0       0       0       0       0
0.8924 0.0024 0.0054 0.0756 0.0006 0.0021 0.0195 0.0020
0.8911 0.0021 0.0076 0.0736 0.0008 0.0024 0.0203 0.0020
0.8908 0.0020 0.0088 0.0722 0.0011 0.0026 0.0204 0.0021
0.8886 0.0019 0.0146 0.0665 0.0015 0.0030 0.0215 0.0023
0.8848 0.0016 0.0257 0.0555 0.0027 0.0037 0.0236 0.0024
0.9300 0.0700       0       0       0       0       0       0
0.8916 0.0029 0.0060 0.0739 0.0006 0.0025 0.0203 0.0022
0.8911 0.0028 0.0077 0.0720 0.0009 0.0027 0.0204 0.0023
0.8905 0.0025 0.0101 0.0696 0.0012 0.0030 0.0205 0.0025
0.8882 0.0021 0.0160 0.0640 0.0018 0.0035 0.0217 0.0027
0.8847 0.0016 0.0250 0.0545 0.0031 0.0043 0.0239 0.0029
0.9300 0.0700       0       0       0       0       0       0
0.8911 0.0029 0.0085 0.0711 0.0007 0.0029 0.0203 0.0026
0.8904 0.0025 0.0100 0.0692 0.0013 0.0033 0.0206 0.0027
0.8898 0.0021 0.0132 0.0654 0.0022 0.0038 0.0207 0.0028
0.8881 0.0021 0.0168 0.0612 0.0029 0.0042 0.0219 0.0028
0.8846 0.0016 0.0249 0.0536 0.0035 0.0047 0.0240 0.0031

Dc =

  1.0e-003 *

-0.0757 -0.1449       0 0.1811       0 0.0034 0.0400 0.0034
-0.0067 -0.0013 0.0062 -0.0004 0.0002 0.0009 0.0024 0.0009
-0.0027 -0.0001 0.0081 -0.0086 0.0014 0.0007 0.0012 0.0001
-0.0019 -0.0001 0.0050 -0.0051 0.0007 0.0004 0.0011 0.0000
-0.0017 0.0000 0.0049 -0.0050 0.0004 0.0002 0.0009 0.0000
-0.0013 -0.0004 0.0034 -0.0029 0.0003 0.0003 0.0008 0.0000
-0.0706 -0.1410 -0.0000 0.1950 0.0004 0.0037 0.0398 0.0037
-0.0188 -0.0065 0.0147 -0.0047 0.0011 0.0017 0.0070 0.0010
-0.0068 -0.0009 0.0177 -0.0159 0.0011 0.0007 0.0035 0.0001
-0.0032 -0.0004 0.0064 -0.0060 0.0006 0.0004 0.0020 0.0000
-0.0020 -0.0003 0.0060 -0.0054 0.0003 0.0003 0.0008 0.0001
-0.0006 -0.0000 0.0038 -0.0029 0.0002 0.0001 -0.0000 0.0000
-0.0821 -0.1399 0.0240 0.1439 0.0011 0.0049 0.0435 0.0043
-0.0077 -0.0023 0.0079 0.0003 0.0019 0.0018 0.0043 0.0010
-0.0038 -0.0007 0.0062 -0.0071 0.0008 0.0007 0.0025 0.0005
-0.0032 -0.0006 0.0071 -0.0063 0.0004 0.0004 0.0021 0.0003
-0.0021 -0.0004 0.0063 -0.0055 0.0003 0.0003 0.0008 0.0002
-0.0000 -0.0004 0.0033 -0.0024 0.0001 0.0002 0.0000 0.0001
-0.0832 -0.1397 0.0190 0.1469 0.0017 0.0064 0.0442 0.0049
-0.0106 -0.0044 0.0146 -0.0006 0.0020 0.0018 0.0062 0.0018
-0.0041 -0.0008 0.0085 -0.0100 0.0013 0.0010 0.0024 0.0001
-0.0030 -0.0005 0.0064 -0.0068 0.0007 0.0007 0.0017 0.0000
-0.0017 -0.0002 0.0056 -0.0054 0.0005 0.0003 0.0008 0.0001
-0.0000 0.0000 0.0008 -0.0012 0.0004 0.0000 0.0000 0.0000
-0.1035 -0.1993 -0.0000 0.2433 0.0004 0.0046 0.0536 0.0046
-0.0136 -0.0081 0.0098 0.0007 0.0011 0.0017 0.0073 0.0010
-0.0053 -0.0015 0.0142 -0.0121 0.0007 0.0010 0.0025 0.0007
-0.0039 -0.0003 0.0108 -0.0099 0.0004 0.0008 0.0018 0.0006
-0.0034 -0.0002 0.0089 -0.0084 0.0005 0.0006 0.0018 0.0001
-0.0010 -0.0002 0.0059 -0.0067 0.0012 0.0002 0.0007       0
-0.1059 -0.2012 0.0045 0.2375 0.0007 0.0047 0.0542 0.0060
-0.0138 -0.0050 0.0124 0.0002 0.0012 0.0024 0.0083 0.0001
-0.0026 -0.0008 0.0083 -0.0093 0.0013 0.0012 0.0011 0.0000
-0.0018 -0.0003 0.0065 -0.0072 0.0010 0.0008 0.0007 0.0004
-0.0034 -0.0002 0.0095 -0.0094 0.0008 0.0006 0.0018 0.0001
-0.0011 -0.0002 0.0043 -0.0043 0.0008 0.0002 0.0008 0.0000
-0.1125 -0.2008 0.0089 0.2318 0.0002 0.0064 0.0604 0.0061
-0.0061 -0.0029 0.0109 0.0004 0.0015 0.0018 0.0016 0.0010
-0.0027 -0.0008 0.0105 -0.0118 0.0016 0.0012 0.0005 0.0007
-0.0031 -0.0010 0.0105 -0.0102 0.0012 0.0011 0.0007 0.0006
-0.0033 -0.0005 0.0085 -0.0086 0.0010 0.0007 0.0019 0.0002
-0.0006       0 0.0017 -0.0027 0.0006 0.0002 0.0008 -0.0000
-0.1130 -0.1992 0.0175 0.2234 -0.0000 0.0066 0.0593 0.0073
-0.0081 -0.0060 0.0112 0.0005 0.0023 0.0032 0.0036 0.0010
-0.0032 -0.0020 0.0111 -0.0145 0.0040 0.0023 0.0008 0.0005
-0.0027 -0.0001 0.0097 -0.0114 0.0022 0.0012 0.0007 0.0000
-0.0029 -0.0000 0.0063 -0.0066 0.0007 0.0005 0.0019 0.0000
-0.0015 -0.0006 0.0041 -0.0026 0.0000 0.0000 0.0007 0.0004

r =

  1.0e-003 *

0.1181 -0.0491 -0.0118 -0.0063
0.0007 0.0031 -0.0012 0.0004
-0.0022 0.0057 0.0005 -0.0009
-0.0012 0.0036 0.0002 -0.0005
-0.0013 0.0035 0.0001 -0.0004
-0.0006 0.0023 0.0000 -0.0003
0.1203 -0.0536 -0.0149 -0.0054
0.0026 0.0086 -0.0008 -0.0006
-0.0039 0.0119 -0.0002 -0.0016
-0.0012 0.0046 0.0001 -0.0005
-0.0015 0.0039 -0.0001 -0.0005
-0.0010 0.0021 -0.0002 -0.0004
0.1056 -0.0279 -0.0104 -0.0073
0.0021 0.0040 -0.0006 0.0000
-0.0012 0.0052 0.0004 -0.0002
-0.0013 0.0050 -0.0001 -0.0004
-0.0015 0.0040 -0.0002 -0.0005
-0.0007 0.0018 -0.0002 -0.0003
0.1070 -0.0299 -0.0096 -0.0061
0.0027 0.0071 -0.0018 -0.0001
-0.0020 0.0067 0.0008 -0.0008
-0.0014 0.0048 0.0004 -0.0005
-0.0015 0.0038 0.0000 -0.0005
-0.0003 0.0006 0.0002 -0.0002
0.1604 -0.0663 -0.0149 -0.0091
0.0052 0.0054 -0.0005 -0.0003
-0.0028 0.0094 -0.0004 -0.0008
-0.0027 0.0074 -0.0004 -0.0004
-0.0022 0.0061 -0.0001 -0.0006
-0.0017 0.0041 0.0006 -0.0009
0.1593 -0.0625 -0.0149 -0.0089
0.0042 0.0067 -0.0010 -0.0004
-0.0022 0.0059 0.0008 -0.0009
-0.0020 0.0046 0.0004 -0.0005
-0.0024 0.0067 0.0001 -0.0008
-0.0010 0.0029 0.0003 -0.0006
0.1594 -0.0568 -0.0147 -0.0080
0.0011 0.0039 -0.0015 -0.0004
-0.0033 0.0073 0.0006 -0.0008
-0.0028 0.0069 0.0002 -0.0008
-0.0019 0.0061 0.0003 -0.0007
-0.0005 0.0016 0.0005 -0.0003
0.1554 -0.0515 -0.0158 -0.0078
0.0032 0.0046 -0.0006 -0.0004
-0.0031 0.0083 0.0025 -0.0017
-0.0032 0.0068 0.0011 -0.0013
-0.0015 0.0048 0.0003 -0.0005
-0.0006 0.0024 -0.0004 -0.0001

[ Last edited by fspdlh on 2009-4-16 at 16:09 ]

赞一下(10人)

回复此楼

5楼2009-04-16 10:27:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fspdlh

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 345.2
红花: 3
帖子: 522
在线: 4小时
虫号: 578902
注册: 2008-07-04

呵呵，搞错了，留个楼层吧

[ Last edited by fspdlh on 2009-4-16 at 16:20 ]

赞一下

回复此楼

6楼2009-04-16 16:06:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fspdlh

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 345.2
红花: 3
帖子: 522
在线: 4小时
虫号: 578902
注册: 2008-07-04

★ ★ ★ ★ ★
sunxiao(金币+5,VIP+0):谢谢参与，欢迎常来仿真编程版 4-16 23:42

昨天一晚加今天一天，终于把该编的都编好了（金币方面应该不会亏待吧，呵呵），lz下载后在exe的开头改各变量的初值，运行
>>[K E P R]=exe;
后看看四个R值，越大说明越好，我试了一些，只有第一个方程拟合得比较好好，其它都不怎么样。

下面是我搞出来的一组。
K =

  1.0e+009 *

0.0000 0.0000 0.0000 -0.0000 0.0000 0.0000 2.6948 0.0000

ans =

  6.4865e-004

ans =

2.5796

ans =

0.0010

ans =

-319.7096

ans =

0.0563

ans =

0.4559

ans =

  2.6948e+009

ans =

0.0029

E =

  1.0e+006 *

-0.0069 0.0192 0.0086 1.2377 0.0052 -0.0251 0.1072 0.0047

ans =

-6.9280e+003

ans =

  1.9214e+004

ans =

  8.5666e+003

ans =

  1.2377e+006

ans =

  5.2418e+003

ans =

-2.5112e+004

ans =

  1.0721e+005

ans =

  4.6763e+003

P =

  Columns 1 through 10

2.0000 1.1694 0.2000 2.0000 2.0000 0.9804 0.2037 0.6591 1.7216 1.4183

  Columns 11 through 17

2.0000 0.2000 2.0000 2.0000 0.5356 0.6219 1.0147

R =

0.9989 0.8851 0.5456 -0.7254

[ Last edited by fspdlh on 2009-4-16 at 17:03 ]

赞一下(10人)

回复此楼

7楼2009-04-16 16:52:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 zjl0421 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 调剂0855-288 +5	x熊二a 2026-04-03	5/250	2026-04-04 00:19 by 猪会飞
[考研] 化学调剂求助 +6	LULONG1 2026-04-03	6/300	2026-04-03 23:13 by qzxyhcsy
[考研] 387求调剂 +3	爱吃片豆土 2026-04-03	4/200	2026-04-03 23:05 by 小小树2024
[考研] 322求调剂 +6	FZAC123 2026-04-03	6/300	2026-04-03 22:23 by 科研小专家
[考研] 材料科学与工程339求调剂 +12	hyz0119 2026-03-31	13/650	2026-04-03 18:33 by ls刘帅
[考研] 五邑大学土木工程招调剂生2026 +3	wyutj 2026-03-31	4/200	2026-04-03 18:21 by zengxj_7201
[考研] 266分，求材料相关专业调剂 +13	哇呼哼呼哼 2026-03-30	15/750	2026-04-03 15:24 by arrow8852
[考研] 0854求调剂 +3	assdll 2026-04-03	3/150	2026-04-03 14:32 by fxue1114
[考研] 338求调剂，一志愿能源动力，外语是日语203 +5	zzz，，r 2026-04-02	5/250	2026-04-03 09:45 by 蓝云思雨
[考研] 一志愿山东大学化学与化工学院材料与化工专硕，360分求调剂 +4	不愿透露姓名的� 2026-04-02	4/200	2026-04-03 09:29 by 遗忘消失的灆
[考研] 327求调剂 +9	汲亦昊 2026-03-29	9/450	2026-04-02 11:44 by 运气yunqi
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +10	哇呼哼呼哼 2026-04-01	11/550	2026-04-01 21:48 by chyhaha
[考研] 379求调剂 +3	?苦瓜不苦 2026-04-01	3/150	2026-04-01 20:09 by 暮云清寒
[考研] 085600，320分求调剂 +5	大馋小子 2026-04-01	6/300	2026-04-01 19:40 by 唐沐儿
[考研] 一志愿南京航空航天大学，080500材料科学与工程学硕 +7	@taotao 2026-03-30	7/350	2026-04-01 14:30 by chenqifeng666
[硕博家园] 考研调剂 +5	骆驼男人 2026-04-01	5/250	2026-04-01 14:28 by syjjj0321
[考研] 求调剂，一志愿北林食品与营养095500，301分，已过六级，有科研经历 +4	快乐储蓄罐 2026-03-31	4/200	2026-04-01 09:26 by JourneyLucky
[考研] 生物考研337分求调剂 +4	cgxin 2026-03-30	6/300	2026-03-31 14:18 by 记事本2026
[考研] 332求92调剂 +8	蕉蕉123 2026-03-28	8/400	2026-03-29 10:46 by 周梓丹
[考研] 312，生物学求调剂 +3	小译同学abc 2026-03-28	3/150	2026-03-28 15:32 by 落睿可思