24小时热门版块排行榜    

查看: 921  |  回复: 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的回帖

约束陶瓷装甲

金虫 (正式写手)

W是啥?L是啥?
2楼2017-04-17 11:23:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

quarkhe

新虫 (初入文坛)

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

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

发自小木虫Android客户端
6楼2017-04-18 05:10:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 quarkhe 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 298求调剂 +4 上岸6666@ 2026-03-20 4/200 2026-03-21 17:14 by 学员8dgXkO
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
信息提示
请填处理意见