| 查看: 1221 | 回复: 4 | |||
shuxiao0811金虫 (小有名气)
|
[求助]
用MATLAB求解一个一元二次方程，满足一定的误差条件
|
|
本人MATLAB几乎小白，有个计算需要用MATLAB求解，试了两天还是没整出来，求帮忙补充一下程序，多谢！ 已知：数据文件ITO-R.xlsx和ITO-T.xlsx（文件见附件）。其中ITO-R.xlsx的第一列是波长Lamda，第二列是对应的反射率测量值R（测量时增益为100）；ITO-T.xlsx的第一列是波长Lamda，第二列是对应的反射率测量值T（测量时增益为100）； 试用MATLAB求解 1. K=-1*Lamda*ln（1-R-T）/（4*300）; 2. 由B=（（n-1）^2+K^2）/（（n+1）^2+K^2），求出n，要求abs（B-（（n-1）^2+K^2）/（（n+1）^2+K^2））<0.001； 3. 由K，n的值求 real=n^2-K^2, image=2nK 运行后输出K，n，real，image的值。 着急，在线等，多谢多谢~ 开头读取xlsx中数据的程序 clc clear all Num1=xlsread('GZO-R.xlsx'); %读取文件1数据 Num2=xlsread('GZO-T.xlsx'); %读取文件2数据 Lambda=Num1(:,1); %文件第一列值为Lambda R1=Num1(:,2); %文件1第二列值为反射 T1=Num2(:,2); %文件2第二列值为透射 R=R1/100; T=T1/100; d=300; K=-1/(4*d).*Lambda.*log(1-R-T); %计算K值 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : ITO-R.xlsx
- 附件 2 : ITO-T.xlsx
2016-08-26 20:05:28, 42.34 K
2016-08-26 20:05:34, 41.85 K
» 猜你喜欢
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
-
国自然结果
已经有6人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有9人回复
-
好奇怪的filecode
已经有5人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有21人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有10人回复
-
据悉今年马上要出结果了
已经有8人回复
-
关于代码变化问题，想知道的进来
已经有23人回复
-
FileCode能看出啥？
已经有6人回复
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
2楼2016-08-27 06:40:07
3楼2016-08-27 12:48:07
whqs8426212
铜虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 13.2
- 散金: 375
- 红花: 7
- 帖子: 782
- 在线: 685.7小时
- 虫号: 2572504
- 注册: 2013-07-30
- 专业: 数论
4楼2016-08-27 17:56:24
shuxiao0811
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 73.5
- 散金: 316
- 红花: 2
- 帖子: 161
- 在线: 67.9小时
- 虫号: 2396969
- 注册: 2013-04-02
- 性别: GG
- 专业: 光学
5楼2016-08-28 15:34:03
100