| 查看: 563 | 回复: 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 程序报错，估计是要取正实根的问题，求怎么改。 另外如果直接拿原始方程编程该怎么写呢？求指导。 |
回复此楼» 猜你喜欢
一志愿沪9，生物学326求调剂
已经有3人回复
-
药学求调剂
已经有10人回复
-
307中医考研调剂
已经有3人回复
-
271求调剂
已经有31人回复
-
291求调剂
已经有4人回复
-
山东省基金2026
已经有11人回复
-
一志愿华中农业071010，320求调剂
已经有10人回复
-
材料工程281还有调剂机会吗
已经有43人回复
-
291求调剂
已经有5人回复
-
085801电气专硕272求调剂
已经有18人回复
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
