24小时热门版块排行榜    

查看: 925  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

quarkhe

新虫 (初入文坛)

[交流] 未定义变量W,求助!!!谢谢 已有2人参与

CODE:
function DP=f2(t,P)
global A N
format long g;
N=15;
S=zeros(1,N);
[m,l]=gausspoint(N-2);
S(1)=-1;
S(N)=1;
for i=2:N-1
    S(i)=m(i-1);
end
S;
X=zeros(1,N);
X(1)=0;
X(2)=3*10^(-4);
X(N-1)=1-5*10^(-5);
X(N)=1;
for i=3:N-2
   X(i)=0.5*S(i)+0.5;
end
%求权系数
A=zeros(N,N,4);
XX=ones(1,N);
for i=1:N
    for j=1:i-1
        XX(i)=(X(i)-X(j))*XX(i);
    end
    for j=i+1:N
        XX(i)=(X(i)-X(j))*XX(i);
    end
end
for i=1:N
    for j=1:i-1
        A(i,j,1)=XX(i)/XX(j)/(X(i)-X(j));
    end
    for j=i+1:N
        A(i,j,1)=XX(i)/XX(j)/(X(i)-X(j));
    end
end
QM=-sum(A(:,:,1)');
for i=1:N
    A(i,i,1)=QM(i);
end
for r=2:4
    for i=1:N
        for j=1:i-1
            A(i,j,r)=r*(A(i,i,r-1)*A(i,j,1)-A(i,j,r-1)/(X(i)-X(j)));
        end
        for j=i+1:N            
            A(i,j,r)=r*(A(i,i,r-1)*A(i,j,1)-A(i,j,r-1)/(X(i)-X(j)));
        end
    end
   HM= -sum(A(:,:,r)');
   for i=1:N
       A(i,i,r)=HM(i);
   end
end
r1=0;r2=40;v=0.25
for i=1:N
    for k=1:N
        if i==2
        K(i,k)=A(i,k,2)+v/(X(i))*A(i,k,1);
        elseif i==N-1;
        K(i,k)=A(i,k,2)+v/(X(i))*A(i,k,1);
        else
          K(i,k)=A(i,k,4)+2/(X(i))*A(i,k,2)-1/(X(i)^2)*A(i,k,2)+1/(X(i)^3)*A(i,k,1);  
        end
    end
end
K1=-K(2:N-1,2:N-1);%B矩阵
K2=eye(2*N-2);
K=zeros(26,26);
for i=1:N-2
    for j=1:N-2
  K=K1;
    end;
end;
for i=N-1:(2*N-4)
    for j=N-1:(2*N-4)
        if i==j
  K(i,j)=1;
        end;
    end;
end;
DL=K1*W;
DW=K2*L;
Q=[W;L];
P=[L;W];
DP=K*Q;
end

[ Last edited by jjdg on 2017-4-17 at 22:33 ]
回复此楼

» 猜你喜欢

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

quarkhe

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by 约束陶瓷装甲 at 2017-04-17 19:35:39
既然是未知量,怎么能放到赋值号右边?...

新手,求具体怎么写?十分感谢!

发自小木虫Android客户端
6楼2017-04-18 05:10:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

quarkhe

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 约束陶瓷装甲 at 2017-04-17 11:23:08
W是啥?L是啥?

是两个未知待求解的列矩阵
3楼2017-04-17 18:02:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

约束陶瓷装甲

金虫 (正式写手)

★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
jjdg: 金币+1, 感谢参与 2017-04-17 22:34:11
引用回帖:
3楼: Originally posted by quarkhe at 2017-04-17 18:02:33
是两个未知待求解的列矩阵...

既然是未知量,怎么能放到赋值号右边?
4楼2017-04-17 19:35:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)

★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
jjdg: 金币+1, 感谢参与 2017-04-17 22:34:24
引用回帖:
3楼: Originally posted by quarkhe at 2017-04-17 18:02:33
是两个未知待求解的列矩阵...

可以定义为全局变量,在调用函数之前赋初值。

发自小木虫Android客户端
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
5楼2017-04-17 20:05:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +13 冬十三 2026-03-15 14/700 2026-03-21 22:10 by peike
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +10 枫桥ZL 2026-03-18 12/600 2026-03-21 22:02 by peike
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 12/600 2026-03-21 17:47 by ColorlessPI
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
信息提示
请填处理意见