2025认证杯数学建模C题成品论文.docx

PAGE1

PAGE1

C

2.1

2.1

2.2

4.1

4.1

4.2

1

5.1

5.1

5.2

5.2.1

5.2.2

5.3

5.4

5.4

5.5

5.6

6.16.2

6.1

6.2

6.2.1

y(t)

5

(t)0 i,xi(t)

(t

)

i1

6.2.2

6.36.3.16.3.2

6.3

6.3.1

6.3.2

6.3.3

10

10

6.46.4.16.4.2

6.4

6.4.1

6.4.2

6.4.3

11

11

6.5

6.5.1

6.5.2

6.5.2

1PAGE2

1PAGE2

6.6

7.1

7.2

1

7.2.1

7.2.2

7.2.3

7.2.3

7.3

7.3.1

7.3.27.47.4.1

7.3.2

7.4

7.4.1

7.4.2

7.5

7.5

1PAGE6

1PAGE6

7.5.1

7.5.2

7.5.2

7.5.3

8.1

8.1

1PAGE7

1PAGE7

8.2

8.2

8.2.1

8.2.2

1PAGE8

1PAGE8

8.2.3

1PAGE9

1PAGE9

8.3

8.3

8.3.1

8.3.2

20

20

8.4

8.5t 70

8.5

t 70

*

2PAGE1

2PAGE1

9.1

9.1.1

9.1.29.1.39.29.2.1

9.1.2

9.1.3

9.2

9.2.1

9.2.2

9.2.3

10.110.2

10.1

10.2

10.3

[1]

[2]

[3]

[4]

[5]

[6]

[6]

importpandasaspd

importmatplotlib.pyplotaspltimportnumpyasnp

#

df=pd.read_csv(OUT_Table.csv)so2_data=df[SO2]

h2s_data=df[H2S]

# 0

so2_corrected=so2_data.clip(lower=0)h2s_corrected=h2s_data.clip(lower=0)

)# 2x2

)

fig,axes=plt.subplots(2,2,figsize=(12,12)

0.7den-# SO

0.7

den-

lw1.2axes[0,0].hist(so2_data,bins=30,alpha= ,color=#1f77b4,sity=True)

lw

1.2

axes[0,0].axvline(0,color=red,linestyle=--, = )axes[0,0].set_title(r(a)$\mathrm{SO_2}$ ,fontsize=12)axes[0,0].set_xlabel( ,fontsize=10)

10axes[0,0].set_ylabel( ,fontsize=)

10

30# SO

30

].axvline(redaxes[0,1].hist(so2_corrected,bins=density=True)

].axvline(

red

0,alpha=0.7,color=#2ca02c,

axes[0,1 0,color= ,linestyle=--,lw=1.2)

axes[0,1].set_title(r(b)$\mathrm{SO_2}$ ,fontsize=12)

axes[0,1].set_xlabel( ,fontsize=10)

H# S

H

axes[1,0].hist(h2s_data,bins=30,alpha=0.7,color=#1f77b4,den-sity=True)

axes[1,0].axvline(0,color=red,linestyle=--,lw=1.2)axes[1,0].set_title(r(c)$\mathrm{H_2S}$ ,fontsize=12)axes[1,0].set_xlabel( ,fontsize=10)

axes[1,0].set_ylabel( ,fontsize=10)

# HS

axes[1,1].hist(h2s_corrected,bins=30,alpha=0.7,color=#2ca02c,density=True)

axes[1,1].axvline(0,