| 查看: 636 | 回复: 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 程序报错，估计是要取正实根的问题，求怎么改。 另外如果直接拿原始方程编程该怎么写呢？求指导。 |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有42人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有29人回复
-
关于Filecode分析方法
已经有12人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有11人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
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
