SimWe仿真论坛 » H20：MathCAD Maple » 请有2.0以上版本的朋友帮个忙

Starriver 发表于 2008-6-27 14:52

请有2.0以上版本的朋友帮个忙

花了不少时间写了程序，却因为手上没有2.0版而无法运行。。。甚是不爽啊

请shamohu版主或其他热心的朋友帮忙运行下、看能否运行即可

[[i] 本帖最后由 Starriver 于 2008-6-27 14:57 编辑 [/i]]

20wangz 发表于 2008-6-27 21:51

你的程序可能有问题。

====== 结果 ======

迭代数: 16
计算用时(时:分:秒:毫秒): 00:00:02:383
计算中止原因: 达到收敛判定标准
优化算法: 全局麦夸特法(Global Levenberg-Marquardt)
目标函数值(最小): 0
delta0_0: 0
delta0_1: 0
delta0_2: 0
delta0_3: 0
delta1_0: 0
delta1_1: 0
delta1_2: 0
delta1_3: 0
delta2_0: 0
delta2_1: 0
delta2_2: 0
delta2_3: 0
delta3_0: 0
delta3_1: 0
delta3_2: 0
delta3_3: 0

====== 计算结束 ======

shamohu 发表于 2008-6-28 00:13

论坛应以答疑解难、探讨问题为主，不应沦为代人做题之处，下不为例！

代码改了一下：
Constant M=3;    //Number of Nodes - 1
Constant R(0:M) = [40, 0, 0, -40];
Constant d(0:M, 0:M)= [ 0, 10, 20, 30,
                       10,  0, 10, 20,
                       20, 10 , 0, 10,
                       30, 20, 10 , 0 ];
Constant Pct = 0.0982, Pcr = 0.1125;
Constant kapa = 3.5, S = 200, Po = 0.398e-6;
Parameter delta(0:M,0:M)=[0,];
Algorithm = SM3;
Minimum;
StartProgram;
Const N=3;         //Conresponding to M
Var i,j : integer;
    Pir, Pit : Double;
    tr, tt : Double;
    Pavg, temT : Array [0..N] of double;
    Pt : Array [0..N] of Array [0..N] of double;
Begin
    for i := 0 to N do
       for j := 0 to N do
          Pt[i,j] := Po* Power(d[i,j],kapa);
    for i := 0 to N do
      Begin
        Pir := 0;
        Pit := 0;
        for j:= 0 to N do
           Pir := Pir + Pcr*delta[j,i];
        for j:= 0 to N do
           Pit := Pit + (Pct + Pt[i,j])*delta[i,j];
        Pavg【i】 := Pir + Pit;End;
    for i := 0 to N do
      Begin
        tr := 0;
        tt := 0;
        for j:= 0 to N do tr := tr + delta[j,i];
        for j:= 0 to N do tt := tt + delta[i,j];
        temT【i】 := tt - tr;
      End;
    ConstrainedResult := For(i=0:3)(200*temT【i】 = R【i】);
    FunctionResult := Sum(i=0:3,Pavg)(Pavg) ;
End;
EndProgram;

结果：
目标函数值(最小): 0.0539118383063271
delta0_0: 6.21789359974084E-8
delta0_1: 8.50042715271178E-8
delta0_2: 3.0318720584012E-7
delta0_3: 0.199999632621388
delta1_0: 3.65407293103159E-9
delta1_1: 1.24348630762805E-7
delta1_2: 4.19854156095324E-8
delta1_3: 1.65124256754609E-7
delta2_0: 1.18956907396729E-8
delta2_1: 1.10371125332822E-7
delta2_2: 9.76578182821259E-9
delta2_3: 2.26123770565042E-7
delta3_0: 5.26310163185821E-9
delta3_1: 1.53883484618193E-8
delta3_2: 3.217965137624E-9
delta3_3: 4.39351719256673E-8

[[i] 本帖最后由 shamohu 于 2008-6-28 00:16 编辑 [/i]]

Starriver 发表于 2008-8-16 15:57

误会了！不过谢谢版主~

这段时间在洽谈购买1stopt注册版，所以没能过来。其实我没有麻烦版主做题的意思，只是想确认下是否是版本的问题~~在此感谢shamuho版主热心的回复及2楼的朋友！:handshake

查看完整版本: 请有2.0以上版本的朋友帮个忙