如何在matlab中绘制t分布曲线
发布网友
发布时间：2022-04-20 08:08
共3个回答
热心网友
时间：2023-05-18 10:40
试试看这个
clear all;
clc;
L1=1; L2=1; L=1; C1=1; C2=1; C=1; m=0.5; w=50;
a=sqrt(2)/2; b=sqrt(2)/2; d=sqrt(2)/2; c=-sqrt(2)/2;
k=1.381*10^(-23); h=6.63*10^(-34);
T=273;
B=1/(k*T);
M1=1/(L-m); M2=1/(L+m);
r1=1/sqrt(C*(L-m)); r2=1/sqrt(C*(L+m));
a1=0; b1=0;
t=0:0.1:10;
a2=-sqrt(2)*(cos((w+r2)*t)/(2*(w+r2))+cos((w-r2)*t)/(2*(w-r2))-1/(2*(w+r2))-1/(2*(w-r2)));
b2=-(sqrt(2)*M2/r2)*(sin((w+r2)*t)/(2*(w+r2))-sin((w-r2)*t)/(2*(w-r2)));
A1=a1*cos(r1*t)+b1*r1*sin(r1*t)/M1;
B1=b1*cos(r1*t)-a1*M1*sin(r1*t)/r1;
A2=a2.*cos(r2*t)+b2*r2.*sin(r2*t)/M2;
B2=b2.*cos(r2*t)-a2*M2.*sin(r2*t)/r2;
q1=a*B1+b*B2; q2=c*B1+d*B2;
I1=(m*(a*A2-b*A1)+L2*(c*A2-d*A1))/(m^2-L1*L2);
I2=(m*(d*A1-c*A2)-L1*(a*A2-b*A1))/(m^2-L1*L2);
k1=(a^2)*h*M1/(2*r1*tanh(h*r1*B/2))+(b^2)*h*M2/(2*r2*tanh(h*r2*B/2))+2*a*b.*B1.*B2;
k2=(c^2)*h*M1/(2*r1*tanh(h*r1*B/2))+(d^2)*h*M2/(2*r2*tanh(h*r2*B/2))+2*c*d.*B1.*B2;
s1=k1-(q1).^2;
s2=k2-(q2).^2;
subplot(2,1,1);plot(t,I1);hold on; grid on; title('温度为273K时电流I1的热力学平均值随时间t变化的曲线');xlabel('t/时间');ylabel('I1');
subplot(2,1,2);plot(t,I2);hold on; grid on; title('温度为273K时电流I2的热力学平均值随时间t变化的曲线');xlabel('t/时间');ylabel('I2');
热心网友
时间：2023-05-18 10:41
ezplot('pdf(''t'',t,24)',[-3 3])
ax=gca;
td=[-2:0.01:1];
yd=pdf('t',td,24);
td=[-2 td 1];
yd=[0 yd 0];
p=patch('Parent',ax,'YData',yd,'XData',td,'FaceColor',[0 0 1]);
ax.YLim=[0 0.4];
热心网友
时间：2023-05-18 10:41
数据见Excel文件“成绩”工作表。请编制数学成绩（math）分布数列，并画出直方图。
分组为：50-60，60-70，70-80，80-90，90-100
x=load('1.txt')
y=[50:10:100]
histcounts(x,y)
z=histogram(x,y)
数据见Excel文件“成绩”工作表。请同时按班级（class）和性别（gender）进行分类，求delphi和english的平均分，标准差，最高分，最低分。
x=readtable('2.txt')
t=x(:,[2 3 5 8])
grpstats(t,{'class','gender'},{'mean','std','max','min'})
数据见Excel文件“耐用时数”工作表。请计算x的平均数与标准差。
x=readtable('3.txt')
data=table2array(x)
x=data(:,1)
f=data(:,2)
xbar=x'*f/sum(f)
std=sqrt((x-xbar).^2'*f/sum(f))
使用加权公式计算平均收入与标准差
data=[10000 30000 50000 70000 90000 110000;
3 11 36 26 14 10]'
x=data(:,1)
f=data(:,2)
xbar=x'*f/sum(f)
std=sqrt((x-xbar).^2'*f/sum(f))
数据见Excel文件“设备价值”工作表。请画出类似的柱形图。
f=fopen('5.txt')
title1=textscan(f,'%s %s',1)
y=textscan(f,'%s %f')
type=y{1}
data=y{2}
bar(data)
ax=gca
ax.XLim=[0.5 5.5]
ax.XTickLabel=type
ax.XLabel.String='项目'
ax.YTick=[0:10000:140000]
ax.YTickLabel=[0:10000:140000]
ax.YLabel.String='金额'
title(ax,'设备金额')
hold on
scatter([1 2 3 4 5],data','b')
plot(data,'r')
画出y=x2 x2+y2=9
ezplot('x^2+y^2-9',[-3,3])
axis equal
hold on
ezplot('x^2',[-3,3,-3,3])
text(1,-2,'x^2+y^2=9')
text(1.1,1,'y=x^2')
title('')
使用描点法绘制t分布密度函数图，自由度取24。
并求出t变量在区间[-2,1]上的概率，并作出如下示意图
x=[-3:0.01:3]
y=tpdf(x,24)
plot(x,y)
y1=tpdf(-2,24)
y2=tpdf(1,24)
hold on
x0=[-2:0.1:1]
y0=tpdf(x0,24)
fill([-2 x0 1],[0 y0 0],'b')
ax=gca
ax.XLabel.String='t'
title("pdf('t',t,24)")
t=tcdf(1,24)-tcdf(-2,24)
某航空集团1949～1960年的国际航线季度旅客数见Excel文件“旅客数量”工作表。试对其建立数曲线趋势模型。Y为旅客数量，t为时间。
y=load('8.txt')
t=[1:1:48]'
lny=log(y)
data=table(t,lny)
%m1=fitlm(t,lny)
m2=fitlm(data,'lny ~ 1+t')
数据见Excel文件“成绩2”工作表。
1.读取数据
将数据另存为txt文件，注意不要第1行。在Matlab中读取该txt数据到一个矩阵中。
2.然后将type=1的记录筛选出来，放到另外一个矩阵中。
3.对原矩阵zy1到zy4列，全部改为0
4.对原矩阵删除zy1列
x=readtable('9.txt','delimiter','\t','ReadVariableNames',0)
t=table2array(x)
y1=t(t(:,3)==1,:)
t(:,4:7)=0
t(:,4)=[]
