24小时热门版块排行榜    

查看: 588  |  回复: 0

烁烁88

银虫 (正式写手)

[求助] Matlab计算程序计算不出来,求助

n=10;
h=0.1;
A=zeros(n,n);
for i=2:n-2
    A(i,i)=6;A(i,i+2)=1;
    A(i,i-1)=-4;A(i+1,i-1)=1;A(i,i+1)=-4;
    i+1;
end
A(1,1)=7;A(2,2)=6;A(n-1,n)=-2;A(n-1,n-2)=-4;A(n-1,n-1)=5;
A(n,n)=1;A(1,2)=-4;A(n,n-1)=-2;A(n,n-2)=1;A(1,3)=1;
A
%上面就是构筑A矩阵,有方便方法还请指教
syms x1 x2 x3 x4 x5 x6 x7 x8 x9 x10;
u=x1;
for i=2:n
    syms(['x',num2str(i)]);
    u=[u,['x',num2str(i)]];
end
u=u';
%这又是构筑一个变量矩阵
[sigma, k, zita]=deal(1,100,1);
for i=1:n
    F(i)=sigma/(1-u(i))/(log(k*(1-u(i))))^2+zita/(1-u(i))^(5/2);
end
F=F';
f=A*u-F
%最终就是解这个方程,求出u矩阵内的每个变量。
solve(f)
回复此楼

» 猜你喜欢

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 烁烁88 的主题更新
信息提示
请填处理意见