楼主: 冰灵吹不散
1679 0

[问答] 实在是小白帖子了,期望大家不要见笑,matlab中求和修改成求最大 [推广有奖]

  • 0关注
  • 0粉丝

小学生

71%

还不是VIP/贵宾

-

威望
0
论坛币
10 个
通用积分
0
学术水平
1 点
热心指数
1 点
信用等级
0 点
经验
225 点
帖子
8
精华
0
在线时间
3 小时
注册时间
2014-3-18
最后登录
2014-4-14

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
就是将Tval(i)的最后的sum改成max,其余不变,因为是多元,所以不知道怎么运行,这是其中的一部分程序,求大家给个提示什么的,方法

climit = zeros(1, kmax);
    [Tval] = deal(climit);
    runlength = kmax;
    weight = (win(i)-1):(-1):0;
    weight = (1-lambda).^weight';
    wm = repmat(weight, 1, dim)
    Tval(i) = sum((sum(wm.*(rankMat((cidx-win(i)+1):cidx, :)-(cidx+1)/2), 1)).^2)/(sd(i)^2)
    icount = 1;
    nbound=10000;
    permT = zeros(1, nbound);
    permT1 = zeros(1,kmax);
    if i == 1
        while icount <= nbound
            per = randperm(cidx)
            pRank = rankMat(per,:)
            permT(icount) = sum((sum(wm.*(pRank((cidx-win(i)+1):cidx, :)-(cidx+1)/2), 1)).^2)/(sd(i)^2);
            icount = icount + 1
        end
        climit(i) = quantile(permT, 1-alpha)
   else
        p=1
       while p<2
           per = randperm(cidx)
           dsample=rawSample(per,:)
           pRank=rankMat(per,:)
           k=max(1,i-win(i)+1)
           q=1
           while q<2
             if k>i-1;
                 weight = (win(i)-1):(-1):0;
                 weight = (1-lambda).^weight';
                 wm = repmat(weight, 1, dim);
                 climit(i) = sum((sum(wm.*(pRank((cidx-win(i)+1):cidx, :)-(cidx+1)/2), 1)).^2)/(sd(i)^2)
                 p=2
                 break
             else
                 for r = 1 : dim
                     esample = dsample(1:k+m0,:);
                     [b,rdx] = sort(esample(1:k+m0,r));
                     kRank(rdx, r) = 1:k+m0
                 end
                     weight = (win(k)-1):(-1):0;
                     weight = (1-lambda).^weight';
                     wm = repmat(weight, 1, dim)
                     permT1(k) = sum((sum(wm.*(kRank((m0+k-win(k)+1):m0+k, :)-(m0+k+1)/2), 1)).^2)/(sd(k)^2)
                     v=climit
                     if permT1(k)>climit(k)
                        q=1
                        k=k+1   
                     else
                        p=2
二维码

扫码加我 拉你入群

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

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

关键词:MATLAB matla atlab Mat Atl matlab 程序

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-5-1 13:22