| 查看: 673 | 回复: 2 | |||
bustertang新虫 (初入文坛)
|
[求助]
matlab二元方程组求解问题
|
|
已知ab两个数组所有元素，ab的元素个数都为146个，a=n^2-k^2,b=2nk 求对应的nk数组，自己手算下发现n^2-b^2/4n^2-a=0,然后写了 n=[]; for i=1:1:146 n(i)=solve('x^2-b(i)^2/4*x^-2-a(i)','x') end 程序报错，估计是要取正实根的问题，求怎么改。 另外如果直接拿原始方程编程该怎么写呢？求指导。 |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有15人回复
-
申请科研项目纯靠实力能成功吗？
已经有8人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
bustertang
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 20.2
- 红花: 1
- 帖子: 19
- 在线: 13.3小时
- 虫号: 3633651
- 注册: 2015-01-08
- 性别: GG
- 专业: 金属结构材料
|
似乎自己搞定了。。。把程序放下面 photonenergyev1=photonenergyev1(~isnan(photonenergyev1)); photonenergyev2=photonenergyev2(~isnan(photonenergyev2)); epsilon1=epsilon1(~isnan(epsilon1)); epsilon2=epsilon2(~isnan(epsilon2)); photonenergyevfit=1.55:0.01:3; epsilon1fit=interp1(photonenergyev1,epsilon1,photonenergyevfit,'spline'); epsilon2fit=interp1(photonenergyev2,epsilon2,photonenergyevfit,'spline'); n=[]; for i=1:1:146 p=[4 0 -4*epsilon1fit(i) -epsilon2fit(i)]; x=roots(p); for ii=1:length(x) if x(ii)>0 & isreal(x(ii)); n(i)=x(ii); k(i)=epsilon2fit(i)/n(i); end end end output=[photonenergyevfit',n',k',epsilon1fit',epsilon2fit']; |
2楼2016-02-10 12:28:47
xgp326265322
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 48.5
- 红花: 3
- 帖子: 427
- 在线: 130.2小时
- 虫号: 2409989
- 注册: 2013-04-09
- 专业: 化工热力学和基础数据
3楼2020-03-20 15:14:33
10