24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gielcheen

木虫 (正式写手)

应助: 16 (小学生)
金币: 8255
散金: 758
红花: 15
帖子: 675
在线: 560.1小时
虫号: 883648
注册: 2009-10-25
专业: 生物医用高分子材料

[求助] 如何做出拟合曲线求值

我在做PLGA载药微球体外释放时，需要进行体外释放的建模，公式如下所示：等式左边是已知的，当做y，等式右边t也是已知的，当做x，R也是已知的，是常数，那么我现在有一组x和y的数据，我想通过拟合求出D值来，请问用什么软件可以做？那么该怎么做？等式右边的求和不知道能求出什么结果？有人做过吗？谢谢！

回复此楼

» 猜你喜欢

垃圾破二本职称评审标准已经有19人回复
职称评审没过，求安慰已经有53人回复
毕业后当辅导员了，天天各种学生超烦已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复

1楼 2012-02-27 13:44:04

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

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

gielcheen(金币+50): ★★★很有帮助非常感谢，但是编程没有看懂。况且n是到正无穷，如果只到1200，还是不太好！ 2012-02-28 16:44:11

用下列编程模式代码(n= 1200)：

CODE:

Parameter a;

StartRange = [1E-100,1e-10];

Variable  x,y;

StartProgram [Pascal];

Procedure MainModel;

var i, n: integer;

    td: double;

Begin

    for i := 0 to DataLength - 1 do begin

        td := 0;

        for n := 1 to 1200 do

            td := td + exp(-n^2*pi^2*a*x[i]/0.000001);

        y[i] := 1/pi^2*td;

    end;

End;

EndProgram;

data;

0        0.98501

0.21        29.86622

0.42        32.67416

1        36.97943

4        43.80644

7        48.87115

10        51.80471

14        53.65742

18        57.08069

22        59.32729

26        60.3525

30        61.04827

35        62.75492

42        65.76456

49        68.92578

60        74.15131

参数最佳估算
-------------------- -------------
a 1.75957508897304E-14

====== 输出结果 =====

文件：数据文件 - 1
No 目标y 计算y
1 0.98501 121.585420370805
2 29.86622 119.487562839931
3 32.67416 117.454445286095
4 36.97943 112.156261269926
5 43.80644 90.7639936726885
6 48.87115 76.4042344209492
7 51.80471 66.3661814206772
8 53.65742 57.0694131575808
9 57.08069 50.6002923996484
10 59.32729 45.8466101914636
11 60.3525 42.1938718181288
12 61.04827 39.2848251790096
13 62.75492 36.3697677962097
14 65.76456 33.1972842183255
15 68.92578 30.7310720717936
16 74.15131 27.7666907960382