| 查看: 471 | 回复: 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 程序报错，估计是要取正实根的问题，求怎么改。 另外如果直接拿原始方程编程该怎么写呢？求指导。 |
回复此楼» 猜你喜欢
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
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