应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Matlab,解复数方程
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 2074  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小刘8885

木虫 (正式写手)

[求助] Matlab,解复数方程

怎么解这个方程？如下图。其中x1和x2都为复数，q和a为已知的常数。x1和x2为数据文件，如下图所示。怎么导入x1和x2并解出一系列对应的y值，并且将求的的y值以类似x1的文件形式导出来，即分为实部和虚部导出？
公式



数据
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : shuju.txt
    • 2012-05-21 20:01:14, 2.77 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-21 20:01:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★
感谢参与，应助指数 +1
小刘8885: 金币+1, 谢谢参与 2012-05-21 21:10:47
用1stOpt试试，比较简单，下面假设q=3.5, a=0.25，以前10个数据为例计算：
CODE:
Constant q=3.5,a=0.25;
ConstStr x1=a1+b1*i,x2=a2+b2*i;
LoopConstant
a1=[1.43,1.44,1.46,1.45,1.45,1.44,1.45,1.45,1.44,1.44],
b1=[0.02,0.04,0.05,0.02,0.05,0.06,0.06,0.05,0.07,0.09],
a2=[0.009780429,0.019290123,0.023456558,0.009512485,0.023781213,0.028935185,0.028537455,0.023781213,0.033757716,0.043402778],
b2=[0.004890215,0.008768238,0.009773566,0.003658648,0.00849329,0.009645062,0.008917955,0.006994474,0.009377143,0.011421784];
ComplexStr = i;
ComplexPar y;
Function (1-q)*(x2-x1)/x2^a*y^a-y+x1=0;

y.实部        y.虚部
26.0549308935745        -3.39054421631229
21.3634333717575        -2.15713840435354
20.5571962148911        -1.69871196786957
27.2314765798707        -2.70147269442524
20.4345049478304        -1.33862681249657
19.0671427961527        -0.949477270203328
19.371747766277        -0.850707043007433
20.575089935508        -0.966582317014032
18.2424521984911        -0.447405238620689
16.8276617636246        -0.0454317407652095
赞一下 回复此楼
2楼2012-05-21 20:54:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie169

铜虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
小刘8885: 金币+8, ★★★很有帮助 2012-05-22 08:50:14
先列上matlab代码
CODE:
function f=main(name,a,q)
data=load(name);
x1=data(:,1)+data(:,2)*i;
x2=data(:,3)+data(:,4)*i;
x0=zeros(length(data),1);
f=fsolve(@(y)func(a,q,x1,x2,y),x0);
r=real(f);
r(:,2)=imag(f);
save('result.txt','r','-ascii');
function f=func(a,q,x1,x2,y)
f=(1-q).*(x1-x2)./x2^a.*y^a-y+x1;
end
end

首先，在matlab下面新建一个main.m的文件，复制上我写的代码。然后把数据文件拷贝到main.m的目录下。
name就是你的数据文件，里面除了列之间的数据，不能有其他的东西，数据文件的名字英文。a，q不用我多说
调用格式是：main('data.txt',1,5)
这是一个示例data.txt就像相当于你的文件的名字。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

小刘8885
3楼2012-05-22 07:09:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小刘8885

木虫 (正式写手)
送鲜花一朵
引用回帖:
3楼: Originally posted by lijie169 at 2012-05-22 07:09:32:
先列上matlab代码

function f=main(name,a,q)
data=load(name);
x1=data(:,1)+data(:,2)*i;
x2=data(:,3)+data(:,4)*i;
x0=zeros(length(data),1);
f=fsolve(@(y)func(a,q,x1,x2,y),x0);
r=real(f);
r(:,

谢谢！我将你的代码写入之后运行怎么出错了？这是什么原因呢？求解。
>> main('2.txt',1/3,0.96)
??? Error: File: D:\Matlab\work\main.m Line: 6 Column: 11
"identifier" expected, "(" found.
赞一下 回复此楼
4楼2012-05-22 08:49:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小刘8885

木虫 (正式写手)
引用回帖:
3楼: Originally posted by lijie169 at 2012-05-22 07:09:32:
先列上matlab代码

function f=main(name,a,q)
data=load(name);
x1=data(:,1)+data(:,2)*i;
x2=data(:,3)+data(:,4)*i;
x0=zeros(length(data),1);
f=fsolve(@(y)func(a,q,x1,x2,y),x0);
r=real(f);
r(:,

还有就是能不能对每行的程序做个简单的注释？我是初学者，很多都不太明白，求高人赐教。感激不尽。
赞一下 回复此楼
5楼2012-05-22 08:52:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie169

铜虫 (著名写手)
建议使用2楼的方案，matlab解这个需要初始解，初始解给的不好，木有答案！！！！
赞一下 回复此楼
6楼2012-05-22 09:36:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xcj3630385

木虫 (小有名气)

长江11号
function f= f_main(data, a, q)
x1 = data(:,1) + data(:,2)*i;
x2 = data(:,3) + data(:,4)*i;

x0 = zeros(length(data),1);
f = fsolve(@(y)func(a,q,x1,x2,y), x0)
r = [real(f), imag(f)];
save('result.txt','r','-ascii');

function f = func(a, q, x1, x2, y)
f = (1-q).*(x2-x1)./(x2.^a).*(y.^a) - y + x1;
%变形式。
赞一下 回复此楼
长江11号
7楼2012-05-25 12:33:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xcj3630385

木虫 (小有名气)

长江11号
引用回帖:
7楼: Originally posted by xcj3630385 at 2012-05-25 12:33:54
function f= f_main(data, a, q)
x1 = data(:,1) + data(:,2)*i;
x2 = data(:,3) + data(:,4)*i;

x0 = zeros(length(data),1);
f = fsolve(@(y)func(a,q,x1,x2,y), x0)
r = ;
save('result.txt','r','-asci

调用方式如下：

clear all;
clc;

data = load('name.txt');  %name即调用文件的名字
a = 0.2;   %任意值
q = 0.06;  %任意值
f_main(data, a, q);
赞一下 回复此楼
长江11号
8楼2012-05-25 12:37:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小刘8885

木虫 (正式写手)
引用回帖:
8楼: Originally posted by xcj3630385 at 2012-05-25 12:37:44
调用方式如下：

clear all;
clc;

data = load('name.txt');  %name即调用文件的名字
a = 0.2;   %任意值
q = 0.06;  %任意值
f_main(data, a, q);

9楼2012-05-25 12:41:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小刘8885 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307求调剂 +7 超级伊昂大王 2026-04-06 7/350 2026-04-06 17:12 by lijunpoly
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-05 5/250 2026-04-06 15:40 by lin-da
[考研] 308求调剂 +13 倘若起风了呢 2026-04-05 13/650 2026-04-06 14:20 by 蒋皓禹
[考研] 304求调剂 +8 c297914 2026-04-05 9/450 2026-04-06 14:18 by 蒋皓禹
[考研] 0854求调剂 +4 assdll 2026-04-05 4/200 2026-04-06 12:29 by 中飞院空管学院�
[考研] 求调剂 +5 wos666 2026-04-03 5/250 2026-04-06 10:13 by 蓝云思雨
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 292分，材料与化工，申请调剂 +22 程晴之 2026-04-01 26/1300 2026-04-04 22:03 by hemengdong
[考研] 272求调剂 +4 松柏常青5 2026-04-03 4/200 2026-04-04 17:03 by babysonlkd
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +3 蓝笺片 2026-04-03 4/200 2026-04-03 17:58 by Jimmyandyou
[考研] 313求调剂 +3 ～微微凉～ 2026-04-03 3/150 2026-04-03 11:25 by 啵啵啵0119
[考研] 材料340分调剂 +7 夏夜晚风_long 2026-04-02 9/450 2026-04-02 21:20 by dongzh2009
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 348环境工程调剂 +3 吴彦祖24k 2026-04-01 3/150 2026-04-02 09:14 by nanaliuyun
[考研] 296求调剂 +4 汪！？！ 2026-03-31 7/350 2026-04-01 22:04 by 客尔美德
[考研] 求调剂0703 +5 周嘉尧 2026-03-31 8/400 2026-04-01 20:32 by ltltkkk
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 求调剂 +4 DADA怪 2026-03-31 4/200 2026-04-01 14:30 by ZXlzxl0425
信息提示
关闭
请填处理意见
关闭