请问怎么在procedure里面循环定义矩阵，十分谢谢

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

比如，令 x=seqa(1,1,100)，我想在某个procedure内定义一个矩阵M， 使得矩阵的某个元素等于x[j] （这个元素的值随着循环不断改变），其他元素为local 变量或者常数，
即 M= x[j] ~ 20 ~ c (local variable)

程序大体这样：            程序正文
                                    do until j>N;        ( j 为全局变量)
                                      调用 procedure
                                    endo;
                                    end;   (主程序完）

                                    procedure (定义矩阵 M）
                                    local c
                                    M= x[j] ~ 20 ~ c
                                    retp;
                                    endp;

我在前面贴过程序，不知道是大家看晕了，还是对那个问题没兴趣，所以我只能把它简化了再问，如果这个算重复发帖的话，就请版大删了吧

小弟也是不得已而为，两天了，这个问题还不知道怎么debug，真的是很郁闷啊。。。哪位仁兄能帮忙troubleshoot下，小弟感激涕零啊

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Procedure ROCE OCED edu ROC 矩阵 定义 Procedure

C:\222\GDP\gdpp.prg(58) : error G0003 : Indexing a matrix as a vector
c:\gauss7.0\src\pgraph.src(58) : error G0474 : ' 0 = asclabel(xl,yl)' : Illegal creation of global
c:\gauss7.0\src\pgraph.src(59) : error G0474 : 'xl' : Illegal creation of global
c:\gauss7.0\src\pgraph.src(59) : error G0474 : 'xl' : Illegal creation of global
我也是研究了几天，还是没办法解决，非常抱歉帮不了你！

我把它转成matlab程序也不行，真的是太诡异了

不过还是谢谢哈

2# zhangtao

F=    phi1[n]~phi2|
    1~0;
关键问题是出在这里，把这里修改好估计应该可以出结果！
proc getnr(nsets,ncols);
    local nr;

    /* check for complex input */
    if iscplx(nsets);
还有以上三行，也得好好想想，但我估计这三行不会影响出结果。

我也一直想不到gauss里面有什么方法可以避开把phi1[n]当做变量
不过后来我用matlab的fminunc把phi1[n]当做参数，只把mu, var, phi2当做变量求最值给求出来了。。。

不知道gauss里面有没有什么function是可以起到fminunc这种作用的？



4# zhangtao

Gauss中也有同样的函数，好像是fsolve什么的，一时记不起了，但是有，我前几周用过。
我一直是尝试，把
F=    phi1[n]~phi2|
    1~0;
上式中的phi1[n],phi2, 1,0,
全部用矩阵表示，然后在
竖着联结，但是没有成功，
估计是矩阵阶数不匹配。其实我看你的程序是gauss程序，也调用了gauss库函数，
在matlab上运行应该不太容易。

不是，我是把GAUSS程序重写成matlab程序，然后在matlab环境里debug

6# zhangtao