24小时热门版块排行榜    

查看: 222  |  回复: 1

m镜花水月

金虫 (正式写手)

[求助] 跪求大神指点 已有1人参与

clc
clear all
close all
a=0.633*10.^(-6);
b=4000;
c=2.0*10.^(-6);%线性膨胀系数
c0=10*10.^(-3);%长度
d=2.2*10.^(-6);%电场方向系数
d0=5*10.^(-3);%电场方向的厚度
T0=298;%室温
T=233:0.1:333;%温度范围
n=sqrt(4.913-2.78*10.^(-2)*a.^2+(0.1173+1.65*10.^(-8)*T.^2)/(a.^2-(0.212+2.7*10.^(-8)*T.^2).^2));%寻常光折射率与温度的关系
G=(a/(2*b))*((d0+d*(T-T0))/(c0+c*(T-T0)))*(1/(n^3));%电光系数与温度的关系
figure('color',[1,1,1])
plot(T,G)
ylim([-0.1,0.1]);
title('T-G(T)温度与电光系数关系图')
xlabel('T/温度')
ylabel('G(T)/电光系数')
hold on
grid on
figure('color',[1,1,1])
plot(T,n)
ylim([1.4,1.65]);
title('T-n(T)温度与折射率关系图')
xlabel('T/温度')
ylabel('n(T)/电光系数')
hold on
grid on
我想做出“温度与电光系数关系图”和“温度与折射率关系图”但做出来不对,比如“温度与折射率关系图”应该是

但是跑出来却是

按照公式程序跑出来应该是曲线或者是有斜率的直线,不可能是水平线,另外电光系数的也不对,该怎样修改程序,跪求大神指点

跪求大神指点
理论上正确折射率.jpg


跪求大神指点-1
自己的错误折射率.jpg


跪求大神指点-2
自己的错误电光系数.jpg
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
m镜花水月: 金币+6, ★★★很有帮助 2017-03-30 19:36:23
jjdg: 金币+2, 感谢参与 2017-03-31 00:19:02
CODE:
clc
clear all
close all
a=0.633*10.^(-6);
b=4000;
c=2.0*10.^(-6);%线性膨胀系数
c0=10*10.^(-3);%长度
d=2.2*10.^(-6);%电场方向系数
d0=5*10.^(-3);%电场方向的厚度
T0=298;%室温
T=233:0.1:333;%温度范围
n=sqrt(4.913-2.78*10.^(-2)*a.^2+(0.1173+1.65*10.^(-8)*T.^2)./(a.^2-(0.212+2.7*10.^(-8)*T.^2).^2));%寻常光折射率与温度的关系
G=(a/(2*b))   *   (     (  d0+d*(T-T0) )./ (   c0+c*(T-T0))     ).*(1./(n.^3)   );%电光系数与温度的关系
figure('color',[1,1,1])
plot(T,G)
ylim([0.11,0.112]*1e-10);
title('T-G(T)温度与电光系数关系图')
xlabel('T/温度')
ylabel('G(T)/电光系数')

grid on
figure('color',[1,1,1])
plot(T,n)
ylim([1.522,1.529]);
title('T-n(T)温度与折射率关系图')
xlabel('T/温度')
ylabel('n(T)/电光系数')

grid on

跪求大神指点-3
1.png


跪求大神指点-4
2.png

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2017-03-30 11:55:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 m镜花水月 的主题更新
信息提示
请填处理意见