24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yy3780

木虫 (著名写手)

应助: 0 (幼儿园)
金币: 4304.2
散金: 50
红花: 3
沙发: 3
帖子: 1079
在线: 401.9小时
虫号: 653605
注册: 2008-11-13
专业: 造血相关器官（肝脏/脾脏/

[交流] 【求助】MATLAB计算行列式【已解决】

我需要用MATLAB计算一个四阶的行列式，行列式的元素是与m有关的函数。
我MATLAB不熟，总是得不到想要的结果。希望大家帮帮我。谢谢大家了。我的QQ是763416717

[ Last edited by wangen994 on 2009-7-5 at 08:33 ]

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有5人回复
面上没中，邀请各位路过的虫友分析一下分数已经有7人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复
国社科又开始会评了，不知道这次命运如何已经有14人回复
麻烦专家们看看评委们的意见（F口面上）已经有13人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有14人回复
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够已经有17人回复
学科评审组评审是指会评吗？已经有5人回复
小白求助投论文要求的highlights应该如何写已经有4人回复

1楼 2009-07-04 12:24:10

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

hitzhang

木虫 (正式写手)

应助: 0 (幼儿园)
贵宾: 2.15
金币: 1356.7
散金: 969
红花: 8
帖子: 863
在线: 226.5小时
虫号: 390575
注册: 2007-06-02
性别: GG
专业: 无机非金属类电介质与电解

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
yy3780(金币+87,VIP+0):太好了，您真厉害，谢谢您的帮助，也谢谢大家 7-4 23:22
csfn(金币+10,VIP+0):辛苦了辛苦了，希望多来交流指导:-) 7-4 23:53
woshilsh(金币+0,VIP+0):鉴于张兄的热心和专业，邀请兄弟作为顾问，并在你的专贴里奖励Vip和金币，兄弟有意请Pm下我！ 7-4 23:54

引用回帖:

Originally posted by yy3780 at 2009-7-4 12:24:
我需要用MATLAB计算一个四阶的行列式，行列式的元素是与m有关的函数。
我MATLAB不熟，总是得不到想要的结果。希望大家帮帮我。谢谢大家了。我的QQ是763416717

[ Last edited by woshilsh on 2009-7-4 at 17: ...

///////////////////////////////////////////////////////////////
function code:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
function E=likethis(m,g,t)
s=size(g);
E=zeros(s(1),4);
for i=1:s(1)
E(i,=maybe(m,g(i),t);
end

function e=maybe(m,g,t)
A=zeros(4,4);
for i=1:4
for j=1:4
      A(i,j)=-t*D(m-1+i,m-1+j,g);
end
end
A=A+diag(m-g^2+[0:3]);
e=sort(eig(A))';

function d=D(m,n,g)
upk=min(m,n);
k=0:upk;
d=exp(-2*g^2)*sum((-1).^(-k).*(prod1(m).*prod1(n)).^0.5.*(2*g).^(m+n-2*k)./(prod1(m-k).*prod1(n-k).*prod1(k)));

function y=prod1(x)
s=size(x);
y=zeros(s);
for i=1:s(2)
if x(i)==0;
      y(i)=1;
else
      y(i)=prod(1:x(i));
end
end


\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
command code:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

>> g=0:0.01:2;
>> m=0;
>> t=0.5;
>> g=[0:0.01:2]';
>> E=likethis(m,g,t);
>> plot(g,E)
>> xlabel('g'),ylabel('E'),title('m=0, t=0.5')
>> legend({'E0','E1','E2','E3'})
>>

[ Last edited by hitzhang on 2009-7-4 at 23:13 ]