| 查看: 633 | 回复: 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 程序报错，估计是要取正实根的问题，求怎么改。 另外如果直接拿原始方程编程该怎么写呢？求指导。 |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有21人回复
-
面上项目filecode邪修
已经有3人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有13人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
xgp326265322
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 48.5
- 红花: 3
- 帖子: 427
- 在线: 130.2小时
- 虫号: 2409989
- 注册: 2013-04-09
- 专业: 化工热力学和基础数据
3楼2020-03-20 15:14:33
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
100