24小时热门版块排行榜    

查看: 2501  |  回复: 4

truewz

至尊木虫 (文坛精英)

[交流] 【求助】用MAPLE计算一个符号矩阵的行列式的解 已有2人参与

我对MAPLE不熟,想用MAPLE计算一个符号矩阵的行列式的解,如下:

                  x        2*a      b
D=      [    2-3b     a         x  ]   
                  4          5         x

%x是要求的量

L=det(D);  

solve(L);  %代入 a=1,b=2, 求出x的三个解


上面复合MATLAB 的语法,请高手用MAPLE编程实现上面的步骤,十分感谢!

[ Last edited by truewz on 2010-6-14 at 16:40 ]
回复此楼
扫地虫
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maplesim

★ ★
truewz(金币+5): 2010-06-12 18:04:43
robert2020(金币+2):多谢解答,欢迎常来仿真模拟版交流。 2010-06-12 18:20:54
两个方法可以得到结果:
1. 在Maple中实现

首先输入矩阵:(如果不知道格式可以看帮助,或者用左侧的矩阵面板)
CODE:
M := Matrix(3, 3, [[x, 2*a, b], [2-3*b, a, x], [4, 5, x]])

按回车键,然后用鼠标右键点击上面矩阵,从右键菜单中选择:
Standard Operations ->Determinant

Maple会自动显示求行列式的命令:
CODE:
LinearAlgebra:-Determinant(M)

然后再用鼠标右键,选择 solve下面的菜单,求出x的值。

如果熟悉Maple的命令,就不要这么费劲。

2. 可以直接在MATLAB下用Maple的命令实现。就是在MATLAB中直接输入上面Maple的命令,但输入格式稍微不同。
maple('PackageName[CommandName](arguments)');
对应上面的命令就是maple(LinearAlgebra[Determinant](M)');
2楼2010-06-12 16:42:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maplesim


truewz(金币+1): 2010-06-12 18:04:49
robert2020(金币+1):多谢解答! 2010-06-12 18:21:14
补充说明:

1. 在Maple中不能用D作为变量名,这是内部定义的微分算子。

2. 在Matlab中使用Maple函数,需要安装Maple Toolbox for Matlab,或者安装最新的Maple 14,这个工具箱免费集成到Maple中了。因为我们实验室用的正版Maple,所以在机房可以同时用这两个软件。
3楼2010-06-12 16:46:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

truewz

至尊木虫 (文坛精英)

robert2020:提示使用“引用回复该贴”。 2010-06-12 18:21:41
感谢楼上仁兄,我需要算的是一个循环程序,所以不能每次都点右键。
我的意思是把行列式赋给一个变量,在把符号变量代入,求出X值。
能否麻烦给出程序(再追加5个BB,共计25个),谢谢


注:我用MATLAB算的30*30的符号矩阵,算了几组数据后就算不下去了,所以想改用MAPLE

[ Last edited by truewz on 2010-6-12 at 18:06 ]
扫地虫
4楼2010-06-12 18:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maplesim


robert2020(金币+1):多谢耐心解答,欢迎常来仿真模拟版交流。 2010-06-13 10:51:36
truewz(金币+19): 2010-07-22 07:26:39
truewz(金币+1): 2010-11-26 17:24:20
你自己动手在Maple中试试。我上面计算出的行列式,本身就是x=f(a, b, c...)符号变量的表达式。如果代数这些符号变量的值,可以得到x的值。
赋值语句的格式是:a := 2:

如果需要简单的编程,也是比较容易。

Maple的语法与C语言非常相似,Maple内置的函数都是用C和C++写的。

这里举几个简单的例子(下面是从一本书上抄的):

Maple内置完整的程序语言,编写Maple程序实际上是非常简单的,只要在你每天使用的一系列命令前后分别加上proc()及end即可。利用Maple自带的函数程序构造应用程序,是一件极为容易的事情。

例子:
CODE:
prog := proc (x, y)
                    sqrt(x^2+y^2)
                   end proc

CODE:
prog(6, 4)

2*sqrt(13)


if..else 语句
例子:
CODE:
estart;
x := -12;
if `not`(type(x, integer)) then
printf("%a 不是一个整数。", x)
elif x >= 10 then
printf("%a 是一个1位以上的整数。", x)
elif x >= 0 then
printf("%a 是一个一位整数。", x)
else
printf("%a 是一个负整数。", x)
end if

for 循环

例子:
CODE:
for n to 5 do evalf(sqrt(n)) end do

5楼2010-06-13 09:13:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 truewz 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 310求调剂 +3 baibai1314 2026-03-16 3/150 2026-03-21 03:56 by JourneyLucky
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 324求调剂 +5 lucky呀呀呀鸭 2026-03-20 5/250 2026-03-20 22:30 by 促天成
[考研] 316求调剂 +5 梁茜雯 2026-03-19 5/250 2026-03-20 21:26 by 你好你好666
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 一志愿西安交通大学 学硕 354求调剂 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
信息提示
请填处理意见