hrefspace

 找回密码
 立即注册
搜索
热搜: PHP PS 程序设计
查看: 789|回复: 5

Mathematica 代码求助

[复制链接]

585

主题

769

帖子

2007

积分

大司空

Rank: 5Rank: 5

积分
2007
发表于 2023-10-3 08:02:21 | 显示全部楼层 |阅读模式
求助高人,如下的代码:

ClearALl["Global`*"];
a = x I; a' = -a; b = b' = 0; c = c' = x; d = x (1 + I);
d' = x (1 - I); e = x + y I; e' = x - y I;
Duichendian[p_, a_, b_] := (a' b - a b' + p' (a - b))/(a' - b');
Duichendian'[p_, a_, b_] := (a b' - a' b + p (a' - b'))/(a - b);
k[a_, b_] := (a - b)/(a' - b');
g = Duichendian[d, a, e]; g' = Duichendian'[d, a, e];
sol = Solve[{k[g, b] == k[g, f], k[c, f] == -1}, {f, f'}];
Print["f=", f /. Last[sol], ", f'=", f' /. Last[sol]];

为什么第一次运行显示:
f=(((1-I) ((1+I) x^3 y-x^2 y^2-(1+I) x^2 y (0&)+x y^2 (0&)+(1+I) x^2 y x^\[Prime]-x y^2 x^\[Prime]))/(4 x^2 y-2 x y^2-2 x^2 (0&)+2 x y (0&)-y^2 (0&))), f'=-((-2 x^3 y+(1+I) x^2 y^2+2 x^3 (0&)-4 x^2 y (0&)+(2-I) x y^2 (0&)-2 x^2 y x^\[Prime]+(1+I) x y^2 x^\[Prime]+2 x^2 (0&) x^\[Prime]-2 x y (0&) x^\[Prime]+y^2 (0&) x^\[Prime])/(4 x^2 y-2 x y^2-2 x^2 (0&)+2 x y (0&)-y^2 (0&)))

再运行一次则显示:
f=((1-I) x ((1+I) x-y))/(2 x-y), f'=-(((1+I) x ((-1+I) x+y))/(2 x-y))

Mathematica 版本:12.3
平台:win10


本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
回复

使用道具 举报

0

主题

189

帖子

163

积分

关内侯

Rank: 2

积分
163
发表于 2023-10-3 08:03:03 | 显示全部楼层
换成上划线也是一样的:

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
回复

使用道具 举报

0

主题

182

帖子

67

积分

关内侯

Rank: 2

积分
67
发表于 2023-10-3 08:03:27 | 显示全部楼层
经查,将 b = b' = 0; c = c' = x; 换成 b' = b = 0; c' = c = x; 就好使了。

即改成:
ClearAll["Global`*"];
a = x I; a' = -a; b' = b = 0; c' = c = x; d = x (1 + I);
d' = x (1 - I); e = x + y I; e' = x - y I;
Duichendian[p_, a_, b_] := (a' b - a b' + p' (a - b))/(a' - b');
Duichendian'[p_, a_, b_] := (a b' - a' b + p (a' - b'))/(a - b);
k[a_, b_] := (a - b)/(a' - b');
g = Duichendian[d, a, e]; g' = Duichendian'[d, a, e];
sol = Solve[{k[g, b] == k[g, f], k[c, f] == -1}, {f, f'}];
Print["f=", f /. Last[sol], ", f'=", f' /. Last[sol]];

那位达人能告诉我是怎么回事吗?
回复

使用道具 举报

0

主题

200

帖子

2

积分

新手上路

Rank: 1

积分
2
发表于 2023-10-3 08:04:13 | 显示全部楼层
我这里是这样的:

1,第一次运行,与你的"第一次运行显示"相同

2,对(1)再运行一次,与你的"再运行一次则显示"相同

3,对“第一次运行显示”作“化简”,也可以有“再运行一次则显示”

也有可能是需要对“文件”作“退出”“保存”“重启”。

我不行,想起网友 Mathematica是高手,可惜跑了。
回复

使用道具 举报

0

主题

192

帖子

2

积分

新手上路

Rank: 1

积分
2
发表于 2023-10-3 08:04:59 | 显示全部楼层
不要使用附标(下标、上横线、撇)定义变量。
因为那些东西根本不是变量,而是表达式。
你会遇到各种各样神奇的bug,没人想解决这些问题。
回复

使用道具 举报

0

主题

190

帖子

18

积分

新手上路

Rank: 1

积分
18
发表于 2023-10-3 08:05:33 | 显示全部楼层
基础语法问题 一大堆,不止四处。我就不一一罗列了。
感觉最好还是 先快速走一遍 速成教程
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|Archiver|手机版|小黑屋|hrefspace

GMT+8, 2024-11-22 09:59 , Processed in 0.079122 second(s), 22 queries .

Powered by hrefspace X3.4 Licensed

Copyright © 2022, hrefspace.

快速回复 返回顶部 返回列表