| 查看: 1188 | 回复: 10 | |||
| 当前主题已经存档。 | |||
yy3780木虫 (著名写手)
|
[交流]
【求助】MATLAB计算行列式【已解决】
|
||
|
我需要用MATLAB计算一个四阶的行列式，行列式的元素是与m有关的函数。 我MATLAB不熟，总是得不到想要的结果。希望大家帮帮我。谢谢大家了。我的QQ是763416717 [ Last edited by wangen994 on 2009-7-5 at 08:33 ] |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
2楼2009-07-04 12:49:08
yy3780
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 4304.2
- 散金: 50
- 红花: 3
- 沙发: 3
- 帖子: 1079
- 在线: 401.9小时
- 虫号: 653605
- 注册: 2008-11-13
- 专业: 造血相关器官（肝脏/脾脏/
3楼2009-07-04 13:07:45
dfq0730
木虫 (文坛精英)
顺天府府尹
- 应助: 0 (幼儿园)
- 贵宾: 7.89
- 金币: 1752.8
- 散金: 5832
- 红花: 18
- 沙发: 16
- 帖子: 14630
- 在线: 191.4小时
- 虫号: 613132
- 注册: 2008-09-26
- 性别: GG
- 专业: 高分子物理与高分子物理化
4楼2009-07-04 17:34:12
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
5楼2009-07-04 17:53:08
yy3780
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 4304.2
- 散金: 50
- 红花: 3
- 沙发: 3
- 帖子: 1079
- 在线: 401.9小时
- 虫号: 653605
- 注册: 2008-11-13
- 专业: 造血相关器官（肝脏/脾脏/
6楼2009-07-04 19:06:47
hitzhang
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
7楼2009-07-04 20:36:36
yy3780
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 4304.2
- 散金: 50
- 红花: 3
- 沙发: 3
- 帖子: 1079
- 在线: 401.9小时
- 虫号: 653605
- 注册: 2008-11-13
- 专业: 造血相关器官（肝脏/脾脏/
8楼2009-07-04 21:27:20
9楼2009-07-04 21:56:12
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
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
|
/////////////////////////////////////////////////////////////// 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 ] |
10楼2009-07-04 23:11:21
200